





THE JOURNAL 


OF 


CHEMICAL PHYSICS 





VoLuME 9 


SEPTEMBER, 1941 


NuMBER 9 





The X-Ray Diffraction Pattern of Tricalcium Aluminate 


Lynn J. BRADY AND WHEELER P. DAVEY 
The Pennsylvania State College, State College, Pennsylvania 
(Received May 13, 1941) 


Using a new type of x-ray diffraction apparatus, the diffraction pattern of high purity 
tricalcium aluminate, (3CaO-Al.O3;), has beer. recorded automatically in the form of a graph of 
intensity vs. angle. The unit cube containing 264 atoms has an edge of 15.235A. 


EPEATED attempts have been made to de- 

termine the crystal structure of tricalcium 
aluminate (3CaO-Al,O3) (1-2-3-4—5). Its com- 
plete x-ray diffraction pattern is, therefore, of 
interest. It is known to be optically isotropic, so 
its crystals must have cubic symmetry. No 
method is known by which crystals can be made 
large enough to permit goniometric measure- 
ments. It is the purpose of this paper to give 
the diffraction pattern of tricalcium aluminate 
in a more complete form than seems to be given 
so far in the literature. This is necessary before a 
solution of its structure can be found. 

Although previous workers have agreed on the 
position of the strongest lines in the diffraction 
pattern, no two workers have agreed on the in- 
tensity, the position, or the number of the 
weaker lines, nor have they agreed on the 
intensities of the stronger lines. 

Harrington! found a number of lines in his 
powder photographs which could not be ac- 
counted for on the basis of a cubic crystal having 
an edge of 7.62A. He assumed, therefore, that 
this compound was only pseudo-cubic. Examina- 
tion of his published data shows, however, that 
all his lines fit a cube whose edge is 15.24A. 





'E. A. Harrington, Am. J. Sci. 13, 467 (1927). 
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Steele and Davey? found, in addition to the 
lines characteristic of a body-centered cube, 
three faint lines in their photographs which they 
were able to account for qualitatively on the 
basis of a unit cube of edge 7.62A. On this basis 
they proposed a structure which obviously cannot 
be correct if the edge of the unit cube is 15.24A. 

Kondo and Yamauchi* and Solacolu‘ while 
proposing no structure agreed that the length of 
the edge of the unit crystal was approximately 
15.2A. 

Lagerqvist, Wallmark, and Westgren® found so 
many lines in their photographs that they had 
to assume the length of the unit cell to be 
15.22A in order to account for all of them. They 
found their lines to match portions of a simple 
cubic pattern, but they have not, as yet, pro- 
posed a structure. 

It is evident from the above that, if we are to 
arrive at a tenable solution of the structure of 
3CaO-Al,03, specimens must be of high enough 


2F. A. Steele and W. P. Davey, J. Am. Chem. Soc. 51, 
2283 (1929). 

3S. Kondo and T. Yamauchi, J. Japan. Ceram. Assoc. 
42, 313 (1934). 

4S. Solacolu, J. Japan. Ceram. Assoc. 42, 319 (1934). 

5 K. Lagerqvist, S. Wallmark, and A. Westgren, Zeits f. 
anorg. allgem. Chemie 234, 1 (1937). 
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Fic. 1. Recorder circuit. 


purity to justify precise determination of the 
diffraction patterns. 


PREPARATION OF SAMPLES 


Early in the course of this work attempts were 
made to make tricalcium aluminate using the 
technique of Lagerqvist, Wallmark, and West- 
gren. These workers fused 75 mole percent CaO 
and 25 mole percent Al,O; together, using the 
oxyhydrogen blowpipe. The fused mass was then 
ground and heated at 1000°C for a period of ten 
days. Although their published’ procedure was 
followed carefully, the petrographic microscope 
showed that anisotropic material was present at 
the end of the twelfth day of heating at 1000°C. 

Shepherd, Rankin, and Wright® showed that 
when 3CaQO-Al.O; is formed from the melt, 
crystals of the tricalcium aluminate frequently 
surround particles of CaO, with the consequent 
local formation of 5CaO-3Al1,03. They further 
showed that a period of 21 days of heating at 
1400°C is necessary to secure complete elimina- 
tion of CaO and 5CaO-3Al.03. 

In the present work, the starting materials 
for the formation of all the calcium aluminates 
were prepared by separately igniting precipitated 
chalk and ‘‘Bakers”’ analyzed aluminum oxide in 
platinum dishes at 1000°C for six hours. The 
resulting CaO and Al,O3; were ground and 
thoroughly mixed in the ratios of 3-1, 5-3, 1-1, 
and 3-5. These mixtures were fused at least once 
and were then ground to fine powders. Each of 
these specimens was subjected to several cycles 


6 E. S. Shepherd, G. A. Rankin, and F. E. Wright, Am. 
J. Sci. 28, 293 (1909). 
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of sintering and grinding, as in the work of 
Steele and Davey. A final purification was made 
on the basis of density, using methylene iodide 
and mesitilene. The resulting four purified prod- 
ucts were examined for homogeneity under a 
petrographic microscope and their refractive 
indices determined. The identities of these speci- 
mens were considered to be established when 
they showed the accepted densities and refractive 
indices. 


APPARATUS 


The diffraction pattern was measured by 
means of a Geiger-Miiller counter. The Geiger- 
Miiller circuits such as Gingrich, Evans, and 
Edgerton’ have developed are well adapted to 
x-ray measurements where the x-ray tube is 
operated under constant, or nearly constant, 
conditions. However, to make a circuit uni- 
versally adaptable it should be made _ inde- 
pendent of fluctuations in the source due to 
changing line voltages on the x-ray filament 
transformer. Such a modified circuit, completely 
a.c. operated, was constructed which records 
the average intensities of the x-ray beams with a 
recording galvanometer.® 

In this circuit pulses from the Geiger-Miiller 
counter are amplified by two 6J7-G tubes in 
series to secure larger, more nearly uniform, 
pulses. These pulses are then converted to a 
uniform size by means of the relaxation oscillator 
type circuit. Pulses from the uniform pulse 
generator are amplified and fed into a so-called 
‘“‘tank”’ condenser having a fixed high resistance 
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Fic. 2. Power circuits. 





7™N. S. Gingrich, R. D. Evans, and H. E. Edgerton, 
Rev. Sci. Inst. 7, 451 (1936). 
8L. J. Brady, Phys. Rev. 55, 1133 (1939). 





on, 








TRICALCIUM ALUMINATE 665 


leak. The potential of the condenser is thus a 
measure of the average rate of arrival of x-ray 
quanta at the Geiger-Miiller counter. This poten- 
tial is then recorded by a recording galvanometer 
in terms of the current leaking off the condenser 
through the resistance. 

The complete circuit is given in Fig. 1. 

It will be seen that the circuit of Gingrich, 
Evans and Edgerton has been altered to include 
an additional tube, (7-6), in the vacuum tube 
voltmeter circuit. This is added to make the 
average output current independent of fluctu- 
ating line voltages on the x-ray filament trans- 
former. It will be observed that, as the current 
through the x-ray tube increases, the grid bias 
on the tube (7-6) increases proportionally with 
a consequent decrease of current delivered per 
pulse by this tube to the electrical tank, (C-20). 
Since the rate of the Geiger-Miiller pulses in- 
creases nearly linearly with increasing current 
through the x-ray tube it is apparent that by 
adjusting R-19, and by suitably choosing the 
sizes of the input condensers to the tube (7-6), 
and (7-7), a setting can be reached such that 
fluctuations of the x-ray filament transformer 
line voltages will not cause appreciable variations 
in the output current. Two electric “eyes” 
(7-4), and (7-5) have been placed in the circuit 
in the manner shown. These tubes serve as 
scales, so that settings of R-19 and the input 
condensers which have given the desired results 
can be duplicated at will. 

The power circuits and the Lifschutz high 
potential circuit are shown in Fig. 2 and Fig. 3, 
respectively. 

The fact that the output current has been 
made independent of the current through the 
x-ray tube removes the last objection to a com- 
pletely automatic x-ray diffraction apparatus. 
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Fic. 4. Balanced filter diffraction—NaCl crystal. e—ZrO2; 
x—SrO. Mo x-ray tube—43 kv peak, 10 ma. 


The Geiger-Miiller counter and the first amplifier 
are mounted on the movable arm of an x-ray 
spectrometer. This arm is driven over the desired 
range by means of a synchronous motor and a 
system of worm gears. This keeps the rotation 
of the arm rigorously in step with the recording 
galvanometer. Obviously, then, the galvanom- 
eter gives a record of x-ray intensities vs. angle. 

Molybdenum Ka radiation was used in con- 
junction with balanced filters of ZrO, and SrO. 
These filters were constructed jointly with P. H. 
Bell of The Pennsylvania State College. The 
close balance obtained is shown in Fig. 4, in 
which the shaded area represents the difference 
between the transmission of the two filters. 

Six diffraction patterns of tricalcium aluminate 
were obtained using filters of ZrO. and SrO on 
alternate runs. Because of the randomness of the 
x-ray quanta registered by the counter, and 
because the apparatus used has an exponentially 
decaying memory, the galvanometer record of 
x-ray intensities shows considerable statistical 
deviations. To minimize this as far as is practical 
the width of the slit system might be increased. 
But previous workers have reported many lines 
relatively close together in the diffraction pat- 
tern, so it is obvious that the slit width should be 
kept as small as possible. If the number of 
readings at a given angle is increased, the actual 
intensity of the diffraction may be obtained by 
the method of averages. The same results may 
be obtained, however, if the angular rate of 
rotation of the spectrometer arm is made very 
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Fic. 6. The diffraction pattern of 3CaO-Al.O;. (The 


slow and a smooth line drawn through the 
recorded diffraction pattern. It is this method 
which is used in this work. 

A very short portion (angular length 20=1 
degree) of a typical diffraction pattern is shown 
in Fig. 5. Reference to Fig. 6 shows that Fig. 5 
actually represents portions of the diffraction 
peaks corresponding to h?+?+P2?=126 and 128. 

The average intensities, as measured by the 
heights of the diffraction record above the zero 
of ordinates, are determined for both sets of 
filters at intervals of one minute of arc 26. The 
differences of the averages are plotted point by 
point on a scale for 26 much smaller than the 
original record and a smooth line is drawn 
through the points. On drawing this smooth line 
apparent ‘‘peaks’’ which obviously do not have 
the characteristics demanded by the slit widths 
are ignored. 

The complete diffraction pattern is shown in 
Fig. 6. The left end of the pattern up to about 
20=4 degrees shows the shadow effect of the 
slits. It is apparent that many portions of the 
pattern contain unresolved lines, the most notable 
of these being in the vicinity of 26=14 to 15 
degrees. Similar studies of the diffraction pat- 


BRADY AND W. P. 


DAVEY 


Fic. 5. A short por- 
tion of a typical diffrac- 
tion record. 
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numbers directly above the peaks refer to h?+k?+/?.) 


terns of the three other calcium aluminates 
shows that no peak in the diffraction pattern 
of Fig. 6 for tricalcium aluminate can possibly 
belong to the diffraction pattern of one of the 
other compounds. In addition to the unresolved 
and partly resolved lines just mentioned, Fig. 6 
shows 22 clearly resolved lines within an angular 
range of 20=42°. A simple cubic structure would 
require, within an angular range of 26=10°, 12 
additional lines. Any structure proposed ior 
tricalcium aluminate must therefore account 
quantitatively for the absence of the 12 ‘“‘for- 
bidden”’ lines, as well as for the presence of the 
22 lines of Fig. 6. The edge of the unit cube 
demanded by the lines shown in Fig. 6 is 15.235A. 
A smaller unit of structure cannot possibly 
account for a large number of the lines in the 
diffraction pattern, while no lines present require 
a larger unit cell. On the basis of its density and 
the size of the unit crystal the chemical formula 
of this compound is (3CaO-Al.O03)04. This means 
there are 72 Ca, 48 Al, and 144 O, a total of 
264 atoms, in the unit cell. 

A highly probable structure has been found 
which accounts for the x-ray dots and which © 
predicts the peculiar properties of 3CaO- Al.Os. 
It will be published elsewhere. 
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The Raman Spectra of Mono- and Dichlorobenzenes{ 


H. SPONER AND J. S. Krrpy-SMITH* 
Department of Physics, Duke University, Durham, North Carolina 


(Received June 9, 1941) 


Raman spectra were taken of gaseous mono- and dichlorobenzenes. Polarization measure- 
ments were made for the dichlorobenzenes in the liquid phase. Interpretations of the stronger 
lines are given in terms of modes of vibration. The experimental results and assignments of 


the frequencies are collected in Tables I-IV. 





INTRODUCTION 


HE Raman spectra of the different chloro- 

benzenes have been thoroughly studied in 
the liquid state by a number of investigators.! 
Polarization measurements have been carried out 
in the case of monochlorobenzene by Simons.’ In 
the course of an analysis of the ultraviolet ab- 
sorption spectra of mono-* and the dichloro- 
benzenes‘ it became desirable to have their 
vibrational frequencies in the gaseous phase. For 
this reason these compounds have been rein- 
vestigated in the vapor phase. To facilitate inter- 
pretation of the observed frequencies polarization 
measurements were made for the dichloroben- 
zenes in the liquid state. 


EXPERIMENTAL 
A. Measurements in the gas 


The substances used in this investigation, 
monochlorobenzene, ortho-, meta-, and para- 
dichlorobenzene, were kindly furnished by the 
Chemistry Department.t They were extremely 
pure, and were carefully introduced into the 
Raman tubes by vacuum distillation. 


{ Presented before the Division of Physical and Inor- 
ganic Chemistry of the American Chemical Society at 
Detroit, September, 1940. 


*Present address: National Institute of Health, 
Bethesda, Maryland. 

1A. Dadieu, A. Pongratz and K. W. F. Kohlrausch, 
Monats. f. Chem. 61, 426 (1932); J. W. Murray and D. H. 
Andrews, J. Chem. Phys. 1, 406 (1933). Complete refer- 
ences may be found in the second book of K. W. F. 
Kohlrausch, Der Smekal-Raman Effekt, Erg. Bd. (Julius 
Springer, Berlin, 1938), or in James H. Hibben, The 
Raman Effect and its Chemical Applications (Reinhold 
Publishing Corporation, New York, 1939). 

_*L. Simons, Soc. Sci. Fennica, Comm. Phys.-Math. 6, 
No. 13 (1932). 

*H. Sponer and S. H. Wollman, Phys. Rev. 57, 1078A 
(1940). 

*S. H. Wollman, unpublished results. 

t We are indebted to Dr. Eunice Moore for the courtesy 
of providing us with the substances. 


The apparatus used, except for modifications 
in the assembly of Raman tube and filter jacket, 
has been previously described.’ Certain changes 
have been made necessary in order to attain the 
relatively high temperatures (200°C in the case 
of the dichlorobenzenes) needed to give sufficient 
vapor pressures for adequate scattered inten- 
sities, and at the same time permit suitable 
filters to be employed. This condition was 
attained by placing an extra glass cylinder 
between the filter jacket and Raman _ tube. 
Strong blasts of hot air were directed between 
this cylinder and the Raman tube and the scat- 
tering section of the tube was heated to the 
desired temperature. The space between the 
filter jacket and the wall of the extra cylinder 
served as an insulating air space and prohibited 
an excessive loss of heat to the circulating filter 
solution. The bottom section of the tube was 
heated by resistance coils which were wrapped 
around it. Temperatures were read with a ther- 
mocouple placed in contact with the coolest 
portion of the Raman tube. With this arrange- 
ment temperatures up to 220°C have been 
maintained at fairly constant values (+5°C) 
during all exposure times. _ 

Twelve Hg arcs were used as the light source, 
excitation being by the 4047 and 4358A lines 
together as well as by 4358 alone. Filters were 
aqueous solutions of CoClz. when excitation by 
both lines was desired, and CoCl,+quinine for 
the isolation of 4358. All spectra have been ob- 
tained upon antihalation Eastman Super Panchro 
Press or Ilford Hypersensitive Panchromatic 
plates. Exposure times have been from 6 to 24 
hours. The spectrograph is the Zeiss 3 prism 
instrument used in previous work. 


5 J. S. Kirby-Smith and L. G. Bonner, J. Chem. Phys. 7, 
880 (1939). 
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668 HB. SPONER AND j.°‘S. KiRBY-SMITH 
TABLE I. CsH;Cl. TABLE II. o-CsH4Cle. 
DEPOLARIZATION Liquip 
LIQuID FACTOR ASSIGN- (SWAINE AND DEPOLARIZATION  ASSIGN- 
GAS (ANANTHAKRISHNAN) (Srmons) MENT Gas MuRRAY) FACTOR MENT 
200 (w) 196 (6b) 0.87 (81) 154 (10) 0.8 
270 (1/2) 194 (w) 203 (3) 0.8 
295 (2) 0.86 239 (1b) 0.8 
416 (m) 420 (6) 0.32 a 300 (m) 330 (0b) 
608 (w) 615 (3) 0.89 Bi 434 (vw) 430 (1) 0.8 
703 (m) 702 (6) 0.13 a 469 (1) 
742 (0) 487 (m) 483 (2) 0.45 a} 
790 (0) 664 (m) 658 (5) 0.35 a1 
830 (0) 760 (w) 756 (0b) 
925 (0) 860 (10) 
989 (0) 1020 (1) 
1001 (st) 1003 (10s) 0.06 ay 1035 (st) 1041 (10) 0.20 a 
1013 (00) 1133 (st) 1129 (5) 0.25 a 
1024 (m) 1024 (7) 0.19 ay 1160 (10) dep (B;) 
1080 (m) 1084 (5) 0.10 a 1271 (w) 1274 (1b) 
1134 (w) 1121 (1/2) (a) 1598 (w) 1577 (4) 0.8 By 
1191 (w) } 1158 (1/2) 0.88 (81) 1607 (0d) 
1176 (1/2) 2994 (1) 
1295 (0) 3082 (st) 3073 (10) 0.40 a 
1321 (0) 3150 (w) 3146 (3) 
1372 (0) my Pee 
1443 ‘oS i —— 
1565 (0 a , 
1589 (m) 1584 (4) 0.87 By TABLE III. m-C.H4Cle. 
3008 (5) = = — 
3028 (0) LIQuIpD 
3079 (st b) 3068 (10d) 0.32 a s (SWAINE AND DEPOLARIZATION _ASSIGN- 
3140 (0) Gas MURRAY) FACTOR MENT 
3180 (w) 3165 (1) 178 (3) ag 
——— a 205 (wb) 202 (3) 0.9 
216 (2) 
356 (w) 366 (1) 
B. Polarization measurements 399 (m) 399 (4b) 0.10 a 
428 (2) 0.6 
The apparatus used here was of conventional 613 (vw) 530 (0) 

° : : 665 (m) 666 (4) 0.25 a} 
type, using a sensibly parallel beam of unpo- 998 (st) 999 (10) 0.20 ia 
larized light as the incident radiation. The light 1018 (0) 

" 1056 (m vb) 1070 (3d) 0.30 

source consisted of two small (15 cm long) Hg 1109 (2) 
arcs similar in design to those used in the gas 1129 (m) 1126 (4) 0.10 a 
: . . (1160 Kohl- dep (B:) 
apparatus. Glass cylinders filled with filter senach) 
solution focused and directed a strong beam of 1240 (0) 
filtered light into the scattering tube. The source oo of 
operated in a vertical position, and a totally 1544 (0) 

. ; ‘ 1590 (m) 1579 (55) 0.8 Bi 
reflecting glass prism was used to direct the 1625 (0) 
scattered light into the polarizing system. This 3090 (st) 3076 (10) 0.45 a1 
consisted of the usual Wollaston prism, a lens ened (w) __ (2) = - 
for focusing the two polarized components onto 
the spectrograph slit, and a half-wave plate. RESULTS 


This arrangement gave good pictures in from 3 
to 8 hours exposure time. 

Under present conditions precise measure- 
ments have not been made. Density marks for the 
calibration of plates and the determination of 
relative intensities have not been put on all 
plates. The non-parallel character of the incident 
beam and imperfections in the half-wave plate 
are the other main sources of error. 


The present results including the assignment 
of frequencies as well as previous liquid and 
polarization data are summarized in Tables I to 
IV. The most recently published values have been 
taken.® Intensities have been determined from 
microphotometer curves. However, due to the 
6 J. W. Swaine and J. W. Murray, J. Chem. Phys. 1, 512 


(1933); R. Ananthakrishnan, Proc. Ind. Acad. Sci. 3A, 52 
(1936). 
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relatively greater amount of Rayleigh scattering 
and general background inherent in gases these 
intensities are only estimated as strong (st), 
medium (m), weak (w) and very weak (vw). 

While nothing has to be added to the tables 
for the dichlorobenzenes, it should be mentioned 
that for monochlorobenzene additional extremely 
weak, broad and diffuse lines have been observed 
at 275 and 500 cm~. These lines appear only on 
several plates, and then with such small intensity 
that they could be spurious. 

The measurement of the frequency in the 
neighborhood of 200 cm~! was taken as an anti- 
Stokes transition from Hg 4047. As a Stokes line 
from Hg 4047 it is masked by the strong Hg 4077 
and from Hg 4358 by the Hg II line at 4398.6A. 
The 200 frequency is usually measured from the 
Hg 4358 line in the liquid, as under these condi- 
tions the spark line is too weak for detection. 

A definite diffuse wing on the low frequency 
edge of the strong 3079 cm~ line has been ob- 
served on all plates. Attempts to resolve any 
structure were unsuccessful. 


DISCUSSION OF RESULTS 


For convenience in discussing the frequencies 
of the investigated benzene chlorides Fig. 1 has 
been included which represents the different 
types of vibration in the benzene molecule.’ The 


TABLE IV. p-CeHsCle. 








LIQUID 
m (SWAINE AND DEPOLARIZATION ASSIGN- 
Gas MURRAY) FACTOR MENT 


302 (4) 0.8 Boy 
333 (8) 0.5 pin 
355 (1) 
627 (5) 0.8 Bog 
710 (00) 
748 (100) 0.2 ag 
885 (0) 
942 (0) 
1070 (3) 0.5 (c1,) 
1087 (2) 
1109 (10) 0.4 arg 
1170 (1) 
1217 (0) 
1294 (0) 
1331 (0) 
1379 (0) 
1576 (8) 0.8 Boo 
1630 (1) 
2953 (0) 
3079 (10) 0.4 ar 
3153 (2) 





326 (w) 
386 (w) 


742 (st) 


1050 (m) 
1110 (st) 


1300 (vw b) 
1440 (vw b) 
1570 (m) 


3090 (st) 








‘Vibration types after Langseth and Lord, reference 9. 


symmetry Dg, of benzene is reduced to V, in 
p-dichlorobenzene, and to C2, in o- and m-dichlo- 
robenzene and monochlorobenzene. The sym- 
metries of the different vibrations of the sub- 
stituted benzenes are represented in Tables V 
to VII. 

The symmetry symbols are listed in Placzek’s 
notation.® C2”, etc. refer to twofold axes in the 
y, etc., direction, ¢, denotes a plane of reflection 
perpendicular to the z axis (molecular plane), and 
i refers to the center of symmetry. The + and 
— signs indicate respectively that a given sym- 
metry class is symmetric or antisymmetric to 
the noted symmetry element. The asterisk at the 
symbol of C2, for the ortho-compound indicates 
that the twofold symmetry axis passes here 
between the two carbons having the substituents. 
Small Greek letters have been used for the sym- 
metry symbols to distinguish them from the 
same symbols with capital Latin letters for elec- 
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Fic. 1. Types of vibrations in benzene. 


8G. Placzek, Handbuch d. Radiologie, Bd. V1/2 (Leipzig, 
1934), p. 283; L. Tisza, Zeits. f. Physik 82, 285 (1932). 
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TABLE V. Symmetries of vibrations for para-Ce6H4Cls in 
group Vj. 








SELECTION RULES 
NUMBER OF 
VIBRATIONS 


ESSENTIAL SYMMETRY 
ELEMENTS 


SYMMETRY CY C2? i RAMAN 





Qig 
Alu 
Big 
Biu 
Bag 
Bou 
Bg 
Bau 


bit+iit+ 
} itr ti+it 
Mw wun mH 
b+i+i+i+ 
Tier iset tt 











tronic terms. Langseth and Lord® have discussed 
in the case of the deuterated benzenes which 
symmetry classes of the group Dg», go over into 
the respective classes of C2, and V,,. 


Monochlorobenzene 


As can be seen from Table I, only the strong 
lines of the liquid have been observed in the 
gaseous phase. Most of them are totally sym- 
metrical vibrations as is indicated in the last 
column under “Assignment.” 

The frequencies 416 and 608 correspond to the 
€,+ vibration of 606 cm! in benzene” (No. 6 in 
Fig. 1). According to polarization measurements? 
the 416 is totally symmetrical and the 608 
non-totally symmetrical. As already pointed out 
by Kohlrausch" the non-totally symmetrical 
vibration is very nearly independent of the 
substitution because here the motion of the 
carbon atom with the substituent is small and 
perpendicular to the C—Cl bond. The totally 
symmetric vibration, however, is strongly in- 
fluenced because the carbon atom attached to the 
chlorine participates in the vibration. 

The 703 represents the vibration of the 
chlorine atom towards the ring. It is totally 
symmetrical and originates out of the four 3000 
cm~! hydrogen vibrations in benzene ay, Bix, 
et, and e~, (Nos. 2, 13, 7, 20 in Fig. 1). The last 
two split in the monosubstituted benzenes into a 
totally symmetrical (a@;) and a non-totally sym- 
metrical (8;) part. We have, of course, to con- 
sider here only the a vibrations. It has no 

*A. Langseth and R. C. Lord, Kgl. Danske Vid. Sels. 
Math.-fys. Medd. 16, No. 6 (1938). 

10 Tn his second volume, Der Smekal-Raman Effekt, Erg. 
Bd., 1938, p. 163, Kohlrausch correlates the Raman lines 
416 and 470 with the 606 vibration in benzene. Our assign- 


ment corresponds to his previous one, cf. reference 11. 
1 K. W. F. Kohlrausch, Physik. Zeits. 37, 58 (1936). 
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meaning to correlate the 703 to a particular one 
of the four totally symmetrical vibrations since 
it will in reality have contributions from all of 
them. The observed frequency 3079 cm~ has to 
be assigned to the remaining hydrogen vibra- 
tions. It is sharp on the short wave side and 
diffuse on the other, and may well cover two 
lines. The frequency of 3180 cm is probably 
due to an overtone. 

The strong lines at 1001, 1024, and 1080 cm— 
are all totally symmetrical according to polari- 
zation measurements.? They are believed to 
correspond to the carbon vibrations ay, Bi, 
and probably e~, in benzene (Nos. 1, 12, 19 in 
Fig. 1). One should expect one of the three 
carbon frequencies Bou, «+, and e~, (Nos. 14, 8 
and 19 in Fig. 1) to decrease considerably in 
monochlorobenzene and this is probably the e-,. 
A detailed correlation of the observed frequencies 
cannot be made because they will be linear com- 
binations of all three modes of vibrations a,,, 
81, and €~,, particularly since the frequencies are 
almost identical. 

The line at 1589 cm is probably one of the 
components that arise from the carbon vibration 
et, = 1596 in benzene (No. 8 in Fig. 1) when the 
degeneracy is released. As the line is found highly 
depolarized it must represent the non-totally 
symmetrical mode of vibration. 

The line at 1134 cm! belongs quite likely to a 
hydrogen vibration, as does perhaps the 1191 
line. The latter may correspond to the depolarized 
line 1158 or the 1176 in the liquid. To find an 
interpretation we consider the six hydrogen 
vibrations of benzene ay, Bou, €*, and €~, (Nos. 
3, 15, 9, 18 in Fig. 1). They all have frequencies 
between 1000 and 1200 cm. In C,sH;Cl one of 
them and in C,sH,Cl. two of them must drop 
down toa very low value. In all compounds there 
must originate at least one totally symmetrical 
vibration around 1100 cm and two non-totally 


TABLE VI. Symmetries of vibrations for CsHsCl and meta- 
CoH Cle in group Co». 
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MONO- 


symmetrical vibrations in the same region. It is 
quite possible that in CsH;Cl the line at 1134 
cm is an a; (with large contribution from e*,) 
and the 1190 cm~ a ;. The vibration with the 
low frequency mentioned above is perhaps the 
200 cm~!. It would be then a £; in accordance 
with the depolarization factor. Another possi- 
bility is that it is the 294 cm observed in the 
liquid. 


Dichlorobenzenes 


Here too only the strong lines of the liquid 
have been found in the vapor. The 606 cm™ et, 
vibration in benzene is split into a,=416 and 
8,=615 cm in monochlorobenzene. The same 
types occur in the meta- and ortho-derivatives, 
but here both will be influenced by the substi- 
tution, the values for the frequencies probably 
not being much apart from each other. Indeed, 
no line has been observed in the neighborhood 
of 600 cm~ for meta- and ortho-C,¢H4Cle. There 
are several lines between 300 and 500 cm™ in 
the meta- and ortho-compounds which may be 
chosen for the assignment in question. Polariza- 
tion measurements indicate that the 399 in 
m-CsH4Cle is totally symmetrical and hence 
quite likely is the a; in question. In o-CgHuCl. 
the strong line at 483 has according to our 
measurements a rather low depolarization factor 
and therefore may be the corresponding ay. 
Kohlrausch" has also made empirically a similar 
correlation. We make no proposal for the non- 
totally symmetrical component. In p-CgH4Cl, the 
606 e+, of benzene splits into an ay, and a Bo, 
vibration. The results in the deutero-benzenes?® 
make it probable that the a, is rather low, lower 
than in CsH;Cl. Hence the most reasonable 
assignment seems to be aj,=326 cm. This is 
further confirmed by the polarization measure- 
ments. The 82, should be higher and is perhaps 
similar to the corresponding 8; in CgHsCl which 
leads to the possible assignment $2,=627 cm=, 
again in agreement with the depolarization 
factor. 

The benzene hydrogen vibrations of high fre- 
quency ayy, Bix, €+, and e~, (Nos. 2, 13, 7, 20 in 
Fig. 1) which give rise in CsH;Cl to one chlorine 
valence vibration of 703 cm! produce two such 
vibrations in the di-derivatives. One of these has 
been observed with certainty in all three cases 
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TABLE VII. Symmetries of vibrations for ortho-CsH4Cl2 in 
group C*>s». 








SELECTION RULES 
INFRA- 
RED 


ESSENTIAL SYMMETRY 
_ ELEMENTS 
SYMMETRY Cot 


NUMBER OF 


VIBRATIONS RAMAN 





ay) 

ors 5 - 
Bi - + 
Be 4 - 








in the lines 742 (para), 665 (meta), and 664 cm= 
(ortho) as has been noticed by previous authors. 
Each line is polarized as it should be for a totally 
symmetrical vibration. The vibration has a 
larger contribution from e*, for the para- 
compound. 

The lines at 3082 (ortho), 3090 (meta), and 
3090 cm (para) have to be assigned, as in 
monochlorobenzene, to totally symmetric hy- 
drogen vibrations. Polarization measurements 
support this conclusion. The frequencies of 3150 
(ortho), 3169 (meta), and 3153 (para, observed 
in the liquid) probably represent overtones. 

As in C.sH;Cl the benzene carbon vibrations 
ai, and @,, (Nos. 1 and 12 in Fig. 1) should give 
frequencies in the neighborhood of 1000 cm™, 
the 81, being Raman inactive in the para-com- 
pound. It is furthermore possible that the e~, 
(No. 19 in Fig. 1) carbon vibration of benzene 
drops to a value in the 1000 cm region. The six 
hydrogen vibrations of benzene a, Bou, e+, and 
e~, (Nos. 3, 15, 9, 18 in Fig. 1) will also give, 
besides two very low frequencies, four frequen- 
cies in the 1100 cm region. In the case of 
p-C HCl. two lines of appreciable intensity have 
been found in the gas at 1050 and 1110 cm~. The 
latter is the stronger of the two. It is very likely 
a totally symmetrical carbon vibration. This is 
supported by polarization measurements and the 
analysis of the ultraviolet spectrum.‘ The 1050 
cm line is probably also totally symmetric 
according to the fairly low depolarization factor. 
No preference can be offered as to the assign- 
ment to a carbon or hydrogen vibration. The 
1170 cm line is perhaps connected with the e+, 
in benzene. 

In meta-C.sH,Cl. three fairly strong polarized 
lines have been observed with frequencies of 998, 
1056, and 1129 cm~'. We would like to consider 
the 998 and 1129 cm~ as belonging to totally 
symmetrical carbon vibrations. This assignment 
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gets support again from the ultraviolet spec- 
trum. The explanation of the 1056 cm line 
seems too indefinite at the moment although its 
low depolarization factor points to a totally 
symmetrical vibration. 

In ortho-C,.H,Cl. two lines were found at 1035 
and 1133 and a weak one at 1271 cm. We 
propose to ascribe the two strong polarized lines 
to totally symmetrical carbon vibrations, again in 
agreement with the analysis of the ultraviolet 
absorption spectrum. The interpretation of the 
1271 cm™ seems ambiguous without additional 
information. We are, however, inclined to believe 
that the depolarized 1160 cm line observed in 
the liquid’ of ortho- and meta-CgH,Cl. cor- 
responds to the depolarized line 1158 cm™ in 
C.H;Cl and hence represents in both cases the 
non-totally symmetrical component originating 
mainly from the hydrogen bending vibration 
e+, = 1178 in benzene. 

It has been mentioned before that out of the 
hydrogen vibrations a2, B2., €*, and e~, (Nos. 
3, 15, 9, 18 in Fig. 1) there must originate two 
low frequencies which represent chlorine bending 
vibrations in the molecular plane. In the case of 
para-C,gH,Cl, only one of them will be observable 
in the Raman effect which we believe to be 
B2,=302 cm originating chiefly from e+, and 
a, and found in the liquid. In meta-CsH4Cl. and 
in ortho-C,H,Cl. there have been observed in 
the liquid three low frequencies of 178, 202, 216 
cm! and of 164, 203 and 239 cm~, respectively, 
while we could detect with certainty only one at 
205 cm= (meta) and one at 194 cm (ortho) in 
the gas. Since in the ortho- and meta-compounds 
one of the above-mentioned two low vibrations 
must be totally symmetrical (a1) and the other 
one non-totally symmetrical (81), we were sur- 
prised to find only high depolarization factors for 
the observed low lines. A detailed assignment 
cannot be given. 

In para-, meta-, and ortho-CsH.Cl. lines have 
been found at 1570, 1590, and 1598 cm, 
respectively. As they are depolarized they cor- 
respond in all probability to the 1589 cm™ in 
C.H;Cl and have likewise to be interpreted as 
the non-totally symmetric one of the two com- 


2 This line has not been observed by Swaine and 
Murray in the meta-compound, but Kohlrausch and col-: 
laborators give it as a weak line. We have confirmed their 
observation. 
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ponents into which the e«*+,=1596 cm in 
benzene (No. 8 in Fig. 1) splits in these molecules. 

The vibrations which we have considered so 
far all take place in the molecular plane, either 
as valence or as bending vibrations. We have 
not mentioned any modes which occur perpen- 
dicular to the molecular plane. Of these the 
carbon vibrations B2, and e*, (Nos. 4 and 16 in 
Fig. 1) in benzene will be influenced only slightly 
by the substitution. They are both Raman 
inactive in benzene, and the e*, is also forbidden 
in the para-derivative. They are allowed in the 
other derivatives considered here. Very few of 
them have been observed in the deuterobenzenes 
and these with very weak intensity. Weak lines 
in the neighborhood of 300 cm in the cor- 
responding chlorobenzenes could be due to these 
vibrations. 

There are four hydrogen bending vibrations 
Qo, Br, €*, and e~, perpendicular to the molecular 
plane in benzene (Nos. 11, 5, 17, 10 in Fig. 1). 
They lie between 670 and 1000 cm-. Out of 
these should result one low frequency for 
monochlorobenzene and two low ones for the 
dichlorobenzenes. The a», is Raman inactive in 
benzene but is known from the infra-red to occur 
at 670 cm™!. In monochlorobenzene a strong 
band has been observed in the infra-red" at 
680 cm which is almost certainly the > result- 
ing mainly from ag, in benzene. The e-, gives 
weak to medium Raman lines in benzene and 
the deuterobenzenes. We suppose that one of the 
lowest frequencies, 200-294 cm in C,sH;Cl, and 
at least one of the three lowest in liquid m- and 
o-CsH.Cle represent vibrations of the type dis- 
cussed here. 

The suggested assignments include most of 
the lines in the gas. Many more have been re- 
ported in the liquid phase. It seems too doubtful 
a procedure, however, to try an explanation of 
these. Some of the interpretations presented here 
have already been suggested by Kohlrausch,” 
while others differ from his. Since the whole 
problem is attacked here from a different point 
of view we believe it worth while to give our 
considerations and conclusions. 

In conclusion we want to acknowledge a 
grant-in-aid from the Duke University Research 
Fund which made this research possible. 


13 J. Lecomte, J. de Phys. et Rad. 8, 489 (1937). 
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The paramagnetic susceptibilities of Ti and Zr have been found to increase linearly with 
increasing temperature from 20°K until a crystal structure change occurs at high temperatures. 
Electrical conductivity, specific heat, and thermal expansion studies were made on Zr and the 
evidence indicates a change in the electron configuration. 


VIDENCE for a phase transition in Ti, Zr, 
and Hf has been obtained by Zwikker,! Vogel 
and Tonn,? and by de Boer and his collaborators.* 
The x-ray evidence‘ shows that at low tempera- 
tures the Ti and Zr have a hexagonal close- 
packed lattice, and that above the transition they 
have a cubic body-centered lattice. Probably the 
same change in lattice structure occurs in Hf at 
its transition point. From our magnetic studies on 
Tiand Zr we believe the mechanism of the phase 
transition may be better understood. 


EXPERIMENTAL 


The chemical analysis of the Ti showed about 
0.015 percent Fe (which means that the magnetic 
susceptibility should be unaffected) and about 2 
percent silicon. The specimens used were solid 
slugs taken from the center of a large metallic- 
looking sample. The Zr analysis showed 0.11 
percent Fe: This dissolved iron contributes to the 
paramagnetism of the Zr only at liquid-hydrogen 
temperatures. The sample was originally pre- 
pared by Professor J. H. de Boer* and consisted of 
large single crystals of ductile metal. Presumably 
there is 1 percent Hf dissolved in the Zr. 

The magnetic susceptibility was determined by 
the well-known method of measuring the force 
which a non-homogeneous magnetic field exerts 
on a specimen. The measurements were made 
with a Bitter type high field solenoid (magnet 
No. 1 of the M.I.T. magnet laboratory).* The 


'C, Zwikker, Physica 6, 361 (1926). 

*R. Vogel and W. Tonn, Zeits. f. anorg. allgem. Chemie 
202, 292 (1931). 

*j. H. de Boer, P. Clausing, and J. D. Fast, Rec. Trav. 
Chim. 55, 450 (1936). J. H. de Boer, W. G. Burgers, and 
J. D. Fast, Proc. Akad. Wet. Amst. 39, 515 (1936). 

*W. G. Burgers, Physica 1, 561 (1934). 

(1939) Bitter and A. R. Kaufmann, Phys. Rev. 56, 1044 
°F. Bitter, Rev. Sci. Inst. 10, 373 (1939). 


adaptation for low temperature measurements is 
a modification of the apparatus developed at 
Leiden and elsewhere.?’ The measurements at 
high temperatures were carried out with the 
specimen enclosed in a small quartz capsule 
which had been previously evacuated and sealed 
off. The chamber in the solenoid into which the 
specimens were hung was continuously flushed 
with hydrogen gas. At liquid- and solid-hydrogen 
temperatures, helium was mixed with the hydro- 
gen to prevent condensation. 

Figure 1 shows the magnetic susceptibility per 
gram of Ti as a function of temperature. There 
are three unusual experimental facts to be ob- 
served from this curve: (1) There is a phase 
transition which appears to be completed at 
1300°K, (2) the transition is continuous such as 
to suggest a ‘‘second-order’’ phase transition, 
(3) the temperature gradient of the susceptibility 
is not in accord with the existing theories of 
paramagnetic metals. The maximum error is 
estimated at 1.6 percent for our magnetic meas- 
urements. The transition temperature reported 
for Ti by de Boer is 1158°K and we concede that 
the 2 percent impurity of silicon could make our 
\-point occur higher. 

Figure 2 shows the magnetic susceptibility per 
gram of Zr as a function of temperature. The 
same unusual experimental evidence has been 
obtained as in Ti. The points observed were com- 
pletely reproducible even in the transition region 
—i.e., no hysteresis with respect to temperature. 
The \ shape of the Zr curve appears more sharp 
than in the case of our Ti and the temperature at 
which the transition takes place agrees with other 
workers.’ It will be noticed that in Zr our curve 
turns up again in the region of temperature 


™M. and B. Ruhemann, Low Temperature Physics 
(Cambridge University Press, 1937). 


673 





SQUIRE AND 


A. R. KAUFMANN 


























200 400 


function of temperature. 








T°K 





| 
1400 





below 50°K, as the experimental points at 20°K 
and 14°K show. However, it is just at these low 
temperatures that the 0.11 percent Fe which our 
specimen contains has a paramagnetism sufficient 
to cause the susceptibility of the Zr to increase. 
Only at 14°K does the susceptibility show a 
significant field dependence. At all other higher 
temperatures, and in particular in the transition 
region, the susceptibility is independent of the 
field. We tried to retain some of the high tempera- 
ture phase by quenching from high temperatures 
but the room temperature value of the suscepti- 
bility remained the same. 

Figure 3 shows our results on the electrical 
resistance of Zr as a function of temperature (full 
line). The ductile Zr was swaged and then drawn 
into a wire about 25 cm in length and 0.38 mm in 
diameter. The wire was annealed in high vacuum 
to relieve strains between successive drawings. 
During the measurements the wire was in a high 
vacuum to prevent gas absorption at high tem- 
peratures. There are two unusual experimental 
facts to be observed from the curve: (1) The 
decrease in resistance becomes appreciable in just 
the region where the ‘‘second-order”’ phase trans- 
ition sets in for Zr as observed in the magnetic 
studies, (2) the slope of the resistance curve 
changes well before the transition region. The 
minimum in the resistance corresponds to the 
peak of the A-point in the susceptibility. The 
dashed curve in Fig. 3 shows the deviation of the 


resistance from a linear function extrapolated to 
high temperatures from the region below 500°K. 
Thus the decrease in resistance is not a sharp 
break but the result of a process going on well 
before the actual falling off occurs. 

Figure 4 shows the specific heat as a function of 
temperature for Zr. The high temperature values 
are taken from Jaeger and Veenstra.* The low 
temperature values were taken by us in an 
attempt to find'a hump in the specific heat of Zr 
which would be similar to that which Simon and 
Critescu® found in Hf at 85°K. We found no such 
hump at these low temperatures for Zr and our 
values are only good to +0.2 cal./mole.!® The 
high temperature results of Jaeger and Veenstra 
are interesting because they show an abnormally 
large specific heat in the region of the transition 
at 1115°K. Their technique consisted in heating 
the specimen to high temperature and _ then 
dropping it into a calorimeter. We have plotted 
their results on ductile Zr which they call speci- 
men “A”’ and we refer the reader to their paper 


8F. Jaeger and W. Veenstra, Rec. Trav. Chim. 53, 922 
(1934). 

®F. Simon and S. Critescu, Zeits. f. physik. Chemie B25, 
273 (1934). 

10 The hump in C, which Simon observed for Hf has 
received considerable attention in recent books on the 
solid state. We believe that a small amount (1 percent) o 
Fe might show magnetic interaction at 85°K and so give 
the observed C, hump. Another possibility would be the 
effect of dissolved hydrogen in the Hf (see K. K. Kelley, 
J. Chem. Phys. 8, 316 (1940)). 
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Fic. 2. Paramagnetic sus- 
ceptibility of zirconium as 
a function of temperature. 

















concerning the uncertainty of the highest two 
values. Certainly there is a phase transition and 
at this temperature either the specific heat be- 
comes infinite or has a hump in the neighborhood 
of the transition. That is, either it is a first- or a 
second-order transition. In view of our results in 
Figs. 2, 3, 5, we feel that, with our specimen, the 
specific heat curve should show a hump in the 
manner which we indicate with dashed lines.. The 
question can only be settled by differential 
specific heat measurements in the transition 
region and preferably with specimens of higher 
purity than thus far reported. 

Finally, we have measured the expansion of Zr 
as a function of temperature (Fig. 5). The normal 
linear expansion coefficient is about 6.5 X 10~® per 
degree centigrade, and remains relatively normal 
up to the transition region. At this point the 
metal contracts on heating through the \-point. 
This latter had been observed earlier in a quali- 
tative way by Vogel and Tonn.? Our results are 
exploratory in nature and represent a maximum 
error of 3 percent. The ductile Zr resembles Cu in 
its elastic property so that we believe the large 
specific heat value at, say, 550°K, is not due 
entirely to the difference: C,—C,=TVa?/X, 
where V=13.9cm* per gram atom, a= 19.5 K10-* 
per degree as the volume thermal expansion 
coefficient, and X=9X10-" cm?*/dyne as the 
compressibility coefficient. For in the present 
case this difference is only about 0.1 cal./mole. 


DISCUSSION 


The electronic structure of Ti, Zr, and Hf is 
similar in that they each have four valence 
electrons with two in the s shell and two in the 
d shell. One would, therefore, expect their 
physical properties to be similar. All three form 
the hexagonal close-packed structure at low 
temperatures having lattice constants! which 
increase with the heavier atoms. 

The magnetic property which we have ob- 
served is not only abnormal in the region of the 
transition but is unusual over the entire range of 
temperatures below the A-point. The suscepti- 
bility resembles the situation which one finds for 
those paramagnetic salts exhibiting antiferro- 
magnetism.'* The susceptibility extrapolated to 
O°K is nearly } of the value at the \-point which 
agrees qualitatively with the Van Vleck" theory 
of antiferromagnetism. The hexagonal close- 
packed crystal forms a layer lattice to conform 


with the requirements of antiferromagnetism. 
One big difference is that in these metals we have 
a crystal structure change, whilst in all cases of 
antiferromagnetism thus far studied there is no 
such change at the A-point. Also, we do not find 


uJ. D. Fast, Rec. Trav. Chim. 58, 972 (1939). 

2H. Bizette, C. Squire, and B. Tsai, Comptes rendus 
207, 449 (1938). C. Squire, Phys. Rev. 56, 922 (1939). 

13]. H. Van Vleck, J. Chem. Phys. 9, 85 (1941). See also: 
L. Neel, Ann. de physique 17, 64 (1932); 5, 256 (1936). 
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against temperature and also 
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Fic. 4. Specific heat per 
mole of Zr. The high tem- 
perature values are taken 
from F. Jaeger and W. Veen- 
stra (reference 8). 
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a field dependence of the susceptibility in contra- 
diction to what the theory seems to predict. 
There is an alternative viewpoint which we 
wish to propose. It would seem that the magnetic 
susceptibility as a function of temperature ought 
to be consistent with the electron theory of 
metals. Existing theory gives x = 26°N(£)[1+--- ] 
where 6 is the Bohr magneton and N(é) is the 
density of states at the top of the Fermi distri- 
bution. The higher terms in the series have the 
effect of slightly decreasing the susceptibility at 





higher temperatures. If we assume that, as the 
temperature increases, the function N(é£) in- 
creases—then the susceptibility would become 
larger. The extra specific heat associated with a 
change of N(£) must be spread over the entire 
temperature range. Professor J. C. Slater has 
suggested to us the following possible mechanism 
which might give just the needed increase in 
N(é): The “‘d” and “‘s’’ energy bands overlap in 
the solid but the d band might be split owing to 
the hexagonal ‘crystalline electric field of the 
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Fic. 5. Linear expansion per 
cm length of Zr. 
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atoms. At low temperatures, a plot of the density 
of energy states against energy will give N(£) to 
be small owing to this splitting. If the splitting 
of the d band is sensitive to lattice spacing, which 
it may well be, then an increase of temperature 
might decrease the width of splitting and the 
magnitude of N(£) should increase. We believe 
this would account for the decrease in resistance, 
as it means essentially the creation of a greater 
number of ‘“‘free’’ electrons. We do not see any 
mechanism correlating the change of N(£é) with 
the existence of the phase transition. 

It is difficult to understand why these metals 
should show a second-order phase transition 
instead of the normal first-order change to be 
expected with crystal structure change. Possibly 
the impurities present in our specimens caused 


the observed second-order transition. If our 


magnetic measurements are to be explained on 
the basis of antiferromagnetism, then we might 
expect a second-order phase transition. However, 
this requires an explanation of why there is a 
crystal structure change and why there is no field 
dependence of the susceptibility. If the magnetic 
results are to be explained on the basis of a 
change in the density of electron states, then we 
might expect only a first-order phase transition 
accompanying the crystal structure change. 

One of us (A.R.K.) participated only in the 
magnetic and resistance studies. We hope to 
obtain some pure Hf at a later date and investi- 
gate its magnetic property. We wish to thank 
Mr. G. H. Chambers of the Foote Mineral 
Company, Philadelphia, for the specimens of Zr. 
We are indebted to Professor J. C. Slater for 
valuable criticism and to Professor F. G. Keyes 
and his associates for the liquid hydrogen. 
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Magnetic and X-Ray Studies on Ti and Zr with Dissolved Hydrogen 


J. Fitzwitiiam, A. KAUFMANN, AND C. SQUIRE 
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(Received June 3, 1941) 


Dissolved hydrogen in Ti and Zr has the effect of changing the lattice of these metals from 
hexagonal close-packed to a nearly cubic face-centered lattice. The paramagnetic susceptibility 
of the system Ti-H and also Zr-H has a value at room temperature different from that of the 
respective pure metals. The dependence of susceptibility on the hydrogen concentration has 
been investigated for Zr in this article. A small trace of iron impurity which was dissolved in 
the Zr was found to precipitate out and regain its ferromagnetic property when the Zr was 


allowed to soak up hydrogen. 


SIMON and J. Aharoni! have studied the 
¢ magnetic susceptibility of the system Pd-H 
as a function of the concentration of the dissolved 
hydrogen and found a linear decrease of para- 
magnetism up to about 0.55 hydrogen atoms to 
each palladium atom. It has been deduced? that 
since the palladium shows 0.55 hole per atom in 
the d band, the dissolved hydrogen is completely 
ionized in the palladium lattice. It is known* that 
both Ti and Zr dissolve hydrogen up to quite 
large concentrations and we have, therefore, 
investigated the paramagnetism of the systems 
Ti-H and Zr-H. We have also studied the crystal 
structure change brought about by the addition 
of hydrogen to these metals and we have noted 
the magnetic behavior of a small amount of 
dissolved iron in the Zr-H system. 


EXPERIMENTAL 


The specimens of Ti and Zr metal which were 
used in the preceding paper‘ were allowed to soak 
up hydrogen by heating the specimens in a 
furnace with an atmosphere of hydrogen. It was 
found that a large amount of hydrogen could be 
dissolved in the Ti at 1000°K and that going to 
higher temperatures drove some of the hydrogen 
out.® Similarly, with the Zr, the specimen began 
to take on hydrogen at 900°K. The dissolved 


1F. Simon and J. Aharoni, Zeits. f. physik. Chemie B4, 
175 (1929). See also B. Svensson, Ann. d. Physik 18, 299 
(1933). 

2.N. Mott and H. Jones, Properties of Metals and Alloys 
(Oxford University Press, 1936). 

3G. Haag, Zeits. f. physik. Chemie B11, 433 (1931). 

*C. Squire and A. Kaufmann, J. Chem. Phys. 9, 673 
(1941), this issue. 

5 A. Sievert and E. Roell, Zeits. f. anorg. Chemie 153, 
289 (1926). 
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hydrogen could be readily removed from the 
metal by pumping at 900°K. We verified that the 


hydrogen did not form a compound but was, in 
reality, dissolved in the metal by determining the 
vapor pressure of hydrogen as a function of 
temperature. It was found that the vapor pres- 
sure was also a function of the concentration of 
the dissolved hydrogen. This means that one 
cannot employ the Clapeyron equation to de- 
termine the heat of adsorption. We did not 
attempt to run adsorption isotherms to obtain 
the true value of the heat of adsorption. 

Normal pure Zr at room temperature has 
a hexagonal close-packed crystal structure 
(c=5.137A, and a=3.227A), but at 1115°K goes 
to cubic body-centered (a= 3.61A).® In Fig. 1(a 
our x-ray powder pattern of Zr at room tempera- 
ture gives agreement with Burger. 

The Zr was filled to approximately 66 atomic 
percent hydrogen, and at room temperature the 
powder pattern showed a tetragonal face- 
centered lattice (c=4.58A, a=4.87A). See Fig. 
1(b). These constants are in good agreement with 
Hagg.’ The introduction of the hydrogen into th« 
specimen was accompanied by an 8.2 percent 
increase in volume. By pumping on the Zr-H 
system at 1000°K, the hydrogen was complete! 
removed and x-ray studies showed a return of the 
metal to the original hexagonal close-packed 
structure. See Fig. 1(c). The lattice constants are 
the same within experimental error. The origina! 
Zr was a large single crystal and the treatment of 
adding and then subtracting the hydrogen caused 
the specimen to become quite polycrystalline. 

With 1.2 atoms of hydrogen per atom of Ti, the 





6 W. G. Burger, Physica 1, 561 (1934). 
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STUDIES ON TI AND ZR WITH DISSOLVED H, 


FIG. 1. (a) X-ray powder pattern of normal Zr showing hexagonal close-packed crystal structure. (b) X-ray powder 
tern of the system Zr-H showing tetragonal face-centered lattice. (c) X-ray powder pattern of Zr after removal of 
olved hydrogen. (d) X-ray powder pattern of Ti-H showing a very nearly cubic face-centered structure. 


system Ti-H was found to have a very nearly Ti-H system showed no magnetic field depend- 
ubic face-centered structure with a=4.45A. See ence of the susceptibility, and this should be 
g. 1(d). carried in mind in understanding the difficulties 
x-ray powder patterns were taken in which we encountered with our Zr-H system. 
cameras having diameters of approxi- In the preceding paper we drew attention to 
tely 57.5 mm. The Cu Ka line was used (the the fact that our Zr specimen contained 0.11 
ile r line was filtered out by Ni) and the tube percent Fe. This iron is dissolved in the Zr and at 
rried 15 ma at 30 kv. No attempt was made to room temperature plays no role in the magnetic 
t high precision on lattice parameters. susceptibility. No amount of heat treatment and 
Since both metal-hydrogen systems went to quenching could make this iron play an impor- 
face-centered forms rather than to the tant role save at liquid-hydrogen temperatures 
(see Fig. 2 of reference 4). However, with 
hydrogen dissolved in the Zr and with suitable 
heat treatment, the iron appeared to precipitate 
out of solution with the Zr in that it gave a 
ferromagnetic component to the magnetization. 
Figure 2 shows a typical magnetization curve for 
The pure Ti showed a susceptibility at 300°K a Zr-H specimen containing 0.015 gram of hydro- 
3.2X10-® per gram, and with 1.2 atoms of gen per gram of Zr. It is evident that there are 
zen per atom of Ti, the system showed an) two components in the magnetization— one 
ise of susceptibility to 3.78X10-® The showing saturation after a few thousand gauss 


ected body-centered structures, it appears 
the binding energy is higher in the face- 
centered form when hydrogen atoms are present. 
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(due to the 0.11 percent Fe) and another part 
showing a normal linear paramagnetic mag- 
netization. There is no hysteresis loop as is found 
for pure iron, but the susceptibility of the 
specimen approaches infinity as the magnetic 
field approaches zero. Figure 3 shows the satu- 
ration magnetization of the Fe component as a 


function of temperature for a specimen which 
contained 0.022 gram of hydrogen per gram of 
Zr. The saturation magnetization at 0°K (extra- 
polated from the measurements at 14°K) is of the 
correct order of magnitude to be completely 
accounted for by the small iron content. All 
specimens containing hydrogen showed very 
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STUDIES ON TI 


clearly a Curie point at 585°K. Above this 
temperature there was no saturation component. 
The Curie temperature for iron is much higher 
at 1000°K. 

Whether or not the concentration of hydrogen 
plays a role other than that needed to change the 
crystal structure we cannot say. The pre- 
dominating factor in causing the Fe to precipitate 
out was the holding of the hydrogen-filled speci- 
men at 1000°K for a period of an hour. Phe- 
nomenologically, we feel that it is very interesting 
that the iron atoms, spaced some 10 atoms of Zr 
apart, can be made to precipitate out and restore 
their ferromagnetic characteristic. We _ believe 
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Fic. 4. Magnetic susceptibility of the system Zr-H as a 
function of hydrogen concentration at 295°K. 
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that somehow the high temperature heat treat- 
ment must have allowed the iron actually to form 
clusters. In all specimens the ferromagnetism 
completely vanished when the hydrogen was 
removed by pumping at 900°K. This implies that 
the Fe has gone back into solution with the Zr. It 
occurred to us that the saturated moment might 
be due to the hydrogen—.e., all protons aligned 
would give a magnetic moment of the order of 
magnitude which we observed. This is most 
unlikely unless one assumes an enormous internal 
Weiss field acting on the protons. We wish to 
emphasize that our Ti, having no Fe impurity, 
showed no such behavior on the addition of 
hydrogen. 

We have felt justified in subtracting off the 
saturated magnetic moment of the Fe at 300°K 
and calling what remained the magnetic moment 
of the Zr-H system. From this magnetic moment 
we compute the magnetic susceptibility and 
Fig. 4 shows the susceptibility per gram of Zr asa 
function of the concentration of dissolved hydro- 
gen. There is some resemblance to the Pd-H sys- 
tem explored by Simon! in that here the sus- 
ceptibility decreases and finally remains nearly 
constant at 2 atoms of H per atom of Zr. The 
change of lattice upon the addition of hydrogen 
makes us hesitate to draw theoretical conclusions 
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from the magnetic studies of either the Ti or 
the Zr. 

Finally, we were interested in what happened 
to the magnetic susceptibility of Zr as the 
hydrogen began to enter the metal. Thus, 
starting with pure, ductile Zr, we began measuring 
the susceptibility as a function of temperature, 
while the specimen was exposed to an atmosphere 
of hydrogen. Figure 5 shows the results. Natu- 
rally, at low temperatures the values are exactly 
those of the pure metal. At 900°K the specimen 
began slowly to take on hydrogen and the 
susceptibility began to increase. Beyond the 
\-point the specimen did not appear to gain or 
lose any appreciable quantity of hydrogen, and it 
is apparent that the susceptibility does not 
resemble that which the pure metal Zr shows in 





RAMA RAO 





its high temperature phase. Comparing the 
\-point with that found for pure Zr in reference 4 
(here represented in dashed lines), it is seen that 
the transition occurs about 200° lower. The 
actual structure change is, however, entirely 
different, which is rather surprising. The meas- 
ured curve in Fig. 5 is not reversible with 
temperature (unless one also pumps out the 
hydrogen) as the new phase is preserved on 
lowering the temperature. 
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A study of the recently obtained data on the velocity of sound and density in liquids reveals 


a simple relation between these quantities namely the law, v! M/p=R where v is the velocity of 
sound in the liquid, M the molecular weight, p the density and R a constant independent of 
temperature. It is shown that the velocity of sound cannot be a proper basis of comparison in 
any homologous series. The result of plotting the constant R against the molecular weight of 
members of homologous series leads to the equation R=aM+8 where a is a general constant 
and @ a characteristic constant for any one homologous series. The difference in R for suc- 
cessive members of homologous series is a constant independent of the series and R is an 
additive function of the chemical structure. Values of R for hydrogen, carbon, oxygen, bromine, 


and chlorine are tabulated. 


STUDY of, recently obtained data' on the 

velocity of sound and density in liquids 
reveals a simple relation? between these quanti- 
ties, namely the law, 


vi'M/p=R, 


where v is the velocity of sound in the liquid at 
the temperature 0, M the molecular weight, p the 
density of the liquid at the same temperature and 


1 Freyer, Hubbard, and Andrews, J. Am. Chem. Soc. 51, 
729 (1929); S. Parthasarathy, Proc. Ind. Acad. Sci. 2A, 
497 (1935); 4A, 59 (1936); W. Schaafs, Zeits. f. Physik 105, 
658 (1937). 
2M. Rama Rao, Ind. J. Phys. 14, 109 (1940). 






R a constant independent of the temperature. 
Since M/p is equal to the molecular volume which 
is proportional to o*, where o is the mean 
molecular distance, it follows* that the velocity 
of sound in the liquid at any temperature is 
inversely proportional to ninth power. of the 
mean distance between the molecules at that 
temperature. As the temperature of the liquid 
increases, o also increases, resulting in a de- 
creased velocity of sound. When the pressure 
increases, « decreases resulting in an increased 


3M. Rama Rao, Curr. Sci. 8, 510 (1939). 
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SOUND VELOCITY 


velocity of sound. These results‘ are true in the 
case of all non-associated liquids. Since R is 
independent of temperature, it follows that a 
change in the velocity of sound with temperature 
is determined by the dependence of the density 
on temperature. The law for the variation of v 
with temperature has been shown! to be 


v=v9(1—0/80.)°*, 


where vo = velocity of sound when 6=0, @.=criti- 
cal temperature on the absolute temperature 
scale. 

The chief importance of sound velocity meas- 
urements in liquids has so far been in their 
rendering possible a determination of the adia- 
batic compressibility. From this and the iso- 
thermal compressibility derived from _ static 
measurements, it is possible to calculate the 
ratio of the specific heats of the liquid. This then 
gives us the value of the specific heat at constant 
volume. The most extensive investigation, under- 
taken with the object of comparing the velocity 
of sound of members of homologous series, has 
been carried out by Parthasarathy.® 

In general, aromatic compounds have higher 
velocity than the aliphatics. For the former the 
range is around 1300 m/sec., while for the 
latter 1300 m/sec. is about the upper limit. 
If we consider the velocity of sound to be the 
proper basis of comparison in homologous series, 
we find that in certain series the velocity in- 
creases with the increase in the length of the 
chain, in others it decreases with the increase in 
the length. In the benzene hydrocarbon series 
the velocity at first increases with the molecular 
weight and then decreases with the higher 
members of the series. The introduction of a 
heavier atom is found to bring down the velocity 
of sound. These qualitative observations are 
shown graphically in Fig. 1, where the variation 
of the velocity of sound of successive members 
of homologous series is studied by plotting the 
velocity of sound against the molecular weight. 
It is noticed that for the saturated hydrocarbons 
the velocity increases with the molecular weight, 
the increase being not proportional to increase in 
the molecular weight. In the case of esters there 

* J. C. Swanson, J. Chem. Phys. 2, 689 (1934). 


°M. Rama Rao, Ind. J. Phys. 14, 109 (1940). 
°S. Parthasarathy, Curr. Sci. 5, 11 (1937). 


IN LIQUIDS 
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Velocity of Sound 








(30.1000) 50 70 90 110 


Molecular Wei ht 
ne 


Fic. 1. 1—paraffins; 2—esters of acetic acid; 3—mono- 
hydric alcohols; 4—benzene hydrocarbons; 5—ketones. 


is a diminution in the velocity with the increase 
in the molecular weight. In the case of mono- 
hydric alcohols and ketones it is found that the 
velocity increases but not regularly with the 
increase in molecular weight. In the case of 
benzene hydrocarbons there is a marked maxi- 
mum of the velocity of sound for xylene. For 
higher members of the series, the velocity de- 
creases but not regularly. The graph shown in 
Fig. 1 is drawn at one temperature for all the 
liquids. Hence it is not possible to get any 
insight into the chemical constitution of the 
substance having as our basis of comparison the 
velocity of sound. 

The simple law between the velocity of sound 
in liquids and molecular volume v!M/p=R, 
shows that a comparison of the values of R of 
two liquids is equivalent to the comparison of 
molecular volumes of these liquids at tempera- 
tures at which the acoustic velocity is the same 
in those liquids. The constant R for a liquid 
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TABLE I. 
SeRIES Ligup Formuta Density Meee. R fie 
Paraffins oe Pentane CsHi2 0.6320 1052 1160 7 
Hexane CcHis 0.6580 1113 1356 mi 
Heptane CrHie 0.6812 1165 1545 nt 
Octane CsHis 0.7020 1238 1746 = 
Esters of acetic acid Methyl =~ ~-CH;COOCH, ~—0.9282,—ts«a2' “aes 
Ethyl CH;COOC;H; 0.8989 1187 1037 wi 
Propyl CH;COOC3H; 0.8910 1182 1211 . 
Butyl CH;COOC,Hs 0.8710 1179 1408 a 
Amyl CH;COOC;Hi 0.8600 1168 1594 ™ 
Monohydric alcohols -—-Methy! CHOH  =————si0.7920 11300 
Ethyl C.H,OH 0.7863 1207 624 aia 
Propyl C;H;,0H 0.8010 1234 806 “ 
Butyl C,H»OH 0.8082 1315 1004 we 
Amyl CsHuOH 0.8136 1347 1198 a 
Benzene hydrocarbons Benzene nein GHe —_— 0.8725 a 1310 ; 979 a 
Toluene C;Hs 0.8631 1320 1170 o 
Xylene CsHio 0.8580 1330 1362 ” 
2X 185 
Cymene CioHi4 0.8470 1308 1731 
Ketones Acetone C;H,O me 0.7910 _ 1203 = 781 — ae 
2X 190 
Diethylketone CsH100 0.8132 1314 1161 . 
Methyl hexylketone C;Hi.O 0.8170 1324 1720 ii 
Methylene chloride CHCl  ~—«*1.3360 1064s ; 
Ethylene chloride C.H.Cl. 1.2550 1240 846 i 








The average difference which corresponds to the value of R of the CH2 group is found from the examination of several series to be 195. 


may be regarded as the molecular volume of the 
liquid when the velocity of sound in it is unity. 
It is well known that molecular volumes of 
liquids at the boiling point are an additive 
function of the molecular volume. Since R can 
be regarded as a molecular volume when the 
velocity of sound in the liquid is unity, it is 
reasonable to suppose that the same additivity 
would apply to the constant R. This has proved 
to be definitely the case. In Table I are given 
values of R for members of homologous series. 
It is found that the difference in the value of R 
between successive members of a homologous 


series is a constant and independent of the type 
of compound. 

In Fig. 2, values of R are plotted against the 
molecular weight of successive members of homol- 
ogous series. It is seen that they are all straight 
lines having the same slope. In the first instance, 
the graph shows that for every addition of a CH: 
group the constant R increases by the same 
amount. Since all of them are a set of parallel 
lines, we can represent each one of them by the 
equation 

R=aM+48, 


where a is a constant independent of the series 
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SOUND 


and @ is a constant characteristic for any one 
homologous series. The relation holds well for a 
number of series. In Table II are given the 
values of a and 8 for a number of series. 

By subtracting from the value of R for a 
paraffin (formula C,,Hoe,+2) that for mCH2 groups, 
i.e., 195, the value of R of two atoms of hydro- 
gen can be found as in the case of the parachor’ 
and the same considerations that apply to the 
case of the parachor also hold good for the case 
of R. The average difference for the series is 
thus found to be 185, and so the value of R for 
hydrogen is taken as 92.5. Since the value of R 
for CHe is 195, it follows that the value for the 
carbon atom is 10. Proceeding in this manner 
it is possible, from known values for carbon and 
hydrogen, to determine the values of R for other 
atoms. For example, that for oxygen is obtained 
as 74 from measurements on ethers. It will be 
observed that no value of R has been attributed 
to a single linkage. Although this is not neces- 
sarily correct, the atomic values of R given for 
hydrogen, carbon and oxygen really include the 
contribution of such a linkage. The single bond 
is taken as an arbitrary zero level from which 
the effect of other structures is reckoned. As with 
other additive properties, so with the constant R, 
constitutional factors must be taken into ac- 
count. The same difference is obtained for all 
substances containing a double bond and a 
mean value of 110, for the contribution of R 


TABLE II. 








SERIES a B 


13.97 +155 
14.01 — 206 
14.00 — 32 
14.02 —120 
13.99 — 48 





Paraffin 

Esters of acetic acid 
Monohydric alcohols , 
Benzene hydrocarbons 
Ketones 








It is to be noticed from the above table that @ is nearly the same for 
all the series examined with a mean value of 14, and 8 varies consider- 
ably from one series to the other. It is positive in the paraffin series 

ut has a negative value in other series. 


*Sugden, The Parachor and Valency (G. Routledge & 
Sons, London, 1930). 
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Fic. 2. 1—paraffins; 2—esters of acetic acid; 3—mono- 
hydric alcohols; 4—benzene hydrocarbons; 5—ketones. 


TABLE III. 
Element R 
Carbon 10 
Hydrogen 92.5 
Oxygen 74 
Chlorine 227 
Bromine 245 


Double Bond 110 


is attributed to the presence of such a bond. 
In Table III are given some of the values ob- 
tained for different atoms and structures. When 
more data on the velocity of sound are available, 
it is seen from the above considerations that it 
should be possible to extend these investigations 
to the problems of chemical structure. 

In conclusion it gives me great pleasure to 
record my thanks to Professor A. Venkata Rao 
Telang, for his encouragement throughout this 
investigation. 
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The Thickness of the Rigid Water Film at a Quartz-Water Interface from a 
Measurement of Newton’s Rings 
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The thickness of the water film between a quartz lens and disk was determined by measuring 
the radii of successive Newton’s rings. The squares of these radii are a linear function of the 
order of the rings from the center, the intercept of the extrapolated graph of which is a function 
only of the radius of curvature of the lens and the vertical separation of the lens and plate. 
By comparing intercepts obtained from data taken before and after evaporation of a water 
layer, a value of the order of 100A was obtained for the thickness of the water film of each 


interface. 





HE postulate of rigid water films, though 
frequently made, has not met with uni- 
versal acceptance. However, J. M. Macauley! 
found a greater shear force for wet than for dry 
glass plates, 1.5X10-5 cm apart, which he 
attributed to rigidity in the water film. Also, B. 
Derjaguin? found a rigidity in water films as 
much as 1.5X10-5 cm thick. He determined the 
rigidity of the water layer by its effect on the 
period of a torsion pendulum. In a previous 
paper,’ the authors showed that the assumption 
of a rigid water layer 9A to 63A thick at a solid- 
water interface could account for experimentally 
determined electrokinetic potential-concentration 
curves without assuming large specific adsorption 
potentials for the alkali and halide ions. 

It was decided to investigate film thickness by 
a study of the interference rings formed in 
monochromatic light reflected from a flat disk 
resting on a convex lens. A rigid water film 
should support the small weight of the disk 
indefinitely, whereas a viscous film would not. 

A few qualitative observations serve to illus- 
trate some possibilities and limitations of this 
method of film measurement. The diameters of 
Newton’s rings fluctuate greatly for small differ- 
ences in the separation of the surfaces producing 
them. It was difficult to get dry, clean surfaces 
close enough together to produce a central dark 
spot. The upper disk would slide easily over the 
lens. A slight tapping on the lens support would 
cause the diameter of the central ring to increase 
or decrease in size or even disappear, indicating 

1J. M. Macauley, Nature 138, 587L (1936). 

2 B. Derjaguin, Zeits. f. Physik 84, 657 (1933). 


3 W. G. Eversole and Paul H. Lahr, J. Chem. Phys. 9, 
530 (1941). 


variations in height of the disk of the order of 
1000A. These variations in ring size may have 
been due in part to dust particles, it being 
difficult to get a dust-free dry quartz surface; or 
to imperfections in the surfaces, 1000A being 
within the tolerance allowed in lenses. 

Perhaps, also, an air cushion was present. If 
the lens and plate were assembled inside a glass 
tube and evacuated for 24 hours at 110° with a 
Hyvac and mercury vapor pump, no appreciable 
change could be observed in ring size. However, 
if water vapor was admitted to the tube, it would 
condense spontaneously between the quartz sur- 
faces. On removal of the water by evacuation, the 
disk would make closer contact with the lens than 
before. It would no longer slide easily when the 
lens was jarred, even though exposed to the air. 
And only by forcibly separating the surfaces 
could the initial conditions be obtained. It was 
difficult to remove the water from between the 
lens and disk and it is not certain that the 
prolonged evacuation used resulted in removal of 
all the water. 

An attempt to measure an air film thickness 
was discarded, since the change in ring diameters 
before and after wetting the lens was not con- 
sistent from one trial to the next. On the other 
hand, determination of a water film thickness by 
comparison of ring diameters before and after 
pumping off the water from the lens gave more 
consistent results. 

The formula used follows simply from the 
optics of Newton’s rings.* Thus: 


€n=nd/ 2p cos 8, (1) 


4R. A. Houston, A Treatise on Light (Longmans, Green 
and Co., New York), pp. 144, 145. 
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WATER FILM AT 
where eé, is the vertical distance between a lens 
and plate at the mth ring, \ is the wave-length of 
light, w is the index of refraction of the fluid 
between the surfaces, and @ is the angle of 
incidence of the light on the lens (in a plane at 
right angles to the plane of the diagram). (See 
Fig. 1.) Also: 


2/2R, (2) 


where a, is the height of the arc subtended by the 
chord, 7», and R is the radius of curvature of the 
lens. 

From the diagram it may be seen that the 
minimum separation of the lens and plate ¢ is 
expressed by 


An=?ln 


t=€n—Qy. (3) 


Substituting the values of e, and a, (Eqs. (1) 
and (2)) in Eq. (3) yields the formula 


t=nd/2pu cos 0—7,?/2R 


7,7 =n\R/p cos 0—2Rt. (4) 


From Eq. (4) it is evident that a graph of 7,2 vs. n 
should give a straight line having an intercept on 


the n=0 axis of (—2Rt). Thus ¢ can be calculated 
independently of \, @ and u. 

Both formulas (1) and (2) involve approxi- 
mations which require that R be large compared 
to r, and a,. The fact that the points fall very 
nearly on a straight line for the first ten or fifteen 
rings is taken as sufficient evidence that the above 
conditions were met. Any optical errors intro- 
duced by irregularities in the lens or tube 
employed were eliminated or diminished by 
subtracting, from the ¢ obtained for a water film, 
the ¢ obtained over the same portion of the lens 
with the film removed. 

The apparatus used is shown in Fig. 2. An 
optically polished flat quartz plate 25 mm in 
diameter and 5 mm thick rested on the convex 
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surface of a plano-convex lens of the same 
diameter, 2-meter radius of curvature, and 3 mm 
thick. A hemicylinder of glass was cemented with 
balsam to the bottom of the lens in order to 
deflect the undesirable light reflections from the 
lower surface of the lens. The lens and plate 
rested on a sheet copper bridge which could be 
slid into place in a Pyrex tube. The Pyrex tube or 
jacket could be opened at a ground glass joint, 
had a stopcock at one end for admitting water or 
dry air, and a stopcock at the other end which 
was connected to a Hyvac and mercury vapor 
pump capable of producing a vacuum of 10~* cm. 
The tube rested horizontally across a heater 
element. A slit source of sodium vapor light was 
used to produce the Newton rings, diameters of 
which were measured with a low powered 
traveling microscope set at about a 30° angle with 
the axis of the lens. Dry air was admitted from a 5- 
foot 2-inch glass tube filled with glass beads and 
P.O;. Water was admitted into the tube by a 
pipette drawn out to a capillary at its end so that 
it could be thrust through the open stopcock. 

The general experimental procedure was to 
clean the lens, plate and tube, slide the lens into 
place in the Pyrex tube, admit a few cc of 
distilled water into the tube, evacuate most of the 
air and allow water to condense between the lens 
and plate. Although condensation of water be- 
tween the lens and plate occurred spontaneously, 
the process could be hastened by suddenly heating 
the jacket (with the heater element and hot 
towels) to raise the vapor pressure. Once water 
condensed on the lens, the water was evaporated 
from the tube by prolonged pumping. Dry air 
was admitted from the P.O; tube, the light and 
microscope were adjusted and the diameter of the 
first twenty rings recorded three times. Then the 
water was again admitted, the air evacuated and 
a water film allowed to condense. Air pressure 
was again brought to normal by carefully opening 
the stopcock and the diameter for the first twenty 
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TABLE I. Calculation of At by the method of least squares. 


_ 2(n)z (nr?) —Z(n*)z(r*) 
~  [S(n)P—NZ(n?)  ’ 


ro>=—2Rt, R=200cm, t= —7r,?/400, 


2 





2 2 
__ f0(air) —70(water) 

















At= 400 » At=3.3X10-* cm. 
2 2 2 2 
n 1 (water) NY (water) ’ (air) NY (air) 
zZ 0.0304 0.0607 
3 0.0293 0.0879 .0429 .1288 
4 .0382 .1529 .0560 .2241 
5 .0473 .2366 .0688 3438 
6 .0562 3373 .0818 4908 
7 .0650 4548 .0945 .6615 
8 .0735 .5880 1072 .8576 
9 .0821 .7387 .1199 1.0791 
10 .0909 .9090 -1324 1.3240 
11 .0995 1.0943 
12 -1085 1.3020 
13 A972 1.5236 
14 .1254 1.7556 
Total 0.9331 9.1807 0.7339 5.1704 

















rings in the water film recorded three times. As a 
check, the water film was evacuated, air admitted 
and another recording made with the dried lens. 

The values of 7,” obtained from the data when 
plotted against n, the number of the ring, gave 
very nearly straight lines. The points for the first 
few values of 1 were more variable due to the 
difficulty of measuring the more diffuse central 
rings. The line curved slightly toward the n axis 
for higher values of m due, perhaps, to the 
approximation involved in the formula a=r?/2R 
and to the fact that the rings were observed 
through a non-optical Pyrex tube. Therefore, ¢ 
was determined by the method of least squares 
for the 2nd to the 10th ring for the air layer and 
for the 3rd to the 14th of the smaller water-layer 
rings. This involved readings taken over the same 
section of the lens. Values of At were obtained by 
subtracting ¢ for the dry surfaces from the ¢ 
obtained for the water film. The data and 
calculations for the first trial are shown in 
Table I. The results of five successive experiments 
are shown in Table II. The average At obtained 
was 230A. This involves two interfaces, so that 
the rigid film at a single interface would be At/2 
or about 100A, subject of course to considerable 
experimental error. This value is to be compared 
with the films of 9 to 63A reported by the authors 
from a study of zeta potentials and with 750A or 
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less reported by Derjaguin.*® Figure 3 shows the 
appearance of the two curves of trial 5 when 
extrapolated to the n=0 axis. 

Several factors were considered which were 
thought to have a possible bearing on the validity 
of the results. The intercepts representing 7,” 
were positive in all but one instance, indicating a 
negative value for the uncorrected ¢. Whatever 
the cause of these negative values (such as a flat 
spot on the lens or a slightly concave disk) it was 
assumed to affect ¢ for the wet and dry lens 
equally and hence was eliminated in the values 
of At. 

The surface tension of the water layer between 
the lens and disk is a factor which could not be 
similarly eliminated. However, it can be shown to 
be an influence opposing the formation of a rigid 
film. Water spreads spontaneously on a flat 
quartz surface even though such spreading re- 
quires an increase in the free energy of the air- 
water interface. The spreading of water between 
the lens and disk involved two quartz-water 
interfaces and a relatively smaller increase in the 
water-air surface. Consequently, the effect of 
surface tension was to augment the weight of the 
disk in compressing the water film. 

According to Langmuir,’ a repulsion would be 
expected between quartz surfaces in water due 
to the osmotic pressure of excess ions held be- 
tween the surfaces by a potential resulting from 
“overlapping’’ of diffuse double layers. If the 
separation of the quartz surfaces ¢ is a result of 
such repulsion, ¢ should be a function of concen- 
tration of electrolyte. To test the effect of 
concentration on ¢, the lens and plate were 
mounted above a dark glass plate in a three-inch 
evaporating dish and conductivity water added 
until the meniscus nearly reached the top of the 
upper plate. Diameters of the first ten successive 


TABLE II. Results. 

















TRIAL rocwater) ro(air) AtX 10° cm 
1 0.00349 cm? | 0.00482 cm? 5.3 
Z .00318 .00381 1.6 
3 — .00057 .00034 2.3 
4 .00335 .00463 3.2 
5 .00088 .00141 1.3 
Av. 2.3 











5 Irving Langmuir, J. Chem. Phys. 6, 873 (1938). 
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Newton’s rings were measured and solid potas- 
sium chloride added in sufficient quantity to raise 
the concentration of the solution above 0.1 molar. 
Readings were taken on successive days, care 
being taken to maintain the original water level 
by addition of distilled water. The only change 
after two to four days in each of five trials was a 
slight increase in ¢. 

Because the diameters of the rings were de- 
pendent upon the level of the water in the dish, a 
slightly different procedure was tried. The lens 
and plate were mounted above a beaker and 
water solution was allowed to siphon through a 
capillary tube into a tiny wad of cotton which 
made contact with the water between the quartz 
surfaces. The rate of flow through the capillary 
was slightly greater than the evaporation rate so 
that excess solution dropped from the lens every 
few minutes. The lens and disk were pressed 
together while wet with conductivity water and 
placed in contact with 0.01 molar thorium nitrate 
solution from the siphon. The diameters of the 
first ten interference rings were measured daily 
for four days. A slight increase instead of a 
decrease in ¢ was observed, although evaporation 
of the solution no doubt increased the concen- 
tration of the solution many times. Since no 
decrease in ¢t was found, it is assumed that the 
plates were held apart by other than the 
“Langmuir effect.” 

Several observations add weight toa hypothesis 
of a water film between the quartz surfaces. The 
midpoint of the central dark spot in the Newton 
ring pattern is dark gray when water wets the 
surfaces. When the water is pumped off, the 
center of the dark spot becomes very black as 
though a true quartz to quartz contact were 
obtained. Also, when the two disks rest in a dish 
of water, the diameter of the observed rings is 
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quite sensitively dependent upon the depth of the 
water, and hence upon the buoyant effect and 
surface tension of the water on the upper disk. As 
the water in the dish becomes deeper, the rings 
get smaller, indicating a greater separation of the 
lens and plate. 

It appears from the results of this experiment 
that water adsorbs on quartz in a rigid multilayer 
the minimum thickness of which is of the order of 
100A. Perhaps surface forces and the potential 
gradient at a quartz-water interface acting on 
water dipoles are sufficient to give rigidity to the 
structure of water in this region. 
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The Production of Concentrated Carbon (13) by Thermal Diffusion 
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A thermal diffusion column 74 feet long and employing methane is being used for the con- 
centration of C8, Every 48 hours 308 milligrams of methane are removed in which the C#/C”® 
ratio has been increased 11.5 times. The power consumption is 6490 watts and six liters/minute 
of cooling water is used. An expenditure of 156 kw-hr. of energy thus results in a transport of 
13.5 mg of C" per 24 hours. The measured output is in good agreement with the value com- 


puted from the theory. 





HE thermal diffusion method of Clusius and 
Dickel' has been used in a number of cases? 
to produce methane in which the C® content 
was increased. Two years ago a 24-foot column 
was constructed in this laboratory with which it 
was possible to obtain a satisfactory check*‘ of 
the theory of operation proposed by Furry, 
Jones and Onsager.® 
The question arose as to how this column 
might be extended to become useful for the 
production of methane in which the C" content 
had been increased from a normal 1.1 percent to 
at least 10 times this value. The fact that brass 
tubing is made in lengths of approximately 12 
feet and that there was available a vertical 
space of approximately 36 feet suggested the 
use of a 72-foot column broken into two parts. 
Furry, Jones and Onsager® have shown that 
if +s represents the transport of gas along the 
tube and if C, is the concentration of C"H, at any 
point along the tube and if the rate of removal of 
CH, from the lower end of the column is 
C,“r one may write 


1 dC; 
ibe C1 C)-——"] + Gar, (1) 
2A dz 


where A and H can be computed in terms of 
various constants for the gas, and the pressure, 
temperature and dimensions of the column. 

The solution of (1) subject to the condition 


1K. Clusius and G. Dickel, Nature 26, 546 (1938); 27, 
148 (1939). 

2H. S. Taylor, Nature 144, 8 (1939); T. I. Taylor and 
G. Glockler, J. Chem. Phys. 7, 851 (1939); 8, 843 (1940); 
W. W. Watson, Phys. Rev. 56, 703 (1939); 57, 899 (1940). 

3A. O. Nier, Phys. Rev. 57, 30 (1940). 

4]. Bardeen, Phys. Rev. 57, 35 (1940). 

5 W.H. Furry, R. C. Jones and L. Onsager, Phys. Rev. 
55, 1083 (1939). 


that C; is small is 
Ci"/CyY=(1+n)/{n+exp [—(1+n)2AL]}, (2) 


where n=7/H and L is the length of the column. 
Table I gives the computed concentration and 
transport under various operating conditions for 
a 72-foot column similar to the 24-foot one 
previously investigated. 
In substituting in Eq. (2) 2AZ was taken as 


(6.58/p?) /(1.63+0.208/p%), (3) 


where is in atmospheres pressure. This value is 
three times the empirically measured value*® for 
the 24-foot section. H was computed from Eq. 
(30) of reference 5. 

The concentrations given are all somewhat 
high as they are based upon the assumption that 
C:2C,(1—C,). Only in the cases of the high 
concentrations would this discrepancy be serious. 
It appears that operation at a pressure of one 
atmosphere with a drain of 0.150 g of CH, per 
24 hours would give sufficient concentration with 
reasonable economy. Moreover operation at a 
pressure of one atmosphere or slightly over 
would minimize difficulties should the system 
have any small leaks. 


L 
TABLE I. Approximate separation factor oe for C8 and 


transport in mg of C™ per day for average C,°=0.01. 











PREs- 
SURE 
(ATMOS.) 0.6 0.8 1.0 1.2 
Ci” Trans- Cr” Trans- Ci” Trans- Ce Trans- 
Cy port CG, ~—s«éPort Cy ~—«~iPort Cy ~—«~éport 
g CH, 
— 
per day 
0.05 16.9 5.4 28.1 8.2 22.4 6.9 12.0 4.0 
0.10 9.1 6.2 16.3 10.5 15.7 10.2 10.6 Son 
0.15 6.53 6.9 10.3 10.5 12.2 12.2 9.58 9.8 
30 3.76 8.2 5.82 12.4 7.65 16.0 6.93 14.6 
0.45 2.85 9.3 4.22 13.6 5.70 18.2 6.11 19.4 
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CARBON (13) BY 

From the standpoint of maximum power 
economy it would be desirable to use a column 
for which the average circumference of the 
annular gas-containing space varied approxi- 
mately inversely with the concentration along 
the column.*® Such a column would not be 
practical to build so one would reduce the size 
by a number of stages having progressively 
smaller cross sections from top to bottom. Such 
an arrangement would not give quite the output 
of a continuous large column but would require 
considerably less power. 

The arrangement used in the present work is 
shown in Fig. 1. Table II gives the dimensions 
and other pertinent data concerning the several 
sections. Section 1 is the original one* and 
section 2 was an extension added later. The 
combination operated at a pressure of 0.8 atmos- 
phere with the hot walls at 300°C and produced 
170 milligrams of CH, per 24 hours in which the 
C® content had been increased by a factor of 6.0. 
In the fall of 1940 section 3 was added, the tem- 
peratures of the hot walls raised to approxi- 
mately 380°C and the pressure to 1 atmosphere. 

Section 1 is connected to the 150-liter reservoir 
by means of two leads of ordinary }-inch copper 
refrigerator tubing. A pump P, circulates the 
gas. This consists merely of an iron cylinder 
which moves up and down with a period of 
approximately one second in a brass tube sur- 
rounded by solenoid actuated by a. multivibrator 
oscillator. A small piece of rubber dam over a 


TABLE II. Data on the three sections of column. 











No. 1 No. 2 No. 3 

Length (ft.) 24 14 36 

(cm) 730 427 1100 
Diameters (inches) 
Steel heated surface 
1.375 0.750 0.375 
i. BD. 1.305 0.312 0.277 
Brass cooled surfoce 
O. D. 2.000 1.375 1.000 

: L. D. 1.936 1.291 0.928 

Spacing ‘‘d’’ between 

hot and cold walls 

(cm) 0.712 0.69 0.71 
Average circumference 

“B"” of annular 

_. Space (cm) 12.5 8.13 5.23 
Size Nichrome V 

heating wire (B. 

: and S. gauge) 8 14 14 
Voltage 61.6 61.6 109 
Current (amp.) 55.2 21.9 15.9 
Power (watts) 3400 1350 1740 
Temp. °C. Hot wall 380 385 367 

Cold wall ~27 ~27 ~27 
2AL 1.21 0.78 1.82 
H (g/sec. of CHa) 3.24 X1075 1.89 X10-5 1.36 X10-5 
Volume of upper reservoir 150 liters 
Volume of lower reservoir 500 cc 








°W. Krasny-Ergen, Phys. Rev. 58, 1078 (1940). 
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Fic. 1. Schematic drawing showing arrangement 
of sections of column. 


hole in the iron cylinder and a similar valve in 
an obstruction in the brass tube insure the 
unidirectional flow of gas in the tube. The pump- 
ing speed of several cc per second is sufficient to 
maintain equilibrium between the column and 
reservoir. 

Section 2 is connected directly to section 1 by 
means of a short brass tube (length=5 cm, 
diameter =1.8 cm). Because of the limited ver- 
tical space available section 3 is set next to 1 
and 2 and is connected to 2 by two leads of }-inch 
copper refrigerator tubing one of which contains 
a pump P? similar to P. 

The construction of sections 2 and 3 was 
similar to that of 1 except that the steel tubing 
constituting the hot walls had sufficient thick- 
ness to eliminate the necessity for the extra tube 
used in section 1 to center the heating wire. 

The methane is kindly supplied by the Southern 
California Gas Company and is over 99 percent 
pure. Every 48 hours 550 cc of gas is removed 
from the bottom by attaching a bulb at the 
point marked C between the valves. B is a 
small appendage which contains KOH pellets 
for the removal of CO». As the column is ex- 
tremely efficient in separating heavy impurities 
in the methane the output is only approximately 
75 percent methane. The equivalent steady drain 
of concentrated methane is thus 6.4 mg/hr.: 
The C/C® ratio in this is 11.8 percent. Every 8 
days the contents of the upper reservoir are 
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changed after closing the valves between column 
and reservoir. The C" in the reservoir is never 
depleted by more than 15 percent. The average 
C83/C” ratio in the reservoir is thus 0.0103. 
The average concentration produced by the 
column is hence 11.5 times. 

The C*/C® ratio at the top of section 3 is 
0.030. Sections 1 and 2 together change the con- 
centration by a factor of 2.9 whereas section 3 
produces a change of 4.0 times. 

These data are not sufficient to determine 
empirically both the separation factor exp (2AL) 
and the transport (/7) for the different sections 
of the column. It is of interest, however, to 
compare the observed changes in concentration 
with those expected from theoretical considera- 
tions. The separation factor for the 24-foot 
section can be estimated from the empirical 
expression found earlier, i.e. 4 of (3). 

This expression was obtained when the tem- 
perature difference between the walls was 300°C. 
Since the operating temperature difference is now 
about 350° a slight change is necessary. The 
numerator and the first term in the denominator 
should be more or less independent of tempera- 
ture, while the second term in the denominator 
should vary roughly as (7 /AT)?, where T is the 
mean temperature and AT the temperature 
difference. The value of 2A,L; obtained from (3) 
with this correction from the change in tempera- 
ture is noted in Table I. The same expression 
has been used to estimate the expected values 
for 2A2L2 and 2A3L;3 for the second and third 
sections. The value of A for the third section 
was taken to be the same as the A for the first 
section. A correction for the difference in spacing 
(d) was made in obtaining 2A2L»2 for the second 
section. These estimates assume that the cor- 
rection for asymmetry in the temperature dis- 
tribution is the same for all three columns. 
Without further information about the operation 
of the columns, this is probably as good a guess 
as any. The transport coefficients, H, as com- 
puted from the theoretical expression’ of Furry, 
Jones and Onsager, are also listed for the three 
sections. 

It is not difficult to compute from these 
values and a combine formula similar to (2) the 


7 Because of the numerous approximations made, it was 
not felt justified to correct H, A, and Kp/K to allow for 
the cylindrical construction of columns. See W. H. Furry 
and R. C. Jones, Phys. Rev. 57, 561 (1940). 
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concentration to be expected for any given rate 
of drain. If methane is drawn at a steady rate 
of 6.4 mg/hr., the calculations yield a concen- 
tration of 12.5 times for the whole column and 
of 3.0 for the first two sections, as compared 
with the observed 11.5 and 2.9, respectively. 
The calculated values are thus a little higher 
than those obtained in practice. 

The discrepancy between practice and theory 
could be due to a number of factors. The rate of 
drain is not constant but intermittent. This 
means that each time a sample is removed the 
lower end of section 3 is disturbed abruptly. 
This disturbance is transmitted through the 
entire system. Also, no doubt, due to the action 
of the pumps there is a certain amount of re- 
mixing at the ends of the columns, a factor which 
would cut down the effective over-all length. 
In the previous work it was shown that in section 
1 there was some remixing due to some asym- 
metry. It might well be true that this is different 
for 2 and 3 than for 1. The existing line voltage 
fluctuation of 4 percent would also contribute to 
the remixing. No account has been taken in the 
computations of the impurities. 

This apparatus has been in continuous use 
since the fall of 1940 and has required practically 
no attention aside from time spent in removing 
the samples from the bottom and replenishing the 
gas at the top. It would thus be fair to state that 
including what little maintenance there is, not over 
13 hours per week of a man’s time was needed. 

Starting from scratch, it would probably be 
easier to build up the column by connecting 
standard sections in series and parallel. The 
third section which is 36 feet long may be a con- 
venient choice. With two of these sections in 
parallel feeding a third in series, we compute a 
multiplication of 11.2 times for the same rate of 
drain as used above. With three sections in 
parallel feeding one in series the factor is 14 
times. Probably the values obtained in practice 
would be somewhat less than these theoretical 
values as was true in the case investigated: 

In the design of any new column one would, 
of course, choose the heating wires of such a size 
that an autotransformer would not be necessary 
as was the case for sections 1 and 2 of our set-up. 

The authors are indebted to Mr. Rudolph 
Thorness for his skill in constructing and as- 
sembling the columns. 
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The Heat Capacities of Molecular Lattices 


I. Introduction 


R. C. Lorp, Jr. 
Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 


(Received June 5, 1941) 


The interpretation of the heat capacities of molecular crystals is discussed. The assumption 
is made that the lattice part of the molecular heat capacity can be represented by a single 
Debye function with the appropriate number of degrees of freedom. The reasons for supposing 
that this assumption is valid are considered, and it is shown that one can reasonably expect 
its validity except at temperatures low in comparison with the characteristic temperature @ of 
the Debye function. It is proposed to apply this point of view to the interpretation of the heat 
capacity curves of a number of molecules. The conversion of C, to C, is discussed for molecular 


lattices. 





T has been widely recognized that the theory 

of specific heats developed by Debye for 
atomic lattices has its shortcomings. Nonetheless 
it has been attended with much success from an 
empirical point of view. The number of examples 
of atomic specific heats which depart markedly 
from the Debye curve because of the inexactness 
of the Debye theory is small. Moreover, the 
development of a more exact theory has been 
precluded up to the present by the inherent 
difficulties of the problem, and it seems hardly 
probable that an improved theory is going to 
better the concordance between theory and 
experiment to any great extent. 

In the last twenty years or so a considerabie 
body of accurate heat capacity measurements on 
molecular crystals has accumulated. It appears 
quite likely that a satisfactory analysis of these 
measurements will yield valuable information 
about molecular structures in the same fashion 
that analysis of atomic heat capacities has led, 
in the past, to increased knowledge of atomic and 
lattice structure. Analysis of molecular heat 
capacities in the solid state, of course, must rest 
on some theory of lattice heat capacity. Pending 
the development of a tractable theory to supplant 
Debye’s for atomic lattices, it seems desirable to 
carry out this analysis as far as possible on the 
basis of the Debye theory. It is the purpose of 
the present series of papers to make such an 
analysis of the heat capacity data for a number 
of molecular crystals, and to point out examples 
of the sort of information on molecular and 
lattice structure which can accrue from study of 
heat capacities in the solid state. 


BRIEF SURVEY OF HEAT CAPACITY THEORY 


A crystal composed of N atoms has 3N degrees 
of freedom. If we suppose the interatomic forces 
to be purely harmonic ones, that is to say, 
Hooke’s law forces, all 3N degrees of freedom! 
are associated with harmonic vibrations. If we 
know the interatomic force constants and the 
atomic masses we can then, in principle, calculate 
the 3N normal coordinates and their vibrational 
frequencies. Once these latter are known, the 
heat capacity of the lattice at constant volume 
C, can be immediately calculated from the 
expression 

aN x; exp (xj) 
C.=k>d , (1) 
i=1 (exp (x;)— 1)? 
where k is Boltzmann’s constant and x,=hyv;/kT. 
The »; are the various frequency values. 

The problem of finding C, thus reduces to that 
of finding the »;. The simplest assumption, made 
by Einstein in his original theory of C,, is that 
all the v; are equal, in which case 


C,=3R[x%e7/(e7—1)?]. (2) 


The C,-temperature curve given by (2) deviates 
from observation at low temperatures, but agrees 
not too badly, considering the crudity of the 
assumption, with the general form of the ob- 
served curve. Much better agreement is found if 
one makes the assumption of Debye that the 
number of degrees of freedom N; having fre- 
quency »y; is given by 

N;=ar?’, (3) 





‘We must except the six degrees of translational and 
rotational freedom of the crystal as a whole, but these six 
are negligible in comparison with the other 3N—6. 
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a being a proportionality constant. This assump- 
tion is based on an analogy of the vibrations of a 
discrete lattice with those of an isotropic con- 
tinuum. The principle difference between the two 
cases is that the continuum has an infinite, the 
lattice a finite (i.e., 3N), number of normal 
vibrations. The distribution of frequencies given 
by (3) is therefore cut off at that v value, ymax, 
which gives a total of 3N v values between zero 
and ymax. The C, due to this set of frequencies is 


3 max xtet*dx 
Cm 3x| — f ——| (4) 
0  (e*—1)? 


Xmax 


Detailed dynamical treatment of lattices in an 
attempt to find a rigorous distribution function 
to replace (3) which could be expressed in terms 
of the independently measurable elastic constants 
of the crystal has been carried out by Born? and 
more recently by Blackman.* Such studies show 
that the real distribution function may deviate 
widely from (3), but attempts to determine a 
more exact frequency distribution law which 
yields a workable C, equation have thus far been 
unsuccessful. One reason for this lack of success 
is the difficulty of solving the dynamical problem 
of finding the distribution law. A second is the 
unfortunate one that even when the distribution 
law is known, the integration of Eq. (1) offers 
further trouble. Because of these difficulties and 
of the empirical successes of the Debye theory, 
the Debye equation remains the only generally 
useful one for the interpretation of lattice C, in 
atomic crystals. 

The first step in the extension of the lattice 
theory of specific heats to molecular crystals is 
to suppose that the N atoms of the lattice are 
grouped into molecules of S atoms each in such 
a way that the extramolecular forces are much 
weaker than the intramolecular. As a result of 
this disparity, the 3N degrees of freedom can be 
classified conveniently into 6N/S extramolecular 
and (3S—6)N/S intramolecular degrees. Two 
assumptions are now made: (1) The frequencies 


2M. Born, “‘Atomtheorie des festen Zustandes,”” Encyl. 
Math. Wiss., Vol. 5, part 3. M. Born and M. Goeppert- 
Mayer, ‘‘Dynamische Gittertheorie der Kristalle,”” Hand- 
buch der Physik, Vol. 24, Part 2. 

3M. Blackman, Proc. Roy. Soc. A148, 365 (1934) et seq. 

‘For linear molecules these figures are respectively 
5N/S and (3S—5)N/S. Mutatis mutandis, everything said 
in the text for non-linear molecules applies equally well to 
linear ones. 
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corresponding to the (3S—6)N/S intramolecular 
degrees of freedom are (N/S)-fold degenerate ; 
(2) these (N/S)-fold degenerate frequencies have 
values quite close to the 3.S—6 frequencies of the 
S-atomic molecule in the gas phase. Under these 
assumptions that part of C, due to the intra- 
molecular frequencies is given by Eq. (1), in 
which the summation runs over the 3S—6 
frequencies, the entire sum being multiplied by 
N/S. When a knowledge of these frequencies is 
at hand, from spectroscopic sources for example, 
the intramolecular part of C, can be calculated 
without difficulty. 

To obtain the extramolecular part of C, (here- 
after called “lattice C,’’) we need the 6N/S 
frequency values associated with the extra- 
molecular degrees of freedom. Experience with 
atomic lattices has shown that we cannot expect 
to find these frequencies in any direct fashion. 
It is necessary to make assumptions about the 
apportionment of the 6N/S degrees of freedom 
among various frequency values. The distribu- 
tion schemes are naturally more varied for the 
molecular lattice than for the atomic: 

(1) One can make the assumption, analogous 
to that of Einstein for atomic solids, that all 
6N/S v values are the same. The lattice C, is 
then given by the simple Einstein formula 
(Eq. (2) with 6R in place of 3R). The agreement 
between calculated and observed temperature 
dependence of C, under this assumption is little 
better than that exhibited by the Einstein 
equation for atomic lattices. 

(2) The 6N/S degrees of freedom can be 
classified into 3N/S degrees associated with the 
motion of the center of gravity of the molecule 
and 3N/6 with angular displacements about the 
center of gravity. The former degrees of freedom 
(the “‘translational’””) may be assumed to have 
frequencies distributed in the Debye fashion and 
thus to have a C, given by the Debye function. 
The 3N/S angular degrees of freedom (the 
“‘rotational’’), on the other hand, may be sup- 
posed, as was first suggested by Born® to possess 
frequencies nearly alike. An Einstein function 
will represent, in that case, the C, due to the 
angular vibrations. To calculate the lattice C., 
then, two frequency values are needed: the max 


5M. Born, Dynamik der Kristallgitter (Teubner, Leipzig, 
1915). See particularly Eq. (200), p. 77. 
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Fic. 1. Temperature trend of apparent @ for molecular 
lattices. Rotational degrees of freedom assumed to follow 
a Debye curve. 


for the Debye distribution and the average », Prot, 
for the angular vibrations. The requisite informa- 
tion for the dynamical calculation of the 6 
values corresponding tO Ymax and Pro¢ is ordi- 
narily not available and there is no direct way 
of finding them independently of thermal meas- 
urements. A priori calculation of the lattice C, 
on the foregoing basis is, therefore, not possible. 
By trial and error, however, one can find those 
values Of @mnax and 6,.¢ which give the closest 
representation of C,. It is found, in practice, as 
shown, for instance, by Eucken and his co- 
workers,® that the representation so obtained is 
often very good. 

(3) A third classification assumes that the 
3N/S angular degrees of freedom have fre- 
quencies distributed in the Debye fashion.’ This 
distribution is known when the ymax is known, 
and the calculation of the lattice C, then in- 
volves a knowledge of two @max values. Since this 
information is no more obtainable on a priori 
grounds than were the @ values of the previous 
classification, one cannot do more than compare 
the lattice part of the observed C, with that 
calculated from two Debye functions using arbi- 
trary @ values. The representation of the ob- 
served C, in this function is just as good as that 
obtained with the assumptions of (2). 

The two preceding assumptions undoubtedly 
are Overly simple and are only a rough. approxi- 


° See, for example, A. Eucken and E. Schréder, Zeits. f. 
physik. Cher. . B41, 307 (1938). 

’ This classification was first proposed by D. H. Andrews, 
Proc. K. Akad. Amst. 29, 744 (1926). 
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mation to the actual frequency distribution for 
the extramolecular degrees of freedom. That 
both of them reproduce the lattice C, rather 
accurately is attributable to two causes: (a) Even 
if a frequency distribution function is only 
approximately correct, the C, calculated from 
the distribution may often be very close to the 
true C,.8 (b) The liberty of choice of two param- 
eters (the two 6’s) clearly makes possible the 
better matching of calculation with observation. 

The usefulness of crystalline heat capacity 
data in aiding the solution of molecular structural 
problems is enhanced if we can find a simple 
and satisfactory way of allowing for the lattice 
C,. It seems to the writer preferable to adopt 
for this purpose the assumption (3), with the 
further restriction that the 6,,.x for the rotational 
degrees of freedom be assumed equal to the @nax 
for the translational ones.’ Such a simplification 
reduces to a minimum the number of parameters 
in terms of which the lattice C, is expressed. 
When the @ values for the translational and 
rotational degrees of freedom are equal, the 6N/S 
extramolecular frequencies give rise to a specific 
heat expressed by the usual Debye equation (4) 
except that the 3R of that equation is replaced 
by 6R. 

At first sight it may seem completely un- 
justifiable to make such a restrictive a priori 
assumption about the relative sizes of the @nax. 
It is shown below, however, that as long as the 
two 6’s are not drastically different, the lattice 
C, is rather accurately representable by a Debye 
function in which a kind of mean @Omax is used. 
The two values in many actual molecular crystals 
turn out empirically to be close together. This is 
not surprising when we consider that molecules 
of relatively high molecular weight, on which 
Omax for the translational degrees of freedom 


8 This fact probably accounts for much of the success of 
the Debye theory. It is interesting to note that those 
authors who wish to emphasize the virtues of the Debye 
theory do so by comparing calculated and observed C, 
(e.g., E. Schrédinger, Handbuch der Physik, Vol. X, p. 
305 ff.). Those who are more concerned with its short- 
comings (Blackman, reference 3; R. H. Fowler, Statistical 
Mechanics (Cambridge University Press, 1936), second 
edition, p. 131; W. H. Keesom and C. W. Clark, Physica 
2, 698 (1935), etc.) compare @ values obtained from C, 
at various temperatures with the theoretically constant @ 
demanded by Debye theory. ’ 

® This additional restriction was first made by Andrews 
(reference 7), in his original treatment of C, for crystalline 
benzene. It is particularly valid for that crystal. 
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depends, ordinarily have correspondingly high 
moments of inertia, which determine the rota- 
tional Onax. The most important exception to this 
generalization is the class of molecules whose 
molecular weights are made up mostly of atomic 
weights much heavier than hydrogen but whose 
moments of inertia are determined chiefly by 
hydrogen atoms. The hydrogen halides, water, 
ammonia and methane are good examples of the 
class. For them, however, our fundamental 
assumption of harmonic extramolecular forces is 
very bad and we could scarcely expect to under- 
stand their heat capacities on the basis of any 
theory restricted to harmonic vibrations. Ac- 
cordingly, we shall make no attempt, in the 
present group of articles, to consider any mo- 
lecular crystal whose C, depends on_ higher 
terms than the quadratic in the power series 
expressing intermolecular potential as a function 
of intermolecular distance. 

The accuracy with which the lattice C, of a 
molecular crystal can be represented by a six- 
degree Debye function involving a single @nax 
may be shown as follows. Suppose, for a given 
crystal, the @nax for the translational and rota- 
tional degrees of freedom are equal. If the 
lattice C, of this crystal be utilized to find the 
value of the Debye 6 at various temperatures, 
6 will, of course, be found to have a fixed value 
independent of T. A plot of 6 versus T will show 
a horizontal straight line. If the two @mnax have 
slightly different values, however, the resultant 
lattice C,, when utilized to find 6 as a function 
of 7, will give a 6 value which is not constant but 
which has a temperature-dependent trend. The 
nature of this trend is illustrated in Fig. 1. 
The ordinate gives the value of the apparent 6 
obtained at any temperature by substituting the 
lattice C, in the Debye equation. The lattice C, 
at any TJ is found by adding the C, due to the 
translational degrees of freedom and the C, due 
to the rotational degrees. These C, values are 
found from the Debye equation in which assumed 
values of @max for the respective degrees of 
freedom have been used. 

For example, if a hypothetical molecular 
crystal has a O0max for translation of 150° and 
Omax for rotation of 300°, it will have a lattice C, 
given by the sum of two Debye equations 
(Eq. (4) ] in which the respective 6’s have been 


LORD, JR. 
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Fic. 2. Temperature trend of apparent @ for molecular 
lattices. Rotational degrees of freedom assumed to follow 
an Einstein curve. 


used. This lattice C, may be used in conjunction 
with a six-degree Debye equation [Eq. (4) in 
which 6R replaces 3R | to calculate an apparent @. 
The apparent @ will be somewhere between the 
two original @ values at all temperatures, but 
will exhibit the above-mentioned trend. What we 
want to investigate is the quantitative character 
of this trend for hypothetical values of the @’s 
in order that we may interpret the trends ob- 
served in the apparent @’s in actual molecular 
crystals. 

In order to make the curves in Fig. 1 some- 
what more general, the @ values are plotted as 
fractions of the constant 6 value (referred to as 
6.), which they ultimately reach at high tem- 
peratures. The temperature scale is correspond- 
ingly generalized by plotting temperature as a 
fraction of @,.. This generalization enables the 
comparison of the temperature dependence of 
the apparent @ for various real crystals, regardless 
of whether the actual 6 values are large or small. 
That the apparent 6’s do approach a limiting 
high temperature value can readily be shown by 
using the form for the Debye equation at high 
temperatures: 

1 @ 
c=3Rr(1-——), 5) 
20 T? 


By adding together the C, terms given by two 
of these equations, each with its own 86, and 
equating the sum to a six-degree function [ Eq. 
(5) in which 6R replaces 3R], the limiting high 
temperature @ value in the latter may be found. 
It is 


6-=[(r?+1)/2 }*Otrans, (6) 
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in which @rans is the @ value assumed for the 
translational degrees of freedomand 7 = 6yot/ trans * 
Aor 1S the other @ value. 

In plotting the curves of Fig. 1, various ratios 
of translational to rotational @ values have been 
used. These ratios are indicated on the respective 
curves. Thus the curve marked “1 : 2” 
sents the apparent @ trend for a crystal whose 
respective translational and rotational 6’s are, 
say, 100 and 200. The @, for these two values is 
somewhere in the neighborhood of 150, to be 
exact at 158. To obtain the ordinate points, then, 
every apparent @ was divided by 158, and the 
temperature scale was altered similarly by di- 
viding each temperature by 158. In this way a 
curve is obtained which is independent of the 
initial @ values; precisely the same curve would 
have been found if the initial @’s had been chosen, 
for example, as 500 and 1000. 

The curves of Fig. 1 show clearly that a single 
§ value (corresponding to @,.) can be used to 
represent the lattice C, curve with accuracy 
comparable to that of experiment at tempera- 
tures above 0.3—0.46., even though the transla- 
tional and rotational @ values may differ con- 
siderably. Of course, when the two 6’s are not 
far apart, the representation is accurate to much 
lower temperatures. We can, therefore, express 
the lattice C, in terms of a single @ with sufficient 
accuracy in that temperature range where the 
intramolecular C, usually becomes of interest, 
namely from 0.3—-0.46, up to a few degrees short 
of the melting point. In the latter neighborhood 
the anharmonicity of the lattice frequencies be- 
comes important and no theory based solely on 
harmonic vibrations holds good. 

In addition to calculating the @ vs. T curves on 
the basis of Debye theory (Fig. 1), curves were 
calculated also on the assumption that the rota- 
tional degrees of extramolecular freedom had 
only one frequency, i.e., that the rotational 
degrees contributed to C, according to the Ein- 
stein equation. These curves are given in Fig. 2. 
In order to calculate the limiting value of 6 
when the rotational C, follows an Einstein 
curve, it is, of course, necessary to use the high 
temperature Einstein equation analogous to Eq. 
(5). It is 


repre- 


(7) 
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This equation, coupled with the Debye equation 
for the translational degrees of freedom, leads to 


0.=[(Sr?+3)/6 ]!Orrans- (8) 


The general form of the curves in Fig. 2 is not 
much different from that of the curves in Fig. 1 
so long as we consider only 6r:/@trans ratios 


larger than 1.14. When the ratio is smaller, the 
apparent @ curve rises above the @, asymptote 
at sufficiently low temperatures. As can be seen 
from Fig. 1, such a rise does not occur, regardless 
of the relative values of the 6's, when both 
translational and rotational degrees of freedom 
follow a Debye frequency distribution. Figure 2 
shows that when the rotational @ is well below 
that for translation, a pronounced minimum 
develops in the apparent @ curve. As one pro- 
ceeds to sufficiently low temperatures, however, 
the Debye heat capacity will ultimately become 
larger than the Einstein, and the apparent 6 
curve will, in the end, always rise above the @, 
value. 

The close similarity between the curves in 
Fig. 1 and those of Fig. 2 whose ratios are less 
than 1: 1.14 means that an apparent @ vs. T 
plot for a real crystal will not make clear whether 
the rotational degrees of freedom give rise to a 
Debye or an Einstein specific heat contribution. 
Of course, if, for a real crystal, the Einstein @ 
were smaller than the Debye @ and heat capacity 
data were available to sufficiently low tempera- 
tures, one might be able to decide by the upturn 
in the curve at low temperatures that the rota- 
tional degrees have Einstein heat capacity. It 
does not, however, seem probable physically 
that the rotational 6 should be much, if at all, 
lower than the translational 6. Indeed, in the 
molecular crystals so far studied, the curves 
follow the form of those in Fig. 1. The only con- 
clusion to be drawn from Fig. 2 for our purpose 
is that the lattice C, can be represented at tem- 
peratures above 0.3-0.4 6, by a single six-degree 
Debye function, even when the rotational degrees 
of freedom contribute to the lattice C, according 
to the Einstein equation. 

Curves analogous to those of Figs. 1 and 2 can 
be drawn for molecular crystals whose molecules 
have only two degrees of rotational freedom. 
These curves differ not at all qualitatively, and 
very little quantitatively, from those of Figs. 1 











and 2. It can be easily shown that the equations 
for linear molecules analogous to (6) and (8) are, 
respectively : 


0.= [ (27r?+3)/5 }*Otrans (9) 
and 
@.=[(1072+9)/15}Otrans- (10) 


In the papers following this Introduction, the 
heat capacities of a number of molecular crystals 
will be interpreted by assuming that the lattice 
part of the C, can be represented by a Debye 
function with a single 6, and the appropriate 
number of degrees of freedom. When the lattice 
C, has been subtracted from the total, the 
remainder will be used as a source of information 
concerning other degrees of freedom of the 
molecule. The kinds of information obtainable 
are varied and of considerable interest in con- 
nection with questions of molecular structure. 
For instance, the location with fair accuracy of 
low molecular frequencies which are, for one 
reason or another, spectroscopically inaccessible 
can often be accomplished. The structure of 
molecules in the solid state can sometimes be 
determined—for example, the linearity or non- 
linearity of a molecule can be established. It 
may also prove possible to determine the energy 
of a low-lying electronic level, or the frequency 
of a molecular vibration associated with re- 
stricted rotation within the molecule. 

The paper following this one will deal with a 
case in which the temperature trend of the Debye 
6 for a presumably atomic lattice can be reduced 
considerably by assuming an association of atoms 
into diatomic molecules. We will then consider a 
number of simple molecules for which both the 
specific heat!® and the molecular energy levels 
are accurately known, to illustrate the accuracy 
with which C, can be calculated when all the 
molecular information is available. Among these 
molecules are Cl2, CO2, benzene and benzene-dg. 
Following this demonstration of the applica- 
bility of the method, consideration will be given 
to certain molecules whose heat capacities yield 
hitherto unavailable information. 


0 Tt has been our experience in these calculations that 
accurate specific heat data are quite essential. Fortunately, 
the accuracy of much of present-day thermal data on 
molecules, particularly those studied by Giauque and his 
collaborators, is ample for our purposes. 
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THE DETERMINATION OF C, FROM C, 


The greatest difficulty in the utilization of heat 
capacity data for molecular crystals arises in the 
conversion of experimentally measured C, to C,. 
This conversion can be accomplished at any 
temperature by means of the thermodynamic 
formula 


C,—C,=ae°VT/B. (11) 


Use of Eq. (11), however, requires a knowledge 
of a, the coefficient of thermal expansion, 8, the 
coefficient of compressibility, and V, the molar 
volume, all of which are functions of temperature. 
The necessary information is often available for 
atomic lattices!! but rarely is for molecular 
crystals. Some modification is, therefore, re- 
quired of Eq. (11) which will enable one to use 
it in the absence of complete knowledge of the 7 
dependence of a, 8 and V. 

The important factor in the right-hand mem- 
ber of Eq. (11) is a’, since the temperature de- 
pendences of both 6 and V are small and to a 
large extent compensatory. The thermal expan- 
sion of a crystal whose atoms are held in their 
equilibrium positions by purely harmonic forces 
is zero. The expansion observed in actual crystals 
is due to the anharmonic character of these 
forces. Griineisen,’ has developed an equation 
for the temperature dependence of a which 
shows this explicitly. Extension of Griineisen’s 
reasoning to molecular lattices can be made in 
the following way. 

It can be shown™ by thermodynamic argu- 
ment, combined with the theory of quasi-har- 
monic lattice vibrations, that the thermal expan- 
sion of a crystal is given by 


B 3N 0; 0 ln 6; ? 
HE n(2)(M2), ca 
V i= T/ \a ln V 


where fg is the Einstein function given in Eq. (1) 
and the summation extends, as it does in Eq. (1), 


1 See, for example, R. C. Lord, Jr., J. Chem. Phys. 9, 
700 (1941) (this issue). , 7 

12 E. Griineisen, Handbuch der Physik, Vol. X, p. 22 ft- 

13 Cf. Griineisen (reference 12); see also A. Eucken, 
Handbuch der Experimentalphysik, Vol. VIII, part 1, p. 281. 
The discussion here given of the C,—C, formula tor 
molecular crystals was indicated in outline in Lord, 
Andrews, and Ahlberg, J. Chem. Phys. 5, 651 (1937). A 
lucid discussion of thermal expansion and the equation 0! 
state of solids will be found in J. C. Slater, Introduction to 
Chemical Physics (McGraw-Hill, 1939), Chapter XIII 
Slater derives Eq. (12) (p. 219, Eq. (4.13)). 
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over the 3N vibrational degrees of freedom of 
the crystal. The factor (@ In 6;)/(@ In V) gives 
the dependence of the vibrational frequencies on 
the volume of the crystal. For strictly harmonic 
vibrations this factor is zero. For quasi-harmonic 
vibrations it is an index to the relationship be- 
tween the amplitude of the atomic vibrations and 
the volume of the crystal. The fx-factor gives the 
temperature-rate at which the crystal takes up 
energy, i.e., increases amplitude of vibration. The 
product of these two factors is then an index of 
the temperature-rate of volume change. As we 
have mentioned above, the 7 dependence of 8 
and V is not great, and we will regard their 


quotient as a 7-independent constant. We will, 


in addition, assume that each (0 In 6;)/(d In V) 
isa 7-independent constant. This approximation 
is good as long as the anharmonicity is small. 

To find a workable expression for a, (12) can 
be simplified with the help of our discussion of 
the classification of the 3N vibrations of a 
molecular crystal. The sum in (12) can be 
broken up in the following way: 


x (“)(: In =) NaS y ,. dO In & 
Fist dag on (-)(—-). 
m AT/Nainv]) sim \T/\ainv 


(13) 
The 6N/S extramolecular vibrations have been 
postulated to follow a Debye distribution. The 
implication of this distribution is that all the 
frequencies show the same volume dependence, 
i.e., that the factor (d In 6;)/(8 In V) is the same 
for all of the 6N/S frequencies. The factor can 
thus be removed from the summation sign, 
leaving a sum which is just the Debye C,. 

The right-hand member of (13) gives the 
contribution of the intramolecular vibrations to 
the thermal expansion. The sharper the separa- 
tion between these frequencies and those of the 
crystal lattice, the less likely it is that there will 
be any similarity in the values of (0 In @)/(d In V) 
for the intra- and extramolecular frequencies. 
Moreover, the intramolecular vibrations them- 
selves differ widely in character, and their anhar- 
monicities may be supposed to differ just as 
widely. For this reason no further simplification 
of Eq. (13) is possible without the help of addi- 
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tional assumptions. In order to effect this simpli- 
fication, we make the ad hoc assumption that all 
(3.S—6) values of (0 In 6,)/(0 In V) are the same. 
It is presumed, also, that they are temperature 
independent. When this common coefficient is 
removed from the summation, the remaining 
sum is simply the total internal C, of the 
molecule. 
Equation (12) thus reduces to: 


aV/B=a'C,(lattice)+6’C,(internal), (14) 


where 
a’=(0 In 6;)/(0 In V) and b’=(0 In 6,)/(A In V). 


Substituting for a in Eq. (11) and making our 
usual assumption of the constancy of 6/V, we 
find 


C,—C,=[aC,(lattice) +bC,(internal) #7. (15) 


If some independent way existed for finding the 
a and b constants, our conversion of C, to C, 
could be readily accomplished. Unfortunately, 
there appears, at present, to be no way of ob- 
taining them. Equation (15) can, therefore, be 
used only as an extrapolation formula. In case 
the a and B values are available in a limited 
temperature range, the a and b constants may be 
evaluated inside this range and C,—C, then 
calculated for temperatures outside it. If even 
this information is lacking, Eq. (15) can be 
used in a purely empirical fashion by evaluating 
a and b from the observed C,—C, differences in 
a temperature range over which the calculation 
of C, is free from doubt. C,—C, can then be 
calculated, for purposes of comparing observed 
and calculated C,, in temperature regions where 
the calculation of C, is dependent on assumptions. 

In the following papers, C,—C, will, of course, 
be evaluated by means of the thermodynamic 
Eq. (11) whenever the data are available. 
It will frequently happen, however, that purely 
empirical use of Eq. (15) must be made. Fortu- 
nately, the C,—C, difference is often so small 
that even rather large percentage error in it 
does not interfere with straightforward interpre- 
tation of C,. The only effect is a decrease in 
the closeness of the quantitative agreement 
between theory and experiment. 
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The failure of the Debye equation for atomic lattices to hold for the heat capacity of metallic 
lithium is attributed to the association of lithium atoms in the crystal into diatomic ions. It is 
then shown, with the help of a reasonable assumption as to the value of the diatomic vibration 
frequency, that the heat capacity of lithium can be understood without postulating a break- 
down of the Debye theory. Some implications of this point of view for the theory of the metallic 
properties of lithium are pointed out and further experimental work to decide the question 


proposed. 


ETALLIC lithium provides an interesting 

example of the helpfulness of specific heat 
studies in suggesting solutions to problems of 
structure in solids. In 1935 Simon and Swain! 
reported an accurate investigation of the heat 
capacity of lithium in the range 15°-300°K. 
They recognized that an anomaly of some sort 
existed in solid lithium, because the heat capacity 
failed to follow a Debye curve and did not have 
the usual simple relationship to the coefficient of 
thermal expansion.” Simon and Swain interpreted 
the observed specific heat as made up of two 
components. The first of these was the normal 
lattice specific heat, which was assumed to be 
given by the usual Debye equation (@=505°), 
and to which the thermal expansion bore its 
usual approximate relation of simple propor- 
tionality. The second was attributed to the 
presence, in the atoms composing the lattice, of 
a single electronic energy level which lay some 
140 cm- above the ground level. The sum of the 
contributions from the two sources agreed nicely 
with observed C, over the entire range of the 
measurements. 

This sort of interpretation had been used 
previously* to explain pronounced departures of 
metallic and other specific heats from the Debye 
curve and would have occasioned no surprise 
except for two difficulties. In the first place, the 
lithium atoms in the lattice should presumably 


1F. Simon and R. C. Swain, Zeits. f. physik. Chemie 
B28, 189 (1935). 

2F. Simon and R. Bergmann, Zeits. f. physik. Chemie 
B8, 255 (1930). 

3 Cf. the review by F. Simon, Ergeb. d. exakt. Naturwiss. 
9, 253 (1930). 


be in the form of Lit ions and the existence of a 
sharp electronic energy level of the value pos- 
tulated conflicts with our ideas of the structure 
of Lit. Secondly, the Debye @ needed to reproduce 
the lattice portion of the specific heat curve is 
surprisingly high in comparison with 6’s for the 
other alkali metals. It conflicts, moreover, with 
the @ value found for lithium from the tem- 
perature dependence of electrical resistance and 
from the temperature dependence of x-ray scat- 
tering (vide infra), the latter two methods giving 
about the same value. 

For these and other reasons the suggestion of 
Simon and Swain is probably to be rejected, 
despite the satisfactory quantitative agreement 
between theory and experiment to which it 
leads. It is rejected by Fowler,‘ for example, 
who cites lithium as an “excellent example of a 
simple substance with a C, which departs 
markedly from Debye’s curve” because of the 
incorrectness of Debye’s frequency-distribution 
function. If this be the real reason for the 
departure, we have no alternative but to accept 
the discrepancy without further comment until 
such time as a more satisfactory theory than 
Debye’s is proposed. The nature of the departure, 
however, as shown by the decided temperature 
trend in the Debye 8@,° gives us a clue to a possible 
means of reconciling the Debye theory with the 
C, curve for lithium. Of more importance, it 


4R. H. Fowler, Statistical Mechanics (Cambridge Un! 
versity Press, 1936), second edition, pp. 131-2. 

5 See Fowler, reference 4, p. 131, Table 10.1, or Simon 
and Swain, reference 1, p. 190, Table 1. 


700 





STRUCTURE 


leads us to an interesting conception of the 
structure of metallic lithium. 

Suppose that we were to assume, for lack of 
knowledge to the contrary, that a substance such 
as crystalline chlorine is an atomic solid. We 
should then expect to find that the lattice C, of 
the crystal is representable by the usual Debye 
equation for atomic lattices. If we substitute in 
this equation the observed C, values for solid 
chlorine at low temperatures, we find, however, 
that the @ given by the C, exhibits a marked 
temperature trend, increasing considerably with 
increasing temperature. If we have faith in the 
Debye equation for atomic lattices, we then 
infer that, for some reason, we cannot apply it 
to crystalline chlorine. If the reason be that solid 
chlorine is a molecular crystal, we can apply the 
Debye equation for molecular lattices as outlined 
in the previous paper. 

In actuality, of course, the evidence is that the 
chlorine atoms in crystalline chlorine are in the 
form of diatomic molecules. We should, there- 
fore, expect that the above temperature trend 
in @ would be eliminated if we remove from C, 
the component due to the internal degree of 
freedom and use the remainder to calculate 6 
from a 5-degree Debye equation. A small tem- 
perature trend in @ at very low temperatures 
should remain, for reasons discussed in the 
previous paper and illustrated there by Fig. 1. 

The temperature trends are plotted for mon- 
atomic and diatomic chlorine in Fig. 1. On the 
same axes is plotted the temperature trend for 
the atomic lithium @, 6, being set at 430. It will 
be seen at once that the trend in @ is similar to, 
but not so marked as, that found when crystalline 
chlorine was assumed to have an atomic lattice. 
If one infers from this similarity that lithium is 
not a monatomic solid but is diatomic, one 
should compute from C, the Debye @ for five 
degrees of freedom per gram molecular weight 
of Liz. We expect, then, to find a constant @, 
Save at very low temperatures and at high 
temperatures where the sixth, or internal, degree 
of freedom begins to contribute. The @—T curve 
for Lis is also shown in Fig. 1. 

The assumption of the existence of diatomic 
molecules in metallic lithium clearly carries with 
it serious implications with respect to the proper- 
ties of the metal. First of all, because lithium 7s 
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IZ: Lithium, atomic mode/ 
IZ: Chlorine, molecular model 


IZ: Lithigm, molecular model 











Fic. 1. Temperature trend of apparent 6@’s for lithium and 
chlorine crystals. 


a metal, the form of the diatomic molecule is 
hardly likely to be neutral Lie. We must presume 
that the molecules are ions, most probably with 
an electronic structure closely related to Lig*. 
The chemical binding in the latter in the gaseous 
state has been computed by accurate wave- 
mechanical methods‘ to be slightly stronger than 
that in neutral Lie. The free electrons necessary 
for the metallic properties of lithium are thus 
available, although the number of them is in 
principle reduced by a factor of two.’ 

A question which is of immediate importance 
for heat capacity theory is that of the numerical 
value of the internal vibration frequency which 
is the molecule’s sixth degree of freedom. Since 
the binding energy of List is only slightly 
greater than that of Lis, it may be supposed 
that the force constant for Li,*+ is not drastically 
different from that of Lis. The latter is known 
from band spectra. Making the two assumptions 
that the force constants are about the same and 
are not altered in passing from gas to solid, we 
have at our disposal an approximate value of the 
vibrational frequency of Li,+. The contribution 
of the internal degree of freedom to C, can thus 
be readily computed. 

The frequency of Liz has been observed’ as 

6H. M. James, J. Chem. Phys. 3, 9 (1935). 

7One might regard lithium, if the diatomic model is 
correct, as a crystal intermediate in character between the 
pure molecular, non-metallic lattice of solid Hs and the 
atomic and metallic lattice of Na. Cf. the interesting 
speculations of E. Wigner and H. B. Huntington on the 
possibility of a metallic modification of solid hydrogen at | 
extremely high pressures: J. Chem. Phys. 3, 764 (1935). 

8F. W. Loomis and R. E. Nusbaum, Phys. Rev. 38, 


1447 (1931); G. M. Almy and G. R. Irvin, Phys. Rev. 49, 
72 (1936). 
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TABLE I. Heat capacity of lithium. Diatomic model: 0= 345, 











v=420 cm. 
Co Co Cv oy Cp Cp 
T°K | (LATTICE) | (INTERNAL) | (TOTAL) | (OBS.) | (CALC.) | (OBs.) 
15 0.06 0.06 0.09 0.06 0.09 
20 0.15 | 0.15 0.19 0.15 0.19 
30 0.51 | 0.51 0.55 0.51 0.55 
40 1.13 1.13 1.15 1.43 i.35 
50 1.97 1.97 1.99 1.97 1.99 
60 2.86 0.01 2.87 2.84 2.88 2.86 
70 3.74 0.02 3.76 3.74 3.78 3.76 
80 4.55 0.06 4.61 4.62 4.64 4.64 
90 5.25 0.10 5.30 5.34 5.39 5.38 
100 5.85 0.17 6.02 6.04 6.08 6.10 
120 6.80 0.32 7.12 7.18 7.22 7.28 
140 7.48 0.49 7.97 8.02 8.13 8.16 
160 7.97 0.67 8.64 8.68 8.84 8.86 
180 8.32 0.82 9.14 9.14 9.40 9.36 
200 8.59 0.96 9.55 9.56 9.86 9.84 
220 8.81 1.09 9.90 9.88} 10.27 | 10.30 
240 8.98 1.19 10.17 | 10.18} 10.59 | 10.56 
260 9.10 1.28 10.38 | 10.44| 10.86 | 10.88 
280 9.21 1.36 10.57 | 10.62} 11.10 | 11.12 
300 9.30 1.43 10.73 | 10.78} 11.32 | 11.32 
320 9.38 1.48 10.86 11.50 
340 9.44 1.53 10.97 11.67 
360 9.49 1.57 11.06 11.81 
380 9.53 1.61 11.14 11.95 
400 9.57 1.64 11.21 12.07 
420 9.60 1.67 11.27 12.18 
440 9.63 1.70 11.33 12.30 


























350 cm, and by using this value in conjunction 
with a five-degree Debye function (@=355) to 
calculate C,, fairly satisfactory agreement with 
experiment is obtained at temperatures above 
30°K. The numerical agreement can be some- 
what improved if the frequency value is taken as 
420 cm rather than 350. One would scarcely 
expect the 350 value to carry over exactly, and 
the 420 cm figure probably represents a more 
nearly correct value for the Li.+ frequency in 
the solid.* For this reason we present in Table I 
the values of C, and C, for lithium calculated by 
setting 0p = 345 and >=420 cm. The calculated 
C, is compared with a so-called ‘‘observed’”’ C, 
obtained by converting the C, measurements 
with the help of the thermodynamic equation 
involving the coefficients of compression and 
thermal expansion. The latter of these quan- 
tities has been measured by Simon and Berg- 
mann’ over the entire temperature range in 
which C,—C, is appreciable, and the value of C,, 


® It is of interest to note that, in the analogous molecules 
H2 and H,* the chemical binding is much stronger in the 
former (4.48 vs. 2.65 ev) and the frequency decidedly higher 
(4160 vs. 2170 cm). 

10F, Simon and Bergmann, Zeits. f. physik. Chemie B8, 
255 (1930). 
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since it is computed from C, by a thermodynamic 
expression, may fairly be considered as ‘“‘ob- 
served.’’ The experimental C, values are listed 
beside a set of C,’s calculated from the expression 


Cyp=C,+aC,*T, (1) 


in which a, the empirical constant, is equal to 
1.7110 reciprocal calorie. The sort of agree- 
ment given by an approximate expression like 
Eq. (1) may be seen by comparing columns 6 
and 7, Table I. 

From their study of the thermal expansion, 
a, of lithium, Simon and Bergmann" found that 
no simple proportionality exists between C, and 
a. Indeed, it was this fact as much as the failure 
of the Debye law which led them to postulate the 
existence of an electronic level as a contributor 
to C,. Since such a level should presumably make 
no contribution to a, its share of the heat capacity 
must, perforce, be removed from C, before 
seeking a relationship between C, and a. Simon 
and Bergmann found an accurate proportionality 
between a and the residual heat capacity as 
given by a Debye curve with 6=505. 

For molecular lattices, however, the theory 
does not forecast so simple a_ relationship 
between C, and a. In fact, for a diatomic lattice 
the relation should be ™ 


a=kiCor+keCur, (2) 


where k; and ke are constants, and Coz and Cr; 
are the lattice and internal components of C,, 
respectively. Only in the event that ki=kz» will 
simple proportionality obtain, and, as we know 
from the C,—C, values of such a diatomic 
molecule as Clo, the two constants can have 
quite different values. 

If we assume the thermal expansion of lithium 
to be resolvable into two components, the two 
constants k; and k2 can be evaluated from Cvz, 
Cor and a@ at any two temperatures. It should 
then be possible to represent a as a function of 
temperature over the whole range of the crystal. 
The values of @ calculated in this way in the 
range 15-300°K are listed in Table II. It will be 
seen that from 100° to 260°, the range in which 
data are available, the agreement with experi- 
ment is fair. The values calculated by simple 

1 Cf. (a) Eq. (12) and (b) references 13 of previous paper 


(R. C. Lord, Jr., J. Chem. Phys. 9, 693 (1941) (this issue). 
2 Cf. Eq. (14) of reference 11. 
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proportion from a Debye @ of 515 ® are listed 
for comparison. Since the two methods give 
divergent results below 100°, it would be inter- 
esting to have at hand expansion data in that 
region. 


EVALUATION OF @ BY OTHER METHODS 


While the agreement is satisfactory between 
experimental values of heat capacity and of 
thermal expansion and those calculated on the 
basis of the diatomic model, something more than 
this numerical coincidence is needed to inspire 
confidence in the model. As we have mentioned 
previously, one of the reasons for doubting Simon 
and Swain’s original interpretation was the dis- 
agreement between the @ found by them and that 
obtained from a study of the temperature de- 
pendence of the electrical resistance of lithium. 
This study was made by Griineisen.“ He inves- 
tigated the numerical values of a parameter 6 
which appears in an equation portraying the 
temperature dependence of resistance. Since this 
equation rests on similar foundations and has a 
similar mathematical form to the Debye equa- 
tion, Griineisen expected and found experi- 
mentally for a number of metals, that the @ 
values obtained from electrical studies and from 
specific heats are very nearly the same. At the 
time of his work" the specific heat of lithium had 
TABLE II. Temperature dependence of thermal expansion. 
ki =2.9X10-°, k2=15 10-8 reciprocal calorie. 








T°K | @(LATTICE)} a(INTERNAL) 


0.17 

0.44 

1.47 

3.28 

5.71 

8.29 
10.85 
13.20 
15.23 
16.97 
19.72 
21.69 
23.11 
24.13 
24.91 
25.55 
26.04 
26.39 
26.71 
26.97 


a(CALC.) | a(OBs.) | a(CALc. @=515) 





0.09 
0.27 
0.91 
2.01 
3.79 
6.21 
9.27 
12.5 
15.9 
19.1 
25.1 
30.0 
34.0 
37.4 
40.0 
42.0 
43.7 
45.1 
46.2 


0.17 | 

0.44 

1.47 

3.28 

5.71 

8.44 
11.15 
14.10 
16.73 
19.52 
24.52 
29.04 
33.16 
36.43 
39.31 
41.90 
43.89 
45.59 
47.11 
28.42 




















: See reference 1, p. 191. 
* E. Griineisen, Ann. d. Physik 16, 530 (1933). 
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TABLE III. 6 values for lithium from various sources. 








SOURCE VALUE OF 6 





Temperature dependence of: 
Specific heat, monatomic model 
Specific heat, diatomic model 345 
Electrical conductivity 363 
X-ray intensities 352 


328-420 


Einstein formula 


Lindemann formula: 

1. Constant = 135 

2. Constant = 115 
Elastic constants (theor.) 
Elastic constants (est.) 





not been measured at low temperatures, and he 
was thus unable to compare a Debye @ with the 
6 value he obtained, namely 363. 

The disagreement between this figure and that 
of Simon and Bergmann? led Pankow and 
Scherrer to an investigation of the temperature 
dependence of the relative intensities of various 
lines in the x-ray pattern of lithium. Here again 
the theory involves a parameter @ in a fashion 
similar to the Debye theory, and opportunity is 
at hand to obtain a numerical value of 6. From 
intensity measurements at 90° and 293°K, 
Pankow and Scherrer, found a @ value of 352 
(+12). 


The agreement of these two 6 values with that found 
from specific heats with the help of the diatomic model 
makes us look into the various other means of calculating @. 
There are two equations of semi-empirical character, the 
Einstein'® and the Lindemann" which are, respectively, 

6=1.76X10-*B-IMIV4#, (3) 
6=135T,,i!MV-4, (4) 


in which 8 is the compressibility, M the atomic weight, 
V the atomic volume, and 7,, the melting point. These 
equations were intended to apply to atomic lattices and 
their extension to molecular lattices is a procedure of 
unknown validity. With regard to the diatomic model for 
lithium, there arises the question of what values to use 
for M and V. If we consider the diatomic molecule as a 
single mass unit, formally we should use twice the atomic 


% G. W. Pankow and P. Scherrer, Helv. Phys. Acta 7, 
644 (1934). G. W. Pankow, ibid. 9, 87 (1936). 

146A, Einstein, Ann. der Physik 34, 170 (1911). The 
constant in Einstein’s equation has been corrected by 
the factor 4/3 to make the equation yield a @ suitable 
for use with the Debye equation rather than with the 
original Einstein C, equation. See N. F. Mott and H. 
Jones, Properties of Metals and Alloys (Oxford University 
Press, 1936), p. 8. 

17F, A. Lindemann, Physik. Zeits. 11, 609 (1910). In- 
order to obtain agreement with the observed 6’s for sodium 
and potassium, the constant of the Lindemann equation 
needs to be reduced from 135 to 115. 
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weight for M. Since the intermolecular forces, as indicated 
by compressibility or melting temperature, are approxi- 
mately twice as great, however, on a diatomic molecule 
as on a single atom, the observed @ should be roughly the 
same whether or not the atoms are dimerized into mole- 
cules. It thus seems appropriate in applying these more or 
less empirical relationships to use atomic values for M 
and V for the diatomic model. 

Another procedure for finding @, based on the theory of 
elastic solids!® makes use of 8 and of Poisson’s ratio. 
Although Poisson’s ratio appears not to have been meas- 
ured for lithium or the other alkali metals, an approxi- 
mate calculation of @ may be made by estimating Poisson’s 
ratio as 3, a value which is probably not far from correct. 
The formula of Born leads to 6=336. Here, as in the 
Lindemann and Einstein equations, the atomic weight and 
volume of lithium have been used. 

An equivalent calculation of @ has been made by Fuchs'® 
in the course of a theoretical study of the elastic constants 
of lithium. He does not explicitly calculate Poisson’s 
ratio, but the method of obtaining @ is that of Born. 
Fuchs finds @=339. 


The @’s obtained by the several methods are 
assembled in Table III. It may readily be 
appreciated from this table that the various 
6 values hardly can be used to confirm or disprove 
the diatomic model. They do indicate, however, 
that @ is quite unlikely to have a value so high 
as 500. 


X-RAy STRUCTURAL STUDIES OF LITHIUM 


Until very recently all investigators of the 
structure of lithium by x-ray methods” agreed 
that the patterns obtained indicate clearly a 
body-centered structure, which has also been 
found for sodium and potassium. There was no 


TABLE IV. Low frequency limit of ultraviolet transmission. 








ATOMIC 





ELEMENT NUMBER vo(OBS.) vo(CALC.) 
Cs 55 6.8 X< 104 8.3 10" 
Rb 37 8.3 9.4 
K 19 9. 10.3 
Na 11 14.3 14.3 
Li 3 14.6 19.3 

13.7 (List) 








18 M. Born, ‘“‘Atomtheorie des festen Zustandes,’’ Encyc. 
Math. Wiss., Vol. 5, Part 3 (1923), pp. 642-43. Born’s 
discussion is summarized in Fowler (reference 4, p. 127), 
where the equations for the calculation of @ will be found. 

19K. Fuchs, Proc. Roy. Soc. A153, 622 (1936). 

2¢ A. W. Hull, Phys. Rev. 10, 661 (1917); J. M. Bijvoet 
and A. Karssen, Proc. K. Akad. Amst. 23, 1365 (1922); 29, 
1286 (1926); F. Simon and E. Vohsen, Zeits. f. Physik. 
Chemie A133, 165 (1928); P. Scherrer and B. Arakatzu, 
Helv. Phys. Acta 3, 48 (1930). See also reference 15. 
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difficulty with the interpretation of the x-ray 
spectra, and on the basis of this evidence alone 
the structure of lithium seems clearly to be 
atomic rather than molecular. Precise study”! of 
the intensities of x-rays scattered by metallic 
lithium powder at various temperatures shows, 
however, that there are complicated features of 
the x-ray spectrum which cannot be explained 
on the basis of a two-atom body-centered 
structure involving scattering centers with 
spherical symmetry. While a unique interpre- 
tation of the experimental data has not yet been 
demonstrated, it appears” that a 7,,*(c) struc- 
ture* which has 16 Li atoms per unit cell gives 
a better agreement with the x-ray intensities, 
which were measured by means of an ionization 
chamber. In this structure the lithium atoms are 
associated in pairs, a situation entirely com- 
patible with the diatomic model. A detailed 
analysis and discussion of these x-ray results is 
to be published by Professor Harris and Dr. 
Griffith in the near future. 


PROPERTIES OF LITHIUM METAL 


Apart from the foregoing thermal and x-ray 
evidence, there is additional experimental ma- 
terial which appears to support the diatomic 
model. One clear implication of the diatomic 
model is a reduction of the number of free elec- 
trons in the metal. If all the atoms are combined 
into List ions, the number of free electrons is 
reduced by a factor of 2. Thus any metallic 
property depending on the number of free 
electrons per unit volume would have an anom- 
alous value for lithium in comparison, let us say, 
with the other alkali metals, which presumably 
have their full complement of free electrons. 

An example of such a property is the trans- 
mission of ultraviolet radiation by films of the 
alkali metals.** There is some definite frequency 
in the ultraviolet region at which each of these 
metals, in thin films, becomes transparent to 


' a -" L. Griffith and P. M. Harris, Phys. Rev. 59, 110A 
1941). 
22 Private communication from Professor P. M. Harris. 
I am greatly indebted to Professor Harris for the privilege 
of hearing of these x-ray results in advance of their pub- 
lication and for interesting discussion of their correlation 
with the heat capacity data. , 
23 Int. Tab. zur Best. von Kristallstrukturen, Vol. 1, p. 333. 
* R. W. Wood, Physical Optics (Macmillan, New York, 
1934), third edition, pp. 558-560. 
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radiation. The frequency depends on the number 
of free electrons per cc.” We can calculate vo, the 
frequency at which transmission sets in, by 
means of the formula given by Zener 


vo? = Ne? / amo. (5) 
Here N is the number of electrons per cc, and e 
and mp are the electronic charge and mass. 
Comparison of calculated frequencies with those 
observed by Wood shows satisfactory agreement 
except for lithium (Table IV). When we use for 
N, however, a value corresponding to the 
diatomic model, agreement comparable to that 
found for the other metals is obtained. The 
theoretical and experimental values are plotted 
as functions of atomic number in Fig. 2. It will 
be seen that the plots follow one another nicely 
provided the diatomic model is assumed for 
lithium.*® 

There are various other properties of metals 
which are dependent on the number of free 
electrons in unit volume, but for none of these 
does comparison of theory with experiment 
appear to offer a clear-cut decision concerning 
the diatomic model. In general, the properties of 
lithium are out of line with those of the other 
alkalis but calculation based on the theory of 
metals is not precise enough to show what the 
source of the anomalies is. As good examples of 
these properties one can cite the energy widths 
of the soft x-ray emission bands,’ the paramag- 
netic susceptibility?* and the pressure coefficient 
of electrical resistance.?® The anomalous behavior 
of these properties of lithium in reference to the 


_ 3C. Zener, Nature 132, 968 (1934); N. F. Mott and C. 
Zener, Proc. Camb. Phil. Soc. 30, 249 (1934). See also A. H. 
Wilson, Theory of Metals (Cambridge University Press, 
1936), pp. 125-126. 

**It has been suggested by Wilson, (reference 25, p. 
126), that the discrepancies in Table IV might well be 
attributed to the effective mass. If we use an effective mass 
of 1.4 mp the agreement improves for Cs and Rb and very 
decidedly for Li. However, it gets worse for Na and K, and, 
more important, there is still no break in the value for Li 
as compared with Na. On the whole it seems to us that 
little improvement is made by using an effective mass in 
place of mo. 

7 See H. W. B. Skinner, Reports on Progress in Physics 
ee 1939), Vol. V, p. 257, and references there 
cited. 

** J. B. Sampson and F. Seitz, Phys. Rev. 58, 633 (1940). 

*? Experimental: P. W. Bridgman, Proc. Am. Acad. 72, 


Guat theoretical: N. H. Frank, Phys. Rev. 47, 282 
35). 


OF METALLIC 






















































LITHIUM 
20r- 
Q 
‘ 
7 r 
» 4 —e— Experiment 
‘ --©-- Theory 
\ . Theory, dialomic model 
a- > 
% Sr \ for lithium 
in - 
1o'* &** 
sec! 
10 
5 1 1 1 L 
° 10 zo 3° 40 kf) 


Atomic Number 


Fic. 2. Low frequency limit of transmission of ultraviolet 
radiation by the alkali metals. 


other alkalis has been explained separately and 
on various grounds for each of them, but it 
seems reasonable to suppose that if the theo- 
retical calculations could be repeated on the 
basis of the diatomic model, they would lead to 
the same sort of agreement with experiment as 
has been found for the other alkali metals. 


FURTHER EXPERIMENTAL WORK 


There are several experiments which might 
enable a decision of the question. One of these 
has already been suggested by Simon and Swain,*° 
namely the extension of the thermal expansion 
data to 15°K or lower. A comparision of the 
values listed in Table II, columns 4 and 6, with 
the experimental a’s in the low temperature 
range would be illuminating. It would also be 
helpful to have accurate C, data from 300°K to 
the melting point. At present, neither a nor C, 
is known very accurately in this range.*! 

In this connection it would also be valuable for 
comparative purposes to have a over the entire 
temperature range for which C, is known for 
sodium and potassium. Simon*® has suggested 
that the heat capacity data for these metals also 
can best be interpreted with the help of an 
electronic level. Here, however, the Debve curve 
seems to fit fairly well at low temperatures and 
the necessity for assuming other components of 
C, not at all convincing. 


30 Reference 1, p. 192. 
3t See, however, C. C. Bidwell, Phys. Rev. 28, 584 (1926). 
% F, Simon, Ber. Berlin. Akad. 33, 483 (1926). 
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In many solids the intermolecular forces are sufficiently 
short ranged so that practically the entire potential energy 
of the system results from interactions between nearest 
neighbors. The thermodynamic properties of a solid with 
respect to a given coordinate, a, (for example, in a ferro- 
magnetic system a might represent the excess unpaired 
electron spin at a given lattice point; or in a substitutional 
binary alloy a might denote which of the two possible 
kinds of atoms are at a given lattice point) can be found 
from the factor of the partition function which averages 
over all possible configurations of @ at all lattice points. 

Here the evaluation of such a factor of the partition 
function is reduced to the calculation of the largest char- 
acteristic value of a linear homogeneous operator equation 
involving the potential energy between two adjacent 
layers of lattice points in the solid. By assuming that all 
possible configurations of a layer are equally probable, 
a lower bound for the partition function for a simple cubic 


ie many physical systems containing a large 
number of molecules, the intermolecular 
forces are sufficiently short ranged to permit the 
total potential energy of the system to be repre- 
sented by the sum of the interactions of each 
molecule with its nearest neighbors. Much 
research has been devoted to applying the 
‘“‘nearest neighbor approximation”’ to the study 
of normal solids,! order-disorder phenomenon? in 
substituted binary alloys, hindered rotation,’ 
ferromagnetism‘ and related topics; but no 
general mathematical technique for the handling 
of these problems (which are abstractly alike) 
has been developed. Usually the partition func- 
tion of the system in question is not found 
directly, but is expressed as a function of a 
parameter which is so chosen to minimize the 
free energy of the system, and often correlations 
(which arise not through direct interaction but 
indirectly through cooperational effects of inter- 


* Sterling Research Fellow in Chemistry. 

1Cf. Born and von Karman, Physik. Zeits. 13, 297 
(1912). Blackman, Proc. Roy. Soc. 148, 365, 384 (1934); 
149, 117, 126 (1935). 

2 For review of this subject see Nix and Shockley, Rev. 
Mod. Phys. 10, 1 (1938). 

3 Cf. Kirkwood, J. Chem. Phys. 8, 205 (1940). 

‘Ising, Zeits. f. Physik 31, 253 (1925). 


lattice is found in terms of the partition function of a 
cube of eight molecules. 

The operator equation for a real crystal becomes quite 
complex, but in some problems it is possible to obtain 
characteristic vectors and characteristic values (and thus 
the partition function) at high and low temperatures. 
For intermediate temperatures a theory of resonance be- 
tween the low temperature ‘‘almost ordered” character- 
istic vectors and the high temperature ‘‘almost disordered” 
characteristic vectors is developed. In cases where the 
operator equation can be solved for microcrystals a few 
layers thick, a perturbation method is discussed in which 
interactions between the individual microcrystals are 
treated as perturbations. 

Some remarks are made concerning phase transitions 
and finally the general theory is applied to a superficial 
treatment of two-dimensional ferromagnetic plates. 


mediate lattice points) between distant neighbors 
are neglected. 

It is the purpose of this paper to introduce a 
method of reducing the evaluation of the par- 
tition function of nearest neighbor systems to the 
solution of linear homogeneous operator equa- 
tions with the hope that the vast store of 
knowledge concerning these equations will be 
helpful in analyzing some of the elusive problems 
of statistical mechanics. The general theory 
proposed in the first part of the paper will be 
applied in the second part to a superficial study 
of two-dimensional ferromagnetic materials. 


I. GENERAL THEORY 


1. Characteristic value theory of the partition 
function 


Suppose a system can be divided into LX 
XN cells as is done in Fig. 1, so that there will 
be L layers each composed of NX M cells® and 
each layer will have the same average behavior 
as any other layer. The molecules in each cell 
may have several or even a continuum of 
possible configurations. Since the same is true of 

5In a simple cubic lattice each cell would contain one 


molecule and each layer, MXN molecules. This is not 
necessarily so in more complicated lattices. 
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each layer, one can define a set {a;} of all those 
configurations available to a physical coordinate, 
a, of a layer. For example, in a binary substi- 
tutional alloy a; might correspond to an ith 
arrangement of the two kinds of atoms in a 
layer; or in a normal solid a; might be one of the 
continuum of simultaneous positions of all the 
molecules of a layer. 

If v(an) is the total potential energy with 
respect to a coordinate a between all neighboring 
cells of a layer in the mth configuration, and if 
V(Q@m, @n) is the total potential energy of inter- 
action between cells of one layer in thé mth 
configuration with their nearest neighbors of an 
adjacent layer in the mth configuration, the 
potential energy factor of the partition function 
with respect to the coordinate’ a is 


Z= 2) +++ exp —[v(ai) +--+ +0(az) 


{a1} {ap} 
+0(a1, a)+-+-+0(ar1, a1) \/kT. (1) 
Making the substitution 
V(ai, a;) =v(a;)/2+v(ai, aj) +v(a;) me (2) 
(1) becomes 
Soot 
Z= pe ak 2 {II exp [— V(ai, ais1)/kT ]} 
v1} ap} t=1 
Xexp (—[v(a1)+0(az) ]/2kT). (3) 


Since v(a1)+z(a@,)/2 is one-half the total poten- 
tial energy of the top and bottom layers of the 
system and can be neglected when the total 
number of layers is large, 


Z=D--- Lit] exp [— Vai, ai4s)/kT J}. (4) 
To simplify this expression let us expand each 
factor as a bilinear combination of the com- 


ponents of a set of linearly independent ortho- 
normal vectors {¥,}. If each cell has v possible 








exp [— V(a4, a) /RT | 
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Fic. 1. 


configurations, each layer has v”* possible con- 


figurations and each vector would have v™% 
components. Writing 
v,(1) 
Yr= ; 
vn(v™%) 
{*,= (y,(1), y,(2), rine le y(v"%)), (5) 
the orthonormality of the ’s demands that 
yMN 
‘tn tn = _ Wn (ai) Wn( ai) _ bua: (6) 
aj=l1 
The desired expansion is 
yMN ypMN 
exp — V(ai, a;)/kT= y & 7. AmnWm(ai)Wn(aj), (7) 
m=1 n=1 


where the coefficients, a,,, result from multi- 
plication of both sides of (7) by yx(ai)Yn(a,), 
summing over all elements of the sets {a;} and 
{a;}, and using the orthonormality conditions: 


yMN ,MN 


Qrn= LD DL vx(ai) exp L— V(ai, aj)/RT Wi(aj), 
aj=1 aj=1 (8a) 
or in matrix form 
t*;, exp (—V/kT)d, =Akhy (8b) 


with 


exp [— V(a4, a2) kT | 


exp (—V/kT) _ exp [— V(as, a1)/kT | exp [— V(aze, as) kT | 


\ 


One might now ask if it is possible to select the 
set {tn} so that a, is zero when h#k. When 
exp(—V/kT) is a symmetrical matrix the 
* Note the implicit assumption of independence of the 
various kinds of coordinates; for example, in a ferromag- 
netic system this would be equivalent to the assertion that 
Spin orientations are independent of lattice vibrations. 





answer is yes, for if the ¢’s are the characteristic 
vector of the matrix equation 


Ad=exp (—V/kT)e, (0) 
then, . 
u*, exp (—V/KT) dn =t* Ant, =An dn 
an 


Anz =AnGnk- 


(10) 
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X, Xs 


Fic. 2. 


It is with this symmetrical case that we shall 
be concerned. The diagonalized expansion of 
exp [— V(ai, a;)/kT |] becomes 


exp — V(ai, aj)/RT = Don AnWx(ai)We(aj). (11) 


' Substituting? (11) in (4) and applying the ortho- 
normality conditions 


Z='5 Me {XT ve(e)¥e(ax)} 
a {a1} {az} 
=> Me EE velon)|?. (12) 
k a= 


Schwarz’s inequality implies that 


»MN yMN 


|X velar) P= X velar) Z 1)=y™", (13) 


a= a\= a\;=1 


Equation (12) has reduced the evaluation of Z 
to solving a characteristic value problem, and, 
moreover, since in a real crystal the number of 
layers, L, is very large, one need only find the 


7An abstractly equivalent process has been used to 
evaluate some of the multiple integrals occurring in the 
theory of imperfect gases. See Montroll and Mayer, J. 
Chem. Phys., in press. 
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largest characteristic value Amax of (9). Thus* 
L—1 
Z= a 2 Wmax(a@1) ]?. (14) 
a) 


Physically | Wmax(ai)|? can be interpreted as the 
relative probability of any layer in the system 
being in the 7th configuration. It is apparent that 
at absolute zero the components of thmax Cor- 
responding to the completely ordered state have 
the value 1/4/o where mp» is the number of con- 
figurations corresponding to complete order—the 
other components are zero. At temperatures 
sufficiently high for thermal agitation to cause 
complete disorder, all configurations are equally 
probable and 


Cmax (1, a Hoey 1)/Jo™9. 


Before investigating the implications and 
technicalities of the results up to this point, let 
us linger a moment to note the formal changes 
required in those problems that permit the con- 
figuration of each cell (or layer) to achieve any 
one of a continuum of energy configurations. In 
the partition function all summation operations 
must be replaced by integrations 


z= | . -{ exp —[v(ai)+++-+v(ar-1, ax) |/kT 
a Xdtar}---dfar}. (1) 
d\a;} is the volume element in {a;} space; for 
example, if {a;} is the set of possible simul- 
taneous positions of the molecules in the ith 
layer, 
d {ai} = (dx;dy,d23) pas (dx uwdy wdoun)- 


All expansions of the potential energy exponen- 
tials must be expressed in terms of the charac- 
teristic functions of the integral equation: 


AV (ai) = ¥(a;) exp [— Vai, ai)/kT Wiajl., 


_ fai} : ) 
2. A lower bound to the partition function 


Since tmax is the vector which maximizes 
* exp (—V/kT)q/y*y=), (15) 
any arbitrary ¢ corresponds to a X that is =\max 
8Since log \21=O(LMN) and log »¥¥=O(MN), in 


taking logarithms to calculate physical properties one can 
neglect the second factor in (14). 
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and therefore would give a lower bound to Z. 
A which permits the calculation of a X=Amax 
without any special appeal to the problem con- 
cerned until the end of the calculation is 


to = (1, 1, 1, Me 1)//v™%, 


for in that case 


1 
\o=— > oF exp -— V(a, a2)/kT. 


pNM a1 a 


Now Aov™* is exactly the partition function of a 
two layer system in which the interactions 
between cells of the same layer is one-half that 
between those of different layers. By standing 
the system on end so that it becomes a system 
of M2XWN layers, the repetition of the ‘“‘com- 
pressing operation”’ (this time with a v*¥-dimen- 
sional vector) reduces the system to one of 
N2X2 layers. Finally, compression of the 
2X2XWN system in the N direction, one obtains 


L—1)(M—1)(N— 
Die = wl L-D(M—-))(N—-1) 


x p~3LMN+ 4(LM+ LN+MN)+4(1—-L—M-N) | (16) 
where wv‘ is the partition function of a system 
of 2X22 cells in which the interaction energy 
between a pair of adjacent cells is only one- 
fourth that in the original system. 

The evaluation of w can be simplified by 
picturing its 2X2 X2 cells (Fig..2) as a square 
with each numbered pair *—i—* in the cube 
replaced by an X; in the square. If o; is the 
configuration of the ith corner of the square, and 
if W(o;,0;) is the mutual potential energy 
between the ith corner of the square with its 
neighboring jth corner (remembering that all 
interactions are only one-fourth their value in 
the original LX MXN system) 


py? y2 
wt= Po --- exp —LW(o1, 02) + W(o2, o3) 
ai=1 o4=1 


(17) 


Developing exp — W(o:i, o;)/RT in terms of its v* 
orthonormal characteristic vectors $; (each with 
v* components), i.e., the characteristic vectors of 


+ W(as, a4) + Wo, o1)/kT ]}. 


wolor)= ¥ (es) exp [— Wor, 03)/kT] 


o2=1 
(o,=1, as te y?) 
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and substituting the resulting expressions in (17), 
one obtains 


(ui= o*;, exp [ —W/kT ]o, 6* 40x). 


In special problems involving a specific knowledge 
of V(aj, a;), it is often possible to guess a better 
\ than Xo (usually it is desirable to choose one 
with an arbitrary parameter so that the vari- 
ational principle can be applied) and thus find 
a better lower bound for Z. This procedure is 
being applied by Professor Onsager to the cal- 
culation of residual entropies of crystals of polar 
molecules. 


3. Resonance between “almost ordered’’ and 
“almost disordered’”’ states 


As has been mentioned before, at a tem- 
perature of absolute zero where a system is 
completely ordered, its characteristic vector 
(unnormalized) has a unit component for every 
ordered configuration and a zero component for 
all others. Then Aax in the ordered state is 


Nor = eXp (— V(or)/kT), (18) 


where V(or) is the total V(a;, a;) in the ordered 
state. At temperatures above absolute zero it is 
possible by the method of successive approxima- 
tions to be described in Part II to express the 
non-ordered components of max as sums of 
exponential functions of (—1/kT) and to write 


Amax =exp [ — V (or)/kT ] 
X(1+>,a,(exp [—),/kT ]) ]. 


The characteristic vectors and values at low 
temperatures that are expressed in this manner 
will henceforth be written as ta. and Ag and 
will be said to represent the ‘‘almost ordered”’ 
state. 

At very high temperatures a system will 
approach complete disorder and all its con- 
figurations will have an equal probability of 
existence causing tmax—>(1, 1, «++, 1)/./v™%, and 
Amax—v™. Any decrease in temperature induces 
the components of ¢ to become linear com- 
binations® of inverse powers of kT and 


Amax = nd { +> Cn(1/kRT) *. 


(18’) 


(19) 


® This point will be discussed in more detail in Part II. 
























These high temperature characteristic vectors 
and characteristic values which will henceforth 
be represented by tla and Aaa will be called 
representations of the ‘‘almost disordered”’ state. 

In a given problem calculations from the ao 
and ad states are valid for quite large tem- 
perature ranges in the neighborhood of the 
extreme cases, but in the intermediate range 
(which usually includes that of a phase transition 
if there be one) neither the ao nor the ad vectors 
give exact results. However it is quite possible 
that a linear combination of these vectors, which 
physically is analogous to resonance between the 
ao and ad states (in the manner that resonance 
enters the theory of chemical binding and 
molecular structure!) gives a more exact charac- 
teristic value than either of the individual ones. 
More explicitly stated, there exist functions A(T) 
and B(T) such that 


t=A(T) aot B(T) haa. (20) 


By combining (20) and (15), the desired A and B 
are those which make 


A®*\got+2A Bai2+ Bra 
A*?+2ABb+B? 





=)\=>max, 
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where 
A12=*.4 Exp (—V/KT) tao 
= t*0 exp (—V/KT) tha, 
b= sathao- 
That is 
F(A, B) =A*\aot2A Baiot+ Brana 
—)(A?+2A Bb+ B*) =max. 


As is well known from calculus of variations, A 
and B must satisfy 


OF 
74 (Aso —A) + B(ai2— Ad) =(, 


_ OF 

2 ca ees =. 
But solutions to these equations exist only when 
d is a root of 





Aso — A wn 0 
Q\2—Ab Nad —A es 
The largest root of this equation is the charac- 


teristic value corresponding to the resonance 
state and has the value 


a Nad tAno ' 
Pts | [mafia 


— ' (21) 










If A and B are chosen so that ¢ is normalized 


and 


(1 —5?) 


A2+24 Bb+B?=1 









2 


_ (Aaa — A)? 
(a12—Nd)* — 2b(a12—Nb) (Naa —D) + (Maa A)?” (22) 
(Aao— A)? -- 
















The ratio A/(A+B) gives an indication of the 
degree of order of a system, and in the neighbor- 
hood of a phase transition one would expect it 
to vary rapidly with a small variation in tem- 
perature. 


4. Remarks on phase transitions 
As one raises the temperature of a physical 
system (containing a large number of molecules) 


Cf. Pauling, Nature of the Chemical Bond (Ithaca, 
1940), p. 8. 





~ (@1g—Nb)*— 26(a12— 0B) (Nao —A) + (Avo)? 


in a given state there often exist temperatures at 
which the system undergoes a radical change in 
state. For example, in a solid with a definite 
crystal structure the lattice is sharply disrupted 
at a temperature called the melting point, and 
above this temperature the liquid state exists. 
There is a temperature at which a ferromagnetic 
material sharply loses its magnetic properties. 
Such changes in state are called phase changes or 
phase transitions. At the transition temperature 
both states are equally probable and can coexist 
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at equilibrium in any proportion. These quali- 
tative properties of phase changes can be trans- 
lated into the language of characteristic value 
problems by ascribing a characteristic vector {; 
to the state below the transition temperature and 
another, ts, to that above it. Obviously \1>Az 
if 7<T trans and A2 >A If T > Trans, for otherwise 
the two would be interchanged (or one would not 
be a tmax at any temperature). Since both states 
are equally probable at Trans, the t's corre- 
sponding to both states must have the same 
characteristic value at 7 = T trans and the criterion 
for a sharp phase transition is 


Ai=A2 at T= T wens 
Ai >Aze if 7< E tune 
Ay <)o if T> T tzane: 


(i.e., Amax is doubly degenerate at Ttrans). 


When M and N are finite in a solid whose cells 
are limited to a finite number » of configurations, 

V(ai, a;)| is always finite and only as either M/ 
or N becomes infinite can | V(a;, a;)|—> «. There- 
fore the elements of the exp (—V/kT) matrix of 
a crystal with finite layers are all positive. Fro- 
benius!" proved that if all the elements of a finite 
matrix are >0, the largest characteristic value 
Amax iS positive and simple (i.e., nondegenerate). 
This theorem implies that Ajax of a finite crystal 
with a finite number of configurations per lattice 
point can at no finite temperature be degenerate, 


5. Perturbation theory of the partition function 
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and therefore sharp phase transitions are im- 
possible. If the crystal is infinite in one direction, 
this can be chosen as the ‘‘L’”’ direction and the 
elements of the exp (—V/kT) matrix would still 
be >0. However, if the crystal is infinite in two 
directions some of the matrix elements can be 
zero and a phase transition can occur. That is, a 
necessary condition for the existence of a sharp 
phase transition is that the crystal be infinite in 
at least two directions. By the same reasoning a 
necessary condition for the existence of a sharp 
phase transition in a two-dimensional square 
lattice is that it be infinite in both directions; and 
no sharp phase transition is possible in a one- 
dimensional nearest neighbor system, finite or 
infinite. 

Although real crystals are finite they have for 
all practical purposes an infinite number of 
lattice points per direction—or more exactly, the 
number of lattice points per direction is suf- 
ficiently large to produce configurations that 
cause some elements in exp (—V/kT) to be as 
close to zero as an experimenter can measure the 
sharpness of a transition. Thus our condition of 
necessity does not prevent the occurrence of a 
phase transition, but merely serves as a warning 
signal to be careful about expecting calculations 
on ‘‘small” crystals to give a correct picture of a 
phase transition, and to remember that the 
properties of a two- or three-dimensional system 
are not always deducible from the study of a one- 
dimensional one.” 


It is apparent that in a real crystal the order of the exp (—V/kT) matrix is so large that the direct 
calculation of its characteristic values is difficult if not impossible. However, one may be able to 
solve analytically (at least at high and low temperatures) the characteristic value problems associated 
with crystals that are only a few layers thick. This immediately suggests the use of a perturbation 
method for the solution of the real crystal problem, for physically a real crystal can be considered as 
a pile of microcrystals with the interactions between neighboring microcrystals as perturbations. Thus 
if the partition function of a crystal s layers thick can be evaluated, that of a real crystal (the subscript 
on Z in this discussion will denote the number of layers in the Z under consideration) can usefully 


be written 
[L/s] 


Z.=d mS >> II exp —[V(aci-1) 415 O(i—1) 04-2) F** * $V (@ie-1, ais) |/RT Xexp (—6(s)/kT), 


@] aL 


i=] 


where 


(23) 


O(s) = V(as, 541) + Vlas, @2e41) °° * + V(ar—s, @1—041)- 


"S. B. Frobenius, Preuss. Akad. Wiss., 514-518 (1909). 
See also R. Oldenburger, Duke Math. J. 6, 357 (1940). 
The author is indebted to Professor Onsager for pointing 
out this theorem and its implications. 


This has been observed before in special cases; for 
example, a one-dimensional gas does not have a sharp 
condensation point. Herzfeld and Goeppert-Mayer, J. 
Chem. Phys. 2, 28 (1934). 
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Using the high temperature expansion of exp — 0(s)/RkT 
[L/s] 


Zr=L'L TL (1—8(s)/kT+0%(s)/2(kT)?— -- -) 
' _— Xexp —[V(ages41, Q(i-1)s42) + om + V(ais—1, ais) | kT 
, — 4) 
tie (L/s—1) , ¢1js—2) 
Pa ag Pv lle 2 _— Z3(L/s—1)(V ot) 
Tr (Vs, s+1)0t ET)? —[Z5( )(Vs, s41)a 
+2(L, s—2)(V 8, 41 V 28,2841) wt (L, ‘s—1)(L, 1s ~ DV nccrdlndt poe (24) 
where 
k+1 


k 
Pr , , CF 
II ¥ «, ns+1 =he x II V oe Qns+1) exp [— | (Q(n—1)s+1 Q(n—1)842) + si +1 (Qns—1 Ons) | kl 
n=1 


Q@(k4+1)s n=1 


k 





= Pn grit +r, k+1 
=lim (—kT) > ee Se | 
— Op Our -OpPk a — aeyiys n=1 
»,-0 
X exp —[V(a(n—1) +15 Q(n—1)s41) + or + V (@ns—1, Ons) +uV (ans, @ns41) \/RT} . (25) 


(ui = 0) 
We shall not exhibit a general expression for these averages, but the method of calculating a par- 
ticular one will be apparent from a discussion of (Vs, 541). Here we must study 
Z2= a1 dade Dias exp — [ V(aya2)+ vr. + V(@s-1, as)+uV(as, 541) 
+ V(as41, Qs42) + -+++V(a2s-1, 2s) | RT, 
which is the partition function of the system in Fig. 3. Now Z> can equally well be calculated on 


the basis of the system being a pile of M layers of 2s N cells with an anemic interaction between 
the sth and the (s+1)st row in each layer: 


Z2= Dai: * Den exp { —[V(B1, Be) +--+ +V (Bus, Bu) J/kT}. 


V’(B1, B2) is the reduced double layer potential (of Eq. 2) between one layer in configuration 1 and the 


other in configuration 2. Proceeding in the manner discussed in Section 1: 
M-1, 


Z2= max t {die Wmax(B) } *, 


where max is the greatest characteristic value and Umax its corresponding characteristic vector of 


(exp [—V (61, B1)/RT] expl—V(B1, B2)/RT] ---) 
exp (—V,;/kT) = lexp [—V(B2, B1)/kT] exp [—V(Bs, B2)/kT] | 
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It is at this point that the perturbation theory is useful, for we can express our characteristic value 
problem as 

A¥(B1) = Lise (Be) exp {- L V o(B182) +H Vi(B1B2) J kT} (26a) 

AV(B1) = DiseW(B2) exp [ — Vo(8182)/kT ](1—wVi(B182)/RT+---). (26b) 
Vo(8182) is the reduced potential energy function for two adjacent layers in the two independent 
systems, and nV (8182) is the weak perturbation resulting from interactions between a pair of layers 
in the left system with the corresponding pair in the system on the right (Fig. 4). As is customary in 
perturbation theory we assume that 


A=AotwrAitwrAct:+-, t=totulityet+---, (27) 


where Ao and to satisfy the equation 


or 


Aotlo = Exp (—Vo kT). (28) 


Since we assumed that the s layer problem has already been solved, A» and wtp are known and Eq. 
(26b) becomes (neglecting all terms of 0(y?) and higher) : 


(No MA1) (ot wth) = exp (—Vo, ‘kKT)(1—nzV, RT) (Wot ut) 
otfo+ Aouryi + wArtyo= exp (— Vo/KT) (o— wVito/kT +p). 
Multiplying by ¢*) and remembering (28) 


Ay= — t*y exp (—Vo/KT)Vito/kT 


or 


ee = (29) 
= — Yar D642 Yo(B1) exp [— Vo(81B2)/RT ]V1(B1B2)Wo(B2)/ RT. 
Therefore, correct to terms of order yp: 
A=*oLexp (—Vo/KT)(1—uVi/kT) Jao. (30) 


For details necessary in the calculation of %, the reader is referred to the chapter on perturbation 
theory in any of the standard quantum mechanics texts. 

The higher order V averages are found by essentially the same method, and (V., 541) can be better 
approximated by taking second perturbations. 

In the low temperature region it seems advisable to handle the perturbations by expanding the 


exponential exp [ — @(s)/kT ] as 
[L/s] 


e~ 8(#) /kT = @—V (or) /kT (4 — > Sua nest oe -), 


n=1 


where V(or) is the reduced potential energy of the perturbing interactions between adjacent micro- 
crystals when they are in a.completely ordered state, and 


fij=exp | -| V (aij) -=v¢or)| AT —1. 


The development would then proceed by combining the high temperature formalism with that of 
Mayer’s treatment of imperfect gases." 


II. APPLICATION TO TWO-DIMENSIONAL FERROMAGNETIC NETS 
1. General remarks 


A problem well suited as an example for the application of the general theory of Part I is the 
Statistics of a two-dimensional ferromagnetic net. Here we shall not be so much concerned with the 
physics or complete solution of the problem as we shall be in using it as a medium of demonstration. 


** Mayer and Mayer, Statistical Mechanics (New York, 1940), p. 277. 
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Wannier and Kramers have studied ferromagnetic nets in detail by a method somewhat similar to 
that described here. They first show that the partition function is the largest characteristic value 
of a set of linear difference equations that are derived by considering the change in the configuration 
of the lattice points of the system by the addition of one more molecule. Their final characteristic 
value problem!" resembles ours rather closely. 

The seat of magnetic properties of a ferromagnetic crystal is the unpaired electron spins in the 
d shells of the atoms at each lattice point in the crystal.!® We shall let — € be the energy of interaction 
of two neighboring atoms whose unpaired spins are parallel and +. that when they are antiparallel. 
Also we shall write J=e/kT and in the operator equation’ exp (—V/kT) =A, the matrix elements 
will be abbreviated 

exp [— V(aia;)/kT ]=[ai, a]. 

In a system in which two possible configurations are available to each lattice point a “‘0’’ may be 
used to represent one configuration and a ‘‘1”’ the other. Then there exists a unique number 1 in the 
binary system (n=2" —1 if there are N lattice points per row) that corresponds to a given con- 
figuration of a row: 

a=0~ 0000: - -00 
a=1~0000---01 
a=2~0000---10 


a=2N%—1~1111---11 
Ferromagnetically speaking a ‘‘0” will represent a lattice point with a + excess spin and a ‘‘1”’ will 


represent a — excess spin. If v(o102) is the potential energy between a lattice point with excess spin 
o, and a nearest neighbor with excess spin o2 


v(1, 1)=2(0, 0), (1, 0) =0(0, 1). 
This symmetry implies the invariance of all potential energy functions with respect to interchanges 
of 0’s with 1’s. A configuration will be called complementary to another if the replacement of 0’s by 
1’s and vice versa makes the second configuration identical with the first. In the binary number 
system the exchange of 0’s and 1’s changes a number a to 2" —1—a because the sum of a number and 
complement is always 

1111---11=2"—1. 
The exp (—V/kT) matrix can be constructed in the symmetrical form 

exp (— V/kT) 
[0, 0)[0, 1] --» (0, 24-1-1] [0, 2¥—1] +++ [0, 24-1] 
(1, 0J[1, 1] -++ (1, 24-1-1] [1,2%¥—1] -++ (1, 24-1] 


[2%-1-1,0] «+--+. [2¥-t—4, 2-1-1} [2-1-1 2N-1] --- (2% 1, 29-1] 
(2*—1,0][2%-1,1] --- [2¥-1,2%-1] | [2"-1,2"-1]_ --- [2%—1, 29-1] 








| 
| 








(29-1, 0] ac aoe [2-1 oN... 1] | (24-1, i 1] i ow (24-1, 2N-1} ) 


4 Unfortunately the present state of world affairs has delayed publication of this work. However, it is abstracted in the 
Proceedings of the Boston Meeting of the American Physical Society, Feb. 1941. See Phys. Rev. 59, 683A (1941). The 
author is greatly indebted to Dr. Wannier for a discussion of his unpublished work; indeed to a large degree it was this 
discussion that inspired the writing of the present paper. in 

16 For reviews of research in this field see Bitter, Introduction to Ferromagnetism, (McGraw-Hill, New York, 1937), 
aa) V; and Van Vleck, The Theory of Electric and Magnetic Susceptibilities (Oxford University Press, 1932), Chap. XII. 

16 Note that in two-dimensional problems the word “‘layer” used so frequently in Part I is to be replaced by the word 
“row.” 
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with the upper left quadrant identical to the lower right and the upper right identical with the lower 
left. This quadrant symmetry follows from the observation that the element in the ith row and jth 
column in the upper left quadrant is [7, 7] and that in the lower right quadrant is [2” —1—1, 
2" —1—j ]. Since the configurations yielding these matrix elements are complementary, the elements 
are equal (the same sort of reasoning applies in the proof of upper right, lower left symmetry). Thus 


exp (—V/kT)=( °): 
[0, 07 --+ (0, 2¥-1-1] 


where 


POE ct kb Keke ke oe 4 
[24-10] +--+. [2%-1-1, 2%-1-1] 

(0, 2%-1] +» (0, 29-1] 

[1, 2-1] cos (1, 24-1) 

[2V-1~— a QN — 1] ai ieee [2V-1— t. 2N-1) 

[0,2%-1]  --- [0,29-"] 





(0, 25-1) ese [2V-1~— is 2N-1} 
The latter equality is the result of the property 
[i, 2¥—1—j7 J=[2" —1-—7, 2N—1—2N+7 ]=[2%-'!—-1 -3, 7]. 


Our operator Eq. (9) now becomes 


(5 Son 


Transforming to a new vector 6 by letting 


me he (1) 


=— )o=Ro; R=R" 
v2 


(1) (—1) 
-1 0 
0 0 -1 


- |? (-—1)= . os . 
J | 0 0 


‘P+Q 0 
( b=r6 
0 P-Q 
and the largest characteristic value of (33) is theflargest root of the secular equation 
|P+Q-—A| =0. (35) 


The (P—Q) portion of the reduced matrix yields smaller characteristic values and may be ignored ; 
so, if {;} is a set of orthonormal vectors in the reduced manifold of @ we need only study 


(P+Q)o=Aé. (36) 


The (P+Q) matrix is still reducible, but its reduction is not a very great simplification and will not 
be performed here. 


we have 
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2. Calculation of the partition function for some narrow strips* 


Using a strictly scalar model of spin interactions, we have the following partition function cal- 
culations. 

(a) Linear chains.—A linear ferromagnetic chain of N equally spaced lattice points has an easily 
calculated partition function. Its exp (—V/kT) matrix is 


ee of 
exp (— V/kT) = ( ) 
re U6 
and the secular equation of this operator shows the largest characteristic value to be Ajax =e” +e’. 
The corresponding characteristic vector is tmax=(1, 1)/v2. Therefore 
Z2n,1= 2" cosh*®—J. (37) 
(b) 2XWN strip.—Every spin configuration of a row of a 2X N ferromagnetic strip corresponds to 


a two-digit number in the binary system, and in the notation of Section 1 of this part, the secular 
equation for this problem is 


a 
[0, 0]+[0, 3]—A 2F0, 1] | 
e-7 +e —r 2 


” | =O. 
2 f+}! 


|P+Q-A| = 


Simple algebra shows the largest root of this equation to be 


ona 
Amax — a 


= 4J 


sinh J 


+a|- $[ (A?—16) coth J+A], 


A?= 16+ (e —e’ +e-/ —e—*”)?= 16+ (sinh 4//cosh J)? 
Ynax = (a, b, a, b)/v2 
a?=8/(A?—A[A?—16]!), 6?=8/(A?+A[A?—16]}). 
[Yea ¥(a) |? = (2a+2b)?/2=2(1+4/A) 


sinh4J A}~! 
Jats] +4/4). 


Zy,2=2 ae 
2sinhJ 2 


(c) 3XWN strips.—The secular equation of a 3X WN strip 
e +e —)d e* +-1 oe +-1 e+e 
ey +1 e*J + e—3J — ey +e e 7 +1 
1 +e2/ e! +e-/ est +e-3/ oun Xr 1 +e- 
ete 1+e-% 1te- el +e-J —y 


|P+Q-A|= 








(and that of wider strips) is of too great a degree to be solved in a direct analytical manner, so one 
must resort to numerical or approximation methods. Fortunately the method of successive approxi- 
mations is ideally suited to the high and low temperature evaluations of Amax- Here one guesses a 
reasonable vector $9 and calculates ‘ 2 

(P+Q)0=od1, (40) 


* The results of this section are not new, having been obtained by Wannier and Kramers; however, the approach, 
especially that in c—e, is new. 
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where the first component of 6; is factored into Aogo(1). This process is continued according to the 
general scheme 


(P+Q)oi=ibinr (6:1) =Fi41(1)), 


until the difference between 6; and 41 is as small as one desires. The largest characteristic value of 
(P+Q) to the same accuracy is” \;. In the completely ordered state the characteristic vector $max 
of (P+Q) is (0,0,0,1) so for low temperature (large —J values) calculations of \ we assume 


(41) 


$0= (0, 0, 0, 1) (note that no attention will be paid to normalization until the end of our calculations 
and that at each step higher powers of e+’ are neglected) and proceed by the method discussed above. 


e®J +e 
e* +1 
1+e* 
ey +e4 

e®4 + e-/ 
eo +1 
1+e7/ 


(P+Q) G0= 


ee +1 
e*7 + e-3/ 
ey’ +e-/ 
1 oe e —2J 
e* +1 
est +e-3/ 


e* +1 
ed +e/ 
eet +e-3/ 
i+e-* 
e* +1 
e-7+e/ 


e7 te~/ 
e~- +1 
i+e-~/ 
e+e 
e! aa e/ 
e241 
i+e-¥ 


(esd 
“ 
= e-5J | 
e3d 
1, 
etd 4 eb 
es + edd 


=(e- + 3et/)| |, 
| ee! + J2eJ 














e/ +e-/ 1 +e-/ 1 +e-2/ 


e+e) 1 ) 
until 


and 


Amax = "ed + 3e* ad + 7e37 +0(e>/ ) 
(ety +-e% + 3a 
1 e! + Veit 4 3¢7/ 
(1+ 2e* + Ve + 22¢')3| e394 + 2e5/ + 3e7/ , 





b20 _ 


1 


The characteristic vector of the P+(Q matrix in the disordered state (unnormalized), ¢* = (1, 1, 1, 1) 
will provide the first approximation at high temperatures. In this range it would be desirable to 
express the partition function in ascending powers of 1/kT (i.e., in powers of J) so the elements of 
the (P+Q) matrix will be written as power series in J. 


2+4J 24+2J 24+2/ 2 1 
2+2J 2 2 2—2J 
2+2J 2 2 2—2J| }1 
L 2 2—2J 2-23 2-43) \1 
2+4J 242F 242/ es +f 8 1 
2+2/ 2 2 2-27) | 1-J 1—J+0(J?) 
= (8+0(J?)) 

2+2/ 2 2 2-23) | i-J 1—J+0(J?) 
2 2-2J 2-2) 2-4/3) \1-—2/ 1—2J+0(J?)) 
2+4J+13J? 24+2J+2/J? 242/427? 2+J* 1 1 
2+2J+2S? 2+9J? 2+ J* 2—2J+2J* | i-J 1—J+J?/2 
24+2J+2/? 2+F* 2+9J? 


2—2F+2J* | 1—J 1—J+J2/2|’ 
243? 2—-2+4+2F? 2-2F4+27? 2-4F41372) 1-2, 1—2J+2J? 


(P+Q)do= =(8+8J) 


i-J 














=(8+14J") 


until 
i Amax ~ 8+ 2052+ 113J4/3+0(J*) 


_ For a detailed discussion of this method see Frazier, Duncan, and Collar, Elementary Matrices (University Press 
Cambridge, 1938), p. 140. 
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TABLE I. 




















aH 


kT 
Fic. 5. 


CUIKAmM EWE O 


and ¢max (normalized) is 
1 
_ 1-4 (J2/2) — 13/6. 2+ 4/24 
ne 1 — T+ (52/2) — 13/69 + J*/24 | [4 —8F+12F? —88/3F2+52I4}! 
1-2 +2J?—16/3F?-+26/3J 





In Table I the exact values of Amax as computed numerically are compared with those given by the 
vectors $ao and gaa. 


heo= Geo EXP (—V/KT)ba03 Aaa= dea xP (—V/KT) dua. 


(d) 2XWN strip by perturbation method.—Since the application of the perturbation method to wide 
strips requires much space, and since at the present we are only interested in principles, we shall be 
satisfied in discussing this method in connection with 2X WN strips. These strips can be imagined as 
two linear chains whose interaction is a perturbation. At high temperatures 


Z2= Dar Var exp { —[v(a1)+0(a2) ]/kT} exp [ —v(aia2)/kT ] 
= da > a2 exp 1 - [v(a1) +(e) | kT} (: ome ee — + * +) 
kT 2(kT)? 
=Z,2— (o(crcts))w , (o*(ara) dav a 
RT kT)? 
(v"(aia2) wv = >a >a v"( aya) exp { - [v(a1)+v(ae) |/kT} 


= am (—kT) "(0"/Ou")> a1 > 22 exp { a: [v(a1) +pv(aya2)+0(a2) |/kT} . 





where 


The quantity 
Z2'(u) = Dar Da: exp { —[v(ar) +0(a2) + pv(arae) |/kT} 


can equally well be written ; 
= Da De exp [— V(B1B2)/kT]=Xmax {de (8) }? 
with the @’s referring to the configuration and 
rye 
——6 @r 


Atv=exp (— V/RkT)¥. 
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Introduction of the transformation of Eq. (34) yields 
o=(P+Q)é, 
ed (e2 4+ ¢-2/) ? 
P+ -( ). 
( Q) 2 e-¥d (e2J 4 e-2/) 


In the unperturbed system, ».=0 and 


co aan 
2 e2t +e-2 


pa’ ‘yo 
“VIN -1/7 


tan corr 0 )o 


0 (e7 —e-/)? 


Now let 


Then 


The characteristic values and characteristic vectors are then 


ok kK 
M=(eo+e-/)2; 8 =(1,0) or $: =(1,1)/v2 
* 
A= (el —e-Y)?; O =(0,1) or (1, —1)/v2. 
ok ok 
Note that ‘ and . are the characteristic vectors and MN and Ae the characteristic values of the 


linear chain matrices, as they must be because we have a linear chain problem when p=0. 
The assumption that 


A=N 4AM +A ye eee, 6= 69+ 6 yt oO pe+ a 
leads to [see Eq. (29) ] 


J 2J —2J 0 1 
mi =4(1, 1)( linia )G) 0 
0 —J(ew~ te) I \1 


J (e747 +e) 0 1 
h2=41,D( )( )=se+e), 
0 —J(e%+e-7%) 7 \1 


1 . ;) 
esa 
v2\1—pJ/2 


at high temperatures. The normalized ¢;“ when »=1 can be written 


1 1 
rear 1) 


which checks with the corresponding vector as found by the method described in Section (c). Inas- 
much as A“ is independent of uy, 


(v(aya2) w=0 
and 


Z.=2P*+0(J%). 


Second perturbations give the term of 0(J?). 
(e) Introduction of external fields—Certain types of external fields do not invalidate the formal 
developments of Part I. For example, let us consider a linear ferromagnetic chain of N equally spaced 


** Cf. Pauling and Wilson, Introduction to Quantum Mechanics (McGraw-Hill, New York, 1935), p. 160. 
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lattice points in a magnetic field H whose direction is parallel to the chain. The interaction between Th 
field and spin when the chain is in a configuration with n, ‘‘+”’ spins and n_ ‘‘—”’ spins, is 
N 
~u(n,—n_)H= —pH & [ny (ai) —n_(ai)]. | - 


i=1 


“=magnetic moment at each lattice point. | 
n(a;)=number of ‘‘+”’ spins at the ith lattice point when it is in configuration a,. by 


n_(a;)=same with ‘‘—”’ spins. 


Then letting 





5 plT Fr 
V(aiaj) =v(aia;) “st [ns(ai)+n4(a;) ]—[n_(ai)+n_(a;) }}, of 
Zy,1(H) = Yai: ++ Deen exp | —[V(aiae) +--+ ++ V(ay-ray) J/kT} ] 
N-1 
= Amex { dia Wimax (a) } ‘, c 
The e~Y/*? matrix is 
ed HH /kT e7 
e-VIkT — ( ). At 
end el tuH/kT 
If Vinax = (a, db); a 2+5? =i, then a(el—#H/kT —y) 4 — =0, ant 
d ae~/ +d(elt#H/kT_)) =0 
an 
Amax =e” cosh pH /kT +(e” sinh? pH/kT+e-*’)}, , 
(Soa Wmax(a))?=a?+b?+ 2ab=1+2ab of 
=14+1/(1+e sinh pH/kT)!. - 


Therefore 


Zy,1(H) =[e’ cosh pH /kT +(e sinh? pH /kT+e-)? }8—"[14+1/(1+e” sinh pH/kT)?}]. 


int 


The magnetic moment of the chain is 
0 log Z 
~ Of Hu/kT] 
and that at saturation is J) =yN. Therefore 
I sinh pH/kT 
. (sinh? pH/ kT +e-*) 8 





This function is plotted in Fig. 5. 


3. Spin contributions to specific heat 


” d log Z c=(= ~) 
Ea 
d(—1/kT) . 
give the average energy and specific heat contribution of spin in terms of the partition function. They 


will be used in the following calculations. 
(a) Linear chain.— 


The well-known formulae 


dJ 
d(—1/kT) — 


Zy,1=2% cosh*—'J; 
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Therefore’ 


- 0 log Zy, 
N=-e— : me wg tee ST 
dJ 


C,/Nk=J? sech? J. 





and 


(b) 2XN strip.—Here the calculations are more tedious. They seem to be most easily performed 
by using the function A defined in Eq. 38. 


Zy,2= (cosh 3J+cosh J+A/2)""1(1+4/A). 
From the equation for the average energy given at the beginning of this section (neglecting terms 
of 0(1/N), the average internal spin energy per lattice point is 
e (3 sinh 3J/+sinh J+A’/2) 
2 (cosh 3J+cosh J+A/2)’ 
C,/2Nk=[10+10 cosh J cosh 3J—6 sinh 3J sinh J+(AA”’ — A’*) /4+(cosh 3J)(A’’+9A) /2 
+ (cosh J)(A’’+A)/2+(3 sinh 3J+sinh J) ]J*/2(cosh 3J+cosh J+A/2)?. 





E/2N= 


At low temperatures 


C,/2Nk~8J°e4 (14+-9e%! + 140 + - - -) 


and at high temperatures 
C,/2NkR~J?(3+9J?+ ---)/2. 


Wannier" has pointed out that the peak in the specific heat curve of plates N X ~ is of order log N. 
The author hopes to treat some of the other nearest neighbor problems mentioned at the beginning 
of this paper at some later date; and in conclusion he would like to thank Professor Onsager for his 
interesting discussions and for the numerous suggestions made throughout the course of this research. 


_' These results differ somewhat in form from those of Ising (see reference 4 or Bitter’s book, p. 145) because if one 
interprets his directed dipoles as spin orientations, the interaction between parallel spins would be zero. 
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Raman displacements, estimated intensities and depolarization factors are listed for di- 
tertiary-butyl ether, (CH3)s;COC(CHs)s, in the liquid state and the results are compared with 
previous data for di-normal-butyl and di-isobutyl ethers. A tentative assignment of the observed 
lines to frequencies of the di-tertiary-butyl ether molecule is made. 


INTRODUCTION 


RICKSON and Ashton’ have recently carried 

out the first successful synthesis of di- 
tertiary-butyl ether, (CH3);COC(CHs)s, and 
have kindly supplied the author with a sample 
of this material for use in a study of its Raman 
spectrum. The results of this investigation are 
given below. 


EXPERIMENTAL 


The spectrograph and the experimental tech- 
nique have been described previously.? The 
Raman spectrum was obtained in the liquid state, 
with excitation by Hg 4358A. The boiling point 
of the liquid was 106.5—107.0°C, d.° was 0.7658 
and mp” was 1.3949. Additional properties are 
listed elsewhere! and will not be repeated here. 
The sample was used as received from Professor 
Erickson, without redistillation. Only slight con- 
tinuous background was obtained, even upon 
long exposure. The film used was Agfa Super 
Plenachrome Press and development was in Agfa 
47 developer for six minutes at 18°C. The de- 
polarization factors were obtained by the method 
previously described,? except that now the 
Polaroid disk has been replaced by a Nicol prism, 
on account of the rather strong absorption of light 
by the Polaroid. It is estimated that this change 
approximately doubles the amount of trans- 
mitted light. 


1J. L. E. Erickson and W. H. Ashton, J. Am. Chem. Soc. 


63, 1769 (1941). 

? Forrest F. Cleveland, M. J. Murray, H. H. Haney and 
Julia Shackelford, J. Chem. Phys. 8, 153 (1940). 

3 Forrest F. Cleveland and M. J. Murray, J. Chem. 
Phys. 7, 396 (1939). 
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RESULTS 


The results are given in the columns at the 
left of Table I, in which Av represents the Raman 
displacement in cm, J represents the estimated 
intensity and p represents the depolarization fac- 
tor. Unusually sharp or broad lines are designated 
by the symbols s and 8, respectively. Data in 
regard to which there is some question are en- 
closed in brackets. The brace joins lines that were 
unresolved on the polarization spectrograms. The 
symbol P means that the line was somewhat 
polarized, but that an estimate of the numerical 
value of the depolarization factor could not be 
made. 


COMPARISON OF THE SPECTRA OF THE 
BuTtyL ETHERS 


The results previously obtained? for di-u-buty] 
and di-isobutyl ethers are included in Table I 
to facilitate comparison of these spectra with the 
present results for di-t-butyl ether. According to 
this table, it will be seen that the spectrum of 
di-t-butyl ether contains strong, polarized lines 
at 327, 567 and 701 and strong, depolarized lines 
at 298 and 902. All of these lines are absent or 
appear with low intensity in the spectra of the 
di-isobutyl and di-n-butyl ethers. Also, the ter- 
tiary ether has more lines in the 1200-1265 cm"! 
region and fewer lines in the 1400-1500 and 2700- 
3000 cm-! regions than the di-isobutyl and di-n- 
butyl ethers. 

Di-isobutyl] ether is distinguished by the group 
of three lines near 450 cm~ and the fairly strong 
lines at 810, 1249 and 1339, while di-n-butyl ether 
differs from the other two by the group of three 
lines near the 844 line and by the line at 1228 
cm—. The strong, polarized line at 840 appears in 
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ail three ethers. The C-H frequencies above 1400 
cm-! are quite similar for the di-isobutyl and 
di-n-butyl ethers, but are quite different in the 
case of di-f-butyl ether. It is interesting to note 
that the spectrum of tertiary-buty] alcohol‘ above 
1400 cm~ is practically the same as that ob- 
tained for the di-t-butyl ether, but that below 
1400 cm~! the two spectra, while similar in some 
respects, are by no means identical. 


ORIGIN OF THE D1-7-ButyL ETHER LINES 


It is not possible with the information at 
present available to make a complete assignment 
of the Raman lines to definite vibrations of the 
(CH3)sCOC(CHs3)3 molecule. A rough approach 
to this problem may, however, be made by con- 
sidering the frequencies to be expected from 
smaller groups within the molecule. The groups 
that make up the molecule may be considered to 
be a pyramidal CH; group, a pyramidal CC; 
group (neglecting the hydrogen atoms) and a 
non-linear R-O-R group, where R represents the 
rigid (CH3)3C group. Obviously, this is only a 
very crude approximation, but if it be accepted 
for the moment, one may then investigate what 
frequencies would be expected for such groups. 

The types of frequencies to be expected for 
such groups are listed in Table II. In the last 
column of this table are listed the values of the 
frequencies® for certain well-known molecules or 
ions which have the same symmetry as the groups 
under consideration. Using this table as a guide, 
one may then choose’ frequencies of the compli- 
cated molecule which correspond approximately 
to these group frequencies. The fundamentals so 
chosen are designated by the symbols va, », etc. 
in Table III. The frequencies va, », vr and vq have 
been selected as corresponding approximately to 
the frequencies 71, v2, v3 and v4 of the CC; group; 
ve, ve and », to the frequencies 7, v2 and v3 of 
the COC group; and », vi, vm and v;, »% to the 
frequencies v1, v2, v3 and v4 of the CH; group. 
Apparently the doubly degenerate frequency v4 
of the CH; group splits into the frequencies v; 


‘M. Magat, Tables Annuelles de Constantes et Données 
Numériques, Effet Raman (Hermann and Company, Paris, 
1937), Vol. 12, p. 32. 

_' Data taken from T. Y. Mu, Vibrational Spectra and 
Structure of Polyatomic Molecules (National University of 
Peking, Kun-ming, China, 1939), pp. 171, 203 and 238. 


OF BUTYL ETHER 723 
and », in the complicated molecule. There is some 
evidence for a similar effect in the corresponding 
frequency of the CC; group, since the frequency 
va, While not resolved into two components, is 
broader than usual on the spectrogram. 

Three lines were observed near the frequency 
chosen as », perhaps due to Fermi resonance 
interactions with some of the combination fre- 


TABLE I. Raman spectra of butyl ethers. 
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Di-TERTIARY Di-Iso D1-NoRMAL 
Av I p Av I p Av I p 
203} 
263 4 0.5 262 2 
276 2 
298 3b (0.9 299 1 294 2 
319 2 
327 3 04 
346 1 
369 3 362 2 
398 1 
424 3 418 3 
450 1 
467 1 
567 3. 04 
701 9s 0.1 
789 } 
810 2 
820 2 
841 7s 0.5 836 6 0.2 844 4 0.4 
882 1 
902 5 0.9 900 1 
930 5 
958 + 0.7 961 3 \ 
979 25 0.9 
994 1 0.9 |} 1001 1 
1026 2 0.9 1026 1 
1068 4 0.7 
1123 3 | 0.7 1120 3 ) 
1137 3f . 1137 3 | 0.8 
[1166] } 1168 3) g7| 1156 3 
1181 3/f ™ 
1210 4 0.9 
1226 4 09 1228 2 
1249 4 0.8 
1261 2 08 
1294 2 0.8 | 1300 5 0.8 
1339 4 0.8 
1392 i Ff 
1445 5 0.9 | 1448 7} 1439 6 ) 0.7 
1461 3 09| 1461 3} 08 | 1459 6f ; 
1478 2} 1484 3 0.7 
2709 2 0.3 | 2718 3 0.4 | 2732 2 0.5 
2793 3 0.4 | 2798 3 0.5 
2837 = 3 2849 3 
2868 2 0.3 | 2866 10 0.2 | 2864 10 } 0.2 
2876 610 
[2891] 1 
2924 10 0.2 | 2905 3 | 0.4 2912 8 ) 0.3 
2952 ¥ ‘in. | = 
2982 10 0.9 | 2963 4 0.8 2967 6 0.9 
2975 7 ™~ 














* See reference 2 in text. 
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TABLE II. Classification of the Raman frequencies. 








FRE- 











QUENCY SymM- SELEC- FREQUENCY 
SyM- DESIG- METRY TION? MOLECULE 
GROUP METRY NATION CLASS RULE or ION Av p 
CH; Cx V1 A, P CH;Cl® 2955 0.05 
v2 A\ Fr 1357 P 
V3 E D 3024 0.86 
V4 E D 1450 D 
CC; Cop V1 A, P IO3;— 779 0.43 
v2 A P 390 — 
V3 E D 826 0.67 
V4 E D 330 6/7 
COC Cop V1 A, P SO» 1151 P(gas) 
V2 A, r 525 Pi(liq.) 
V3 By D 1366 Dv(liq.) 








« P=polarized, D =depolarized. 

+ Only the frequencies which relate chiefly to the CH; group are listed 
here. There are two other frequencies v5 and v¢ in which the Cl atom is 
more directly involved. 


quencies vy+ve vet and va+r,. Fermi reso- 
nance also seems to exist between the funda- 
mental v, and some of the overtones 27;, 2v; and 
2v,.° The frequency v, is considered to be anal- 
ogous to the frequency vs in the methyl halides, 
in which the nucleus of the carbon atom and the 
center of gravity of the hydrogen nuclei in the 
methyl group are displaced in opposite directions, 
perpendicular to the symmetry axis of the methyl 
group. It is the motion of the methyl group which 
would correspond to rotation if the group were 
free. The frequency v¢ is doubly degenerate and 


6 It should be remarked that the 2924 line excited by 
Hg 4347A would fall at 2867 cm. However, it is believed 
that the intensity of the line observed at 2868 is too great 
to be accounted for by attributing it to this exciting line. 


CLEVELAND 


TABLE III. Tentative assignment of di-t-butyl 











ether frequencies. 
FREQUENCY FREQUENCY 
ASSIGN- ASSIGN- 
MENT Av I p MENT Av I p 
vye—Va 203 j Vetvp 1166] } 
Va 298 3b 0.9 ve+va 1210 4 0.9 
vr | 
Vb 327 3 O04 w+ { | 1226 4 0.9 
Va— Vb 369 3 vat ve } | 1261 2 0.8 
Va-—Va 424 } Vi 1392 if 
Ve 567 3 O04 v; 1445 5 0.9 
Va 701 9s 0.1 Vk 1461 3 @3S 
Ve 841 7s 0.5 2»; | 2709 2 0.3 
Dy. | 
wy 902 5 09 WwW 12868 2 0.3 
Vat va | f 994 1 0.9 VI | 2924 10 0.2 
s+ 4 
2 0.9 Vm 2982 10 0.9 


Vat vp J { 1026 


it is not clear whether the frequencies observed 
at 994 and 1026 result from the splitting of », 
into two frequencies in the large molecule or 
whether the two lines result from Fermi reso- 
nance interaction between vy, and one of the 
combination frequencies va+v_ Or vat. 

The four remaining frequencies have only small 
intensities and may be ascribed to combination 
frequencies, as indicated in Table III. Some of 
them may, however, be fundamentals, since for 
the di-t-butyl ether molecule there should be 
(3N—6) =75 fundamentals and only thirteen (or 
fourteen, if both the frequencies 994 and 1026 
are regarded as fundamentals) have been chosen 
above. It should again be emphasized that the 
above assignments are only tentative because of 
the very rough approximations made in arriving 
at them. 
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Letters to the Editor 








HIS section will accept reports of new work, provided 

these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 





Isotopic Exchange Between Carbon Tetra- 
bromide and Bromine 


J. H. HopGeEs aAnp A. S. MICELI 
University of Michigan, Ann Arbor, Michigan 
April 14, 1941 


EIGHED quantities of carbon tetrabromide and 
of bromine containing radioactive isotopes were 
allowed to react in Pyrex vessels of sufficient size to permit 
complete volatilization of the reactants. The radioactive 
bromine was obtained by exposing bromoform to neutrons 
from the cyclotron; it was allowed to age for periods of 
about 24 hr. to allow for decay of the two isotopes of 
shorter life.! The activity of the samples was measured with 
a Lauritsen electroscope. The rate of exchange of radio- 
active bromine could be determined conveniently over the 
temperature range 170° to 200°C. Molar concentrations of 
carbon tetrabromide and bromine were varied from a ratio 
of 10:1 to a ratio of 1:10. The reaction appears to be 
homogeneous, and the rate appears to. be directly pro- 
portional to the concentration of the carbon tetrabromide 
and to the square root of the concentration of the bromine. 
This indicates a mechanism involving bromine atoms. If it 
is assumed that exchange occurs between a carbon tetra- 
bromide molecule and a bromine atom, the energy of activa- 
tion of this bimolecular process is zero. This is obtained by 
plotting log [(measured exchange rate)/(CBr,)(Br)] 
against 1/T. The value of (Br) is calculated from the equi- 
librium constant of the reaction Br2(g)—2Br(g). 

The same procedures were repeated using Pyrex capsules 
of small capacity in order to maintain the reactants in the 
liquid phase. The exchange rate was measured over the 
temperature range 107° to 160°C. The mechanism of the 
reaction seems to be the same as that in the gaseous phase. 
In this case, however, there is an energy of activation of 
about 3 kcal. The values of the activation energies for the 
reaction in the gaseous and liquid phases may be somewhat 
in error, first because of the extrapolation of the equilibrium 
constant expression to the rather low temperatures of our 
experiments, and second because of the inherent limits of 
accuracy in our experimental technique. An additional 
small error is involved in the use of this same equation in 
calculating the degree of dissociation in the liquid phase. 


JOURNAL OF CHEMICAL 


725 





PHYSICS VOLUME 9 


No parallelism can be detected between the intensity of 
the radioactivity of the bromine sample and the value of 
the specific reaction rate calculated for the runs. Further- 
more, no exchange occurs at room temperature. These facts 
seem to indicate a thermal reaction and not one initiated by 
the energy of the nuclear decomposition. 

Calculation of theoretical rates on the basis of the 
collision theory, assuming a bimolecular reaction between 
carbon tetrabromide molecules and bromine atoms gives 
values in both the gaseous and liquid phases which are too 
large by a factor of between 104 and 105. 

1 Willard, J. Am. Chem. Soc. 62, 256 (1940). 


2 Liberatore and Wiig, J. Chem. Phys. 8, 165, 349 (1940); Libby, ibid. 
8, 348 (1940). 





Infra-Red Absorption Spectrum of Acetone 


DONNA PRICE 


Mallinckrodt Chemical Laboratory, Harvard University, 
Cambridge, Massachusetts 


July 29, 1941 


HE infra-red spectrum, from 3—20y, of gaseous acetone 

was obtained with a spectrometer! using fluorite, 
rocksalt, and KBr prisms. The acetone used was Mallinc- 
krodt C.P. which had been refluxed and redistilled in this 
laboratory. 

The spectra obtained for the three prisms are shown in 
Fig. 1. There is, though this is not evident in the figure, a 
shoulder on the low frequency side of the 1363 band. The 
frequencies and their qualitative intensities are given in 




















PERCENT TRANSMISSION 
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FREQUENCY 
Fic. 1. The infra-red absorption spectrum of acetone gas. The trans- 

mission curves were taken in a 30-cm glass cell using NaCl windows 


with gas at the following pressures: a, 5 mm; b, 1.1 mm; c, 16 mm; 
d, 200 mm. 
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TABLE I. The infra-red spectrum of acetone gas. 











FREQUENCY 
IN cm™! INTENSITY FREQUENCY INTENSITY 
547 medium 1964 weak 
774 medium 2109 medium 
896 strong 2234 weak 
1202 740 cm~ 2418 weak 
1363 1007 cm=! 2554 weak 
1430 strong 2984 strong 
1725 1552 cm! 











Table I. By use of Wells’ method,? the intensities of the 
three strongest bands were computed for =100 mm. 
These results are the numerical values for the 1202, 1363, 
and 1725 bands given in Table I. It was assumed that ace- 
tone was at least as complicated as propylene, and hence 
no dilution with an inert gas was necessary. Moreover, the 
region 1.1-5 mm of the curve, In (Do/D) vs. p, was assumed 
practically a straight line. Wells’ curves seem to justify 
this, but since only two points were taken and these deter- 
mined by graphical integration, the values are given only 
as semi-quantitative approximations. Corrections were 
made for the blanks, but the stray effect was considered 
negligible in this region. 

No other data are available for the vapor, but Coblentz 
obtained the infra-red absorption spectrum for liquid ace- 
tone, and Edsall and Wilson? have obtained the most 
recent Raman spectrum. 

1H. Gershinowitz and E. B. Wilson, Jr., J. Chem. Phys. 6, 197 
(1938); E. B. Wilson, Jr., and A. J. Wells, J. Chem. Phys. 9, 319 (1941). 

2A. J. Wells, thesis, Harvard University, Cambridge, Massachusetts 


(1941); E. B. Wilson and A. J. Wells, J. Chem. Phys. 9, 659 (1941). 
3 J. T. Edsall and E. B. Wilson, Jr., J. Chem. Phys. 6, 124 (1938). 





The Exchange Reaction Between Gaseous and 
Combined Nitrogen 


T. H. Norris, S. RUBEN, AND M. D. KAMEN 


Chemical Laboratory and Radiation Laboratory, University of California: 
Berkeley, California 


August 13, 1941 


SING radioactive nitrogen (N') as an indicator, 
Nishina, et a/.! have reported that Ne gas undergoes 
rapid exchange at room temperature with dilute (6.1-0.2M) 
solutions of NaNO3;, NaNOo, NH2OH-HCI, etc.Moreover, 
according to these authors, the rate of exchange decreased 
with increasing concentration of nitrate or nitrite. Since 
these results are both unexpected and at variance with 
other investigations,?:? it has been considered desirable to 
repeat the experiments. 

Charcoal was bombarded in a gas-tight chamber (con- 
taining 100 cc tank N2) with 30u of 8-Mev deuterons. 
Under these conditions ~30 percent of radioactive nitrogen 
formed! is present in gaseous state. An appreciable fraction 
of the gaseous activity is present as combined nitrogen? 
(formed by the recoiling N™ as a result of deuteron impact). 
Since the presence of combined nitrogen “impurity” can be 
construed as evidence for exchange, it is important to 
remove it from the labeled Ne. To this end the gaseous 
activity was pumped out of the target chamber and mixed 
with carrier (inactive) HCN and NOs. The mixture was 


slowly pumped through a heated combustion tube con- 
taining CuO and was then passed through three glass 
spirals immersed in liquid air. By means of two Tépler 
pumps, one at each end of the train, the gas was forced 
back and forth several times through the CuO and cold 
traps. 

After this treatment two portions of ~50 cc of the N»* 
were vigorously shaken for ten minutes with 10 cc NaNO; 
and 10 cc NaNO: solution in separate (previously evacu- 
ated) vessels. The solutions were then gently boiled to 
drive off dissolved Nz and were analyzed for radioactivity 
using methods described elsewhere.’ The results are 
summarized in Table I. 


TABLE I. 











SUBSTANCE CONCENTRATION RADIOACTIVITY? 
TESTED (MoveEs/LITER) (Counts/ MINUTE) 
1.2 <4 
NaNOsz 
ok <4 
1.3 <3 
NaNOz 
| <5 








2 All figures are corrected for decay of N" to the time of counting the 
first sample and are therefore comparable. The N2* had an activity of 
105 counts/minute. 


It is apparent from Table I that <0.01 percent of the 
N.* has exchanged with NO3;" or NO2-. This is in agreement 
with the results of Burris and Miller? on Ne fixation by 
Azobacter using heavy (N') nitrogen gas. 


, — limori, Kubo, and Nakayama, J. Chem. Phys. 9, 571 
1941). 

2 Ruben, Hassid, and Kamen, Science 91, 578 (1940). 

3 Burris and Miller, Science 93, 114 (1941). 

4 The remainder can be removed by completely burning the charcoal 
to carbon dioxide. This was not undertaken due to the short half-life of 
N' (10.5 minutes). 

5 Ruben, Kamen, and Hassid, J. Am. Chem. Soc. 62, 3443 (1940). 





Measurement of Tension in Solutions 


ARTHUR F. Scott AND G. MARSHALL PoUND 
Department of Chemistry, Reed College, Portland, Oregon 
July 18, 1941 


LTHOUGH many properties of solutions have been 
studied, there appears to be no record of any sys- 
tematic study of the maximum hydrostatic tension which 
solutions can withstand under varying conditions of tem- 
perature and concentration. With the publication of Vin- 
cent’s ingenious method for measuring this tension such a 
study seemed possible and we have accordingly undertaken 
it in the hope of obtaining information of value to.the prob- 
lem of the structure of solutions. This letter is a preliminary 
report on the first experiments of this project. 

Vincent’s apparatus can be described briefly. It consists 
of a metallic bellows, closed by a needle valve, which con- 
tains the liquid under test. Inward leakage of air into the 
bellows is eliminated by surrounding the bellows and valve 
by the same liquid. With the bellows filled and the valve 
closed a gradually increasing load is applied to the bellows 
until the liquid breaks. The hydrostatic tension existing 
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TABLE I. 
MAXIMUM 
PERCENT TENSION 
No. Ca(NOs)2 TEMP atmos. 
1 6 20 1.59 
6 0 20 1.53 
1 0 23.5 1.43 
4 21 20 0.82 
3 34 21 0.87 
2 51 21 0.58 
2 51 23 0.58 
2 51 25 0.45 
5 66 19 0.28 








in the liquid prior to the breaking is then calculated from 
the load and the known effective area of the bellows. Our 
apparatus was constructed last winter from designs sent 
us by Vincent. When Vincent’s complete paper was pub- 
lished! recently it was discovered that our procedure and 
experience with the apparatus were much the same as he 
described. Two exceptions, however, should be noted. First, 
in all but the first experiment we used an arrangement 
whereby the bellows could be filled under vacuum. Second, 
we could not detect bubbles escaping from the apparatus 
after breaking as Vincent observed with three organic 
liquids and we found no tendency for the maximum tension 
to increase on repeated breaks. 

For our first measurements we chose aqueous solutions 
of calcium nitrate because they are known? to be stable to 
a high degree of supersaturation. The results of our meas- 
urements with these solutions and pure water are sum- 
marized in Table I. 

In this table the number gives the chronological order in 
which the experiments were carried out. The fact that the 
first and last experiments yielded sensibly the same value 
for the tension of pure water and the fact that the results 
do not vary in the order in which they were obtained would 
seem to establish that no permanent distortion of the bel- 
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lows took place. The recorded values for the tension are 
the means of at least four determinations of the force 
necessary to break the liquid in the bellows. These tension 
values were calculated as described by Vincent and might 
be uniformly in error by possibly ten percent because of an 
uncertainty in the effective cross section of the bellows. 

The most significant fact revealed by these fragmentary 
data is that the maximum tension of solutions of calcium 
nitrate diminishes with increasing concentration. This 
would seem to indicate that the solute weakens rather than 
strengthens the bonds between the solvent molecules. It 
is of interest to note also that in a plot of the maximum 
tension values against values of the apparent molal volume 
of the solute at corresponding concentrations the plotted 
points fall closely on a straight line. Assuming this linear 
relationship and extrapolating to the hypothetical solution 
state at which the tension which the solution can with- 
stand is zero, the value for the apparent molal volume at 
this hypothetical limit of the solution state is found to be 
66 ml. This figure is of the order of magnitude of the volume 
of anhydrous calcium nitrate, measurements of which 
range from 67.5 ml to 73.4 ml. Finally, it may be pointed 
out that the maximum size of bubble which could with- 
stand the observed maximum tension will increase with 
increasing concentration because in the case of this solu- 
tion the surface tension increases while the tension and 
vapor pressure decrease with increasing concentration. This 
is noteworthy because Vincent calculated from his meas- 
urements that the radius of the maximum bubble for alco- 
hol, ether, and mineral oil was approximately constant and 
of the order of 2.0 10-5 cm. 

The authors are greatly indebted to Dr. Fred D. Ayres 
for invaluable help in building the apparatus and conduct- 
ing the experiments. 


1R.S. Vincent, Proc. Phys. Soc. 53, 126 (1941). 
2A. F. Scott and G. L. Bridger, J. Phys. Chem. 40, 461 (1936). 





